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INTRODUCTORY REMARKS 


F. R. Eijirich 
Institute of Polymer Research, Polytechnic Institute of Brooklyn, Brooklyn, N. Y. 


The performance of high polymers in all forms of usage and their be- 
havior during measurements is intimately linked to their structure. Solu- 
tions of long-chain molecules, in particular, show very unusual solubility 
patterns and a readiness to exhibit non-Newtonian viscosity. The papers 
presented in this monograph deal with some new developments in our views 
of these phenomena. 

The first advances in an understanding of polymer solutions were made 
by P. Flory and M. Huggins, in conjunction with the model of the random 
flight chain molecule treated by W. Kuhn and H. Eyring, by considering 
the configurations of chain molecules in solution and calculating the cor- 
responding entropy of mixing. Notwithstanding the great progress rep- 
resented by this approach, the limits of some of the basic assumptions and 
of the agreement with the experiments remained as a stimulus for further 
work. Various refinements were required, especially with respect to the 
lattice model, to the specific rotational characteristics of the chains as 
influencing coil dimensions, and to the energetic interactions of the polymer 
molecules with the solvent and one another. 

Much progress has been made subsequently by Flory and his school by 
the ingenious analysis of long- and short-range polymer interaction in 
terms of the expansion factor a and of the ideal random flight dimensions 
in theta (@) solvents. A survey of the most recent advances is given in 
the paper by Peterlin with special reference to chain structure and to the 
effect and the nature of the coil expansion. The paper by Wall and Mazur, 
then, makes an important contribution to the problem of coil dimensions 
by computer analysis of the configurations as a function of lattice type 
and of bond angle of chains with three levels of energy interaction. The 
influence of the presence of other polymer molecules, that is, the concen- 
tration effect on coil dimensions, is then discussed by Fixman on the basis 
of a Gaussian model for the distribution of chain segments. The paper by 
Kurata makes an important contribution to the application of the lattice 
theory to polymer solutions. While recognizing the general validity of 
Flory’s derivation of the entropy of mixing, the theory is extended to take 
account of the inevitable close connection of coil segments, due to their 
chemical linkages, and the concomitant solvent clustering. Extensive use 
is made of experimental data from the @ state to derive functional forms 
and values for the interaction constants. 

It might be added parenthetically that the formulation of the 6 state is 
not without conceptual difficulties because of the fact that the “ideality” 
of @ solutions (as at the Boyle point) is operational and not real. It is 
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true that the validity of van’t Hoff’s limiting law shows the existence of 
ideal additivity, or of proportionality, between chemical potential and log 
mole fraction, and therefore indicates an ideal entropy. However, as in 
other regular solutions, the AH term is not zero but positive, so that ran- 
dom mixing at temperatures close to precipitation in any but the most 
dilute solutions is most unlikely. In other words, with lowering tempera- 
tures the distribution of the solute must pass a state of clustering on the 
way from complete randomness to the critical point. Although the tem- 
perature range of such clustering is probably narrow, as Debye’ has re- 
cently found, with increasing molecular weight this range must move up- 
ward toward the @ temperature and eventually extend above it. In the 
latter cases, the absence of bias in the segment-solvent or solute distribu- 
tion at the @ point can be only an apparent one. More likely, the “‘ideal- 
ity” will be produced by some state of order that also leads to a cancella- 
tion of the factors contributing to the virial coefficients; to the extent that 
this happens deductions from 6-point data may have to be revised. Re- 
ported discrepancies in the entropy and interaction parameters from 
different methods or between calculated and experimental values, or ab- 
normally large values of the ratio Mw : My may be due to this inadequacy 
of our picture of the polymer distribution. 

A most unexpected behavior of polymer chains that has recently come 
to light and that may well be related to the state of coiling at the critical 
point is the crystallization in thin platelets by way of uniformly folded 
chains. There are many puzzling features to this phenomenon, not the 
least being its widespread occurrence and the apparent disregard for ac- 
quiring a large surface energy. Fischer’s paper presents a thermodynamic 
theory that has found wide support in Europe. A different, kinetic, hy- 
pothesis discussed by F. P. Price is to be published elsewhere.?:* 

The dimensions, density, and rigidity of the polymer coils express them- 
selves in their contributions to the solution viscosity. The second part of 
Peterlin’s paper reviews the recent status of our knowledge in this area. 
In particular he shows that we have now a much better understanding of 
the gradient dependence of the viscosity of dilute solutions based on intra- 
molecular hydrodynamic interaction and coil deformation, but a generally 
satisfactory theory has not been proposed. Zimm, in a brief report, out- 
lines a new approach that seems promising. 

Generally, the papers in this monograph are concerned with the discon- 
tinuous nature of polymer solutions and the interaction of chains with the 
surrounding medium. With increasing concentration the separation of 
the particles disappears more and more, until at high concentrations or in 
the melt a practically continuous medium occurs. At this stage it be- 
comes fruitful to inquire into the general properties of media that exhibit 
the characteristics of polymer systems, especially their viscoelasticity. 
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The mathematical-mechanical analysis of such fluids has made great prog- 
ress under the impetus of work by A. E. Green, R. 8. Rivlin, J. E. Erickson, 
and others, a development that has not received as wide an attention as 
it deserves. The concluding paper by Coleman and Noll affords an excel- 
lent acquaintance with the basic approach and results in this field. 
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SHORT- AND LONG-RANGE INTERACTION IN THE 
ISOLATED MACROMOLECULE* 


Anton Peterlint 
Department of Chemistry, Wayne State University, Detroit, Mich.t 


Preface 


Extensive hydrodynamic and optical studies on dilute solutions have re- 
vealed a great many data on size and other physical properties of isolated 
linear macromolecule. The dependence on molecular weight, temperature, 
solvent, and rate of shear has shown the inadequacy of the original concept 
of the Gaussian coil, which did not take into proper account either the fine 
structure of the chain or the steric, energetic, hydrodynamic, and optical 
interactions between distant sections of the chain. In the last 10 years, 
therefore, more refined models were introduced in order to obtain agreement 
with experiments. Their success makes worthwhile a short survey with 
emphasis on the remaining gaps that it will be necessary to cover in the 
future. 

The actual shape of statistically coiled linear macromolecule depends on 
valency angle between subsequent monomers, steric and energetic hindrance 
of rotation around the chemical bond (short-range interactions), on the 
volume requirement of the monomers, and on energetic interaction between 
any two monomers and between the monomer and the surrounding medium 
(long-range interactions). The short-range interactions determine not 
only the length of the statistically independent segment, its optical and 
hydrodynamic properties, but also the resistance of the chain against the 
rate of shape change of the coil. The long-range interactions, however, do 
not appreciably change the character of the segment but influence the over- 
all coil dimensions by excluding sterically or energetically unfavored con- 
figurations. As a consequence, the root-mean-square (r.m.s) gyration 
radius does not increase as the square root of molecular weight. Propor- 
tionality with M"” is found only exceptionally, that is, in ‘‘ideal’’ solvent. 
The size and shape of the coil in turn influence the hydrodynamic and opti- 
cal interaction. 

The non-Gaussian character of actual coil, the anisotropy of hydrody- 
namic and optical interaction, and the resistance of the coil against rapid 
shape changes influence the gradient dependence of limiting viscosity and 
streaming birefringence. In spite of the fact that the basic character of 


* The investigation reported in this paper was supported in part by a research 
grant from SKNE, Beograd, Yugoslavia. a ‘ : 


_t On leave of absence from the Department of Physics, University Ljubljana, 
Ljubljana, Yugoslavia. 


{ Present address: Physikalisches Institut, Technische Hochschule, Minchen 
Germany. 
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the single contributions is already known, no one has yet succeeded in con- 
sidering them all together in the full gradient range. 


Introduction 

In dilute solution a linear unbranched macromolecule containing P mono- 
mers in the chain assumes a coillike shape and steadily changes its configura- 
tion so that in any experiment mean values over all possible configurations 
are observed. From investigation of monomers one knows, at least in 
principle, their volume v , link length bo , valency angle 180° — a, electric 
charge éy , dipole moment yo , optical polarizability tensor ao , hydrodynamic 
resistance coefficient Aj and soon. By taking into account all kinds of pos- 
sible mutual interactions of monomers and averaging over all permitted 
configurations under the eventual influence of an external field one expects 
to be able to explain the effects observable with dilute polymer solutions 
and their dependence on molecular weight, temperature, solvent, and sol- 
vent composition (ionic strength, pH, mixture of solvents). The interac- 
tions may be due to the geometry of monomers, their space requirement, 
to attractive and repulsive forces between two monomers and between a 
monomer and a molecule of the solvent, and to hydrodynamic forces and 
electric fields originating from free charges and electric dipoles. 

The effects concerned are primarily X-ray and light scattering, dielectric 
polarization, sedimentation in gravitational and inhomogeneous electric 
field, viscosity, electrical and streaming birefringence. What matters is 
always the excess effect of solution over that of solvent divided by con- 
centration and extrapolated to infinite dilution. These speczfic or intrinsic 
effects all depend on the properties of the monomer in given solvent and on 
their mutual geometric arrangement in the chain and at the special condi- 
tions of experiment, that is, on shape and size of the macromolecular coil, 
which, for example, determine the extent of optical and hydrodynamic inter- 
action and the resulting modification of observable intrinsic effects. This 
fact stresses the importance of coil geometry, which therefore must be con- 


sidered first. In order not to complicate the situation too much let us start 
- with solution at rest without any external field.” 


A very rough but useful model for such a randomly coiled macromolecule 
is a chain with Z statistically independent segments of equal length A 
joined in such a way that the orientation of any segment is, at least in first 
approximation, independent on that of its neighbors. Within a polymer 


series and under unchanged experimental conditions Z is proportional to 


P, the proportionality factor ¢ being the number of monomers in a segment. 
Some properties, such as the volume and the charge of the segment, are 
strictly equal to ¢ times the corresponding values of monomer. In general, 
however, they are influenced by the special geometry of monomer arrange- 


ment in the segment and on their interaction so that they are es than 


¢ times the values of monomer. 
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Chain configurations and their probability of occurrence are determined 
by forces among all chain atoms and those of solvent. In this generality, 
the problem is unduly complicated and must be simplified properly by ap- 
plying generally accepted geometric and energetic formulations of above 
mentioned forces. Thus, for example, forces between first neighbors result 
in special valency angles and interatom distances and in a more or less 
fixed monomer structure that is usually known from crystallographic and 
X-ray investigation. The monomer unit may be rigid (cellulose) or flexible 
due to single bonds with more or less free rotation around it (for example, 
vinyl polymers, polyesters, and polyamides). The same forces also deter- 
mine the possible geometric arrangement of subsequent monomer units and 
hence the geometrical properties of statistical segment. 

With monomer units symmetric in chain direction, as for example in 
polyethylene —(CH;-CH:)—, and with those that, although asymmetric, 
require special orientation for regular valency bond between subsequent 
monomers, as for example in polysiloxanes —(SiRiR2-O)— and cellulose 
—(C.H»0;)—, the atom or atom group sequence in unbranched polymer 
chain is absolutely determined and regular. With most vinyl polymers 
such as —(CH,-CHR)—, however, the longitudinal assymmetry gives rise 
to four possible arrangements (head-head, head-tail, tail-head, tail-tail), 
which in polymer chain may occur either in high order (positionally isotactic 
-ht-ht-ht- or syndiotactic -hh-tt-hh-tt- polymer) or distributed at random 
(positionally atactic polymer -ht-hh-th-tt-hh-), depending on polymeriza- 
tion conditions. Moreover with monomers exhibiting two-mirrored d- and 
l-isomers (most vinyl polymers are built from such monomer units 
—CH,:CHR—), the resulting polymer may contain both, either in well- 
ordered arrangement or distributed at random, or only one. In the case of 
positional order (positional tacticity) this monomer isomerism leads to 
sterically isotactic (all monomers are either d- or /-isomers), syndiotactic 
(d-l-d-l sequence) or atactic (random distribution of both isomers) poly- 
mers. A similar situation arises with monomers exhibiting czs- and trans- 
isomery, as for example in polybutadiene and polyisoprene, leading to cis- 
and trans-tacticity, syndiotacticity, and atacticity respectively. 

The above-mentioned polymers with regular chain order are positionally 
tactic, but may show different degrees of sterical tacticity. Polyethylene, 
of course, exhibits complete positional and sterical order (isotacticity). 

The spatial conformation that results with given fine structure of the 
chain, however, still depends on many other factors that may be either 
purely geometrical, due to volume requirement of single constituents and 
valency angles between adjacent groups, or energetic, due to special inter- 
action between any two monomers’ units or between such a unit and the 
solvent. It was soon found that some factors affect chiefly the single seg- 
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ment (short-range interactions) and that others affect the mutual arrange- 
ment of segments in the random coil (long-range interactions). A special 
situation arises with electric charges that influence both. This last effect 
will be only mentioned at the proper place without going into any further 
detail, because this would require an extended review by itself. 


Short-range Interactions 


For sake of simplicity let us first consider the polyethylene chain. One 
can take CH, as monomer, by = 1.53 A, valency angle 109.5° = 180° — « 
(tetrahedral bond with cosa = 14 With completely free rotation around 
the valency bond the mean square ead distance of such a chain with 
P = N + 1 monomers reads* 


i =) = (Xb)) = VE aby) 


= bo [N(1 — cos’ a) — cos a(1 — cos” a)]/(1 — cos a)? 


(1) 


with b; the link joining the jth and the 7 + 1st monomer. But there is no 
experimental possibility for direct determination of h. Experiments yield 


- only the mean square gyration radius* 


rR? 


3 » (tin)/2(N + 1)" = 20 (N + 1 — 8)h'/(N + 1)” 
be {Lt cose N(N +2) _ COs a | N 
1—cosa 6(N + 1) ~ (1 — cos a)! N+1 


ENN He) in Note aay 2 cont” =|} 
(Ne bola Ak eM Pad (WV + 1)? 


rj, vector joining the jth and kth monomer, h, r.m.s. end-to-end distance 
of a chain with s links or s + 1 monomers. With increasing N both ex- 
pressions may be simplified to 


h? = Nbo'(1 + cos a)/(1 — cosa) = ND’ = KM 
R? = Nbo’(1 + cos a) /6(1 — cos a) = Nb’/6 = h’/6 = KM /6 


(2) 


(3) 


The effective link length 
b = bol(1 + cos a)/(1 — cos a)]'” = boctg(a/2) (4) 


is always larger than the true link length b). In the case of tetrahedral 
bond angle one has b = 2”"byp = 1.414bo . 

An additional useful quantity for characterizing the chain and especially 
the short-range interactions on segment length is the persistence length 


a = b/(1 — cosa) (5) 
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5,5a 


introduced by Kratky and Porod. It equals the mean length of projec- 
tion of infinitely long-chain on the direction of first link. 
The model with Z independent segments of length A yields 


hz = ZA’ 

R; =ZA/6 
for mean square of end-to-end distance and gyration radius respectively. 
In order to obtain Z and A from M or P or N and R or h one wants addi- 
tional relationship between Z, A, and an observable quantity. Kuhn and 
Kuhn! introduced the rather arbitrary condition that the model and the 
molecule have the same extended length 

L = ZmAm = Nbo cos (a/2) (7) 


Combination of Equations 3, 6, and 7 then yields Kuhn’s preferential quan- 
tities 


(6) 


A, = W/L = iy Se = b/sin (a/2) = 2a cos (a/2) 
7. RN TRICL = Gomialie Nose are) (8) 
tm = N/Z = 2/(1 — cosa) = 1/sin’ (a/2) 


One must not forget, however, that Equation 7 is rather unrealistic, yield- 
ing too high values for A,,. The actual statistically independent segment 
is shorter: A < An. 

According to Equation 3 both quantities h and # in the limit of large V 
are proportional to the square root of molecular weight. With small N, 
however, the increase of R’ and h’ with M is more than linear. The cor- 
responding effective bond length by starts with bo at N = 1 and approaches 
asymptotically the limiting value b (FIGURE 1), the approach being the 
slower the closer cos a to unity. The ratio 6R’/h’ also appreciably differs 
from unity at low NV. It starts with 1.5 at N = 1, rapidly decreases to a 
minimum and then rather slowly approaches from below the limiting value 
1 (riguRE 2). In that which follows N will be always so large that all these 
effects, characteristic for short chain, that is, for chains having small Z, 
may be neglected and the asymptotic formulas Equation 3 used instead of 
Equations 1 and 2. In such a case the distribution function of interchain 
distances, with exception of those between closest neighbors, is Gaussian. 

Experiments, that is, angular dependence of light scattering,’ combina- 
tion of intrinsic viscosity and sedimentation or diffusion constant® yield 
R’. Proportionality with molecular weight is found in zdeal solvent char- 
acterized by vanishing of second virial coefficient, that is, at the so-called 
Flory’s 9-point. The corresponding gyration radii Re are markedly smaller 
than in good solvents, but they are definitely higher than predicted by 
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by= h/N'/? cos &= 0.95 
4 09 
L 0.8 
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Figure 1. Initial increase of effective link length by = h/N1/? with N from bp at 
N = 1 till b = & [A + cos a)/( — cosa) at N = ~. 


4 
” 
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B. 2. Initial changes of the ratio 6R2/h? for different valency angles 180° — 
Bm a. eee for the shortest chains—the ratio starts with 1.5 at one link with two 
__ beads—and very small cos a the square of gyration radius is slightly less than h?/6. 
The lowest curve corresponds to a wormlike chain®:** with Z = L/2a.. 
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Equation 8. An experimental effective link length be can be calculated by 


applying Equation 3 


heise ONY 


(9) 


The ratio be/b is always larger than unity (TABLE 1). It is 1.2 for poly- 
ethylene and 2.1 for polymethylmethacylate. 


TABLE 1 


Cor, Extension Dus To HINpERED ROTATION AS MEASURED BY 
bo/b = [(1 + (cos ¢))/( — (cos g))|? FOR DIFFERENT 


POLYMERS AT 9-POINT 


be/b Reference 
Polyvinylbromide te ae 1.9 9 
Br 
Polyacrylic acid —CH,—CH— 1.8 10 
COOH 
Polyvinylacetate —CH,—CH— 2.3 11 
COOCH; 
Polystyrene —CH.—CH— 2.22.4 | 12,13 
C,H; 
Poly-4-vinyl pyridine —CH,—CH— 2.2 14 
eH sN 
Polyisobutylene —CH.—C(CHs3)2— 2.3 12 
1.7-1.8 | 13, 15, 16 
Polymethylmethacrylate —CH,—C(CH;)— 2.0-2.2 | 10,17, 18 
COOCH; 1.68-1.89 | 19 
Polyethylmethacrylate —CH,.—C(CH;)— 1.9 20 
COOC2H; 
Polybutylmethacrylate —CH,—C(CH;)— 2.1 21 
COOC,Hy 
Poly-2-ethylbutylmeth- —CH.—C(CH;)— 2.3 22 
acrylate 
OOC,Hs(C2Hs) 
Polyhexylmethacrylate —CH.—C(CH;)— 2.4 23 
COOC.«His 
Polyoctylmethacrylate —CH.—C(CH;)— 2.3 24 
OOCsHi, 
Polylaurylmethacrylate —CH.—C(CHs) 2.8 24a 
COOCi2Hes 
Polydimethylsiloxane —Si(CH;).—O— 162 25 
Poly-cis-isoprene —CH,—CH=C(CH;)—CH.— 15 26 
Poly-trans-isoprene —CH,—CH=—C(CH;)—CH.— 1.3 26 
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In order to explain this additional chain extension one has to consider 
restriction of rotation around the chemical bond caused by geometric and 
energetic interaction of close-by monomers. In the limit of sufficiently 
high N, restricted rotation influence, calculable by the same scheme as dem- 
onstrated in Equation 1, leads to the following modification” of Equation 4 


p? = p2 it cos @ 1 + (cos ¢) 
"eos ot (cos ¢) 


(10) 


with ¢ the angle between the planes through the jth link and the preceding 
and subsequent link respectively. With pure trans (cis) configuration one 
has ¢ = 0 (180°). Restricted rotation therefore may either increase or 
decrease the effective link length and hence R and h according to which con- 
figurations, cvs or trans, are favored by neighbor interaction. Bulky side 
groups have high space requirement and hence prevent chain configurations 
at which the groups would partially overlap, and these are mainly configura- 
tions yielding a shorter segment. Very efficient are phenyl and methyl 
groups. Rather long aliphatic side chains do not matter very much since 
their bulk can be easily accommodated far enough from the main chain. 
There may be a special temperature effect found by Schulz and Kirste” 
polymethylmethacrylate. The ratio be/b steadily rises from 1.68 at 

= —40° C. (2-methylpentanon-4 and methylisovaleriate as solvents) till 
1.89 at +60° C. (isoamylacetate). The authors propose as explanation 
that, with increasing temperature, the side groups have more motional free- 
dom and hence increased space requirements so that the rotational restric- 
tions to the main chain are greater than at lower temperature. 

Forces between adjacent monomers and among neighboring groups, char- 
acterized by a distance-dependent interaction energy u(r), may have vari- 
ous effects, all of them primarily showing a temperature dependence of the 
type exp(—u/kT). With exception of forces between free electrical 
charges in a solvent of low ionic strength they are typical short-range forces 
decaying with a very high power of mutual distance so that-one has to con- 
sider the interaction of but a small number of neighbors.”'”* For the most 
part merely the influence of first neighbor is taken into account. The sol- 
vent is considered only in so far as it changes the coefficients in the inter- 
action energy. In the simple polyethylene chain one has to do only with 
the repulsion of protons favoring trans configuration and hence increasing 
the effective link length. With increasing temperature these limitations 
get less severe, the corresponding decrease of the average (cos ¢) leading 
to a reduction of b and hence of coil diameter and gyration radius. 

Side chains with active groups complicate the situation. An extreme 
case is DNA and other polynucleotides in which the strong hydrogen bonds 
between the sugars lead to double or even triple stranded helical configura- 
tions” that enormously increase b and persistence length a (~300 A). sade 
Such a helix is an extremely rigid structure as compared with single chain, 
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even when the latter exhibits high rotational restriction. As a consequence 
the deviations in coil structure from the limit (Equations 3 and 10) are 
not completely negligible as may be seen from angular dependence of light 
scattering (FIGURE 3). Hydrogen bonds and electrostatic interaction of 
charges in this kind of molecules are very sensitive to solvent composition, 
that is, ionic strength, pH, presence of groups competing with the intra- 
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helical hydrogen bonds (for example, urea), and one actually observes 
“melting” and a net influence of these factors on “melting” temperature 
at which the helices transform in normal chains with high flexibility." A 
limited amount of helical arrangement is also found with isotactic poly- 
styrene where, according to Birstein and Ptitsyn,” sections of left- and 
right-handed helices alternate, the single helical section containing three 
monomer units in the average. As a consequence be of isotactic polystyrene 
is larger than that of atactic one.” 

Net electric charges yield an interaction energy proportional to the in- 
verse distance between the charges. That is a typical long-range force. 
However it also very strongly influences the arrangement of next neighbors 
strongly favoring configurations with largest possible next-neighbor inter- 
charge distance. These are not always the most extended chains because 
in special cases, for example, in positionally isotactic acrylic acid (ht-chain), 
the nearest charges are farther away at cis than at trans configuration so 
that b decreases with ionization and decreasing ionic strength of solution.* 
Hence in every case, the effect of free charges requires a special study of the 
basis of the actual fine structure of the chain and of the charge position. 

Dipole interaction decreases with the inverse third power of distance and 
hence chiefly influences the orientation of next neighbors. When the main 
component of dipole moment is parallel to the chain axis, electrical inter- 
action favors trans configurations with parallel dipole arrangement. As a 
consequence the square of total dipole moment of macromolecule is greater 
than the sum of squares of single moments, the value resulting in a chain 
with tetrahedral valency angle and no dipole interaction.” Typical exam- 
ples are poly-o-chlorstyrene and poly-o-bromstyrene, showing an increase 
by a factor of 1.3 and 1.1 respectively (TABLE 2). When, however, the 
main component is perpendicular to the chain axis the antiparallel czs con- 
figurations are favored, yielding a lower contribution to the mean square of 
total moment (polyvinylchloride 0.87, polyoxymethylene 0.87, polymethyl- 
methacrylate 0.73 to 0.81). The contribution, of course, is still less when 
the dipole-bearing group can at least partially rotate independently on the 


chain (polyparachlorphenylmethacrylate 0.58). Investigations on dipole 


moment in good solvents agree with those at © point, hence demonstrating 
that dipole interaction only influences the short-range structure of the chain 
and is not affected by coil dimensions." 

As already mentioned, the chief experimental method for investigation of 
short-range interaction, that is, of the effective link length, is measurement 
of gyration radius in ideal solvent. The same information may be obtained 
from small-angle X-ray scattering.”""" By plotting (rrcure 4) the rela- 
tive intensity P(d) = I(8)/I(0) multiplied by s over s = sin(d/2) a 
distinct knee is obtained at s* = 2/ka(k = 2m/X), yielding 


=i2/ ks" . (11) 
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The effective link length follows from Equation 4 and 5: 
b = bo(2a/bo — 1)" = bo(4/ks*bo — 1)™” (12) 


This method is independent of coil dimensions and therefore may be ap- 
plied also to nonideal solutions, permitting study of the changes of short- 


TABLE 2 


Ratio or OBSERVED R.M.S. DrpotE MomMENT oF MACROMOLECULE AND THAT 
PREDICTED FOR COMPLETELY UNRESTRICTED MONOMER ORIENTATION 


(wobs?/Pyo2)/? | Reference 
Poly-o-chlorstyrene —CH,—CH— 1.3 36 
Cl 
Poly-o-bromstyrene poe aie sd chan | 36 
ieee 
Polyoxyethylene —CH.—CH,—O— 0.87 37, 38 
Polyvinylchloride —CH,.—CHCl]— 0.87 39 
Polymethylmethacrylate —CH,—C(CH;)— 0.73-0.81 | 38, 40 
OOCH; 
Poly-p-bromstyrene —CH,—CH— 0.71 36, 39 
Br 
Poly-p-chlorphenylmeth- —CH.—C(CH;)— 0.59 41 
acrylate 


boo )-a 


range order in good solvents and the eventual influence of long-range inter- 
action on the structure of short-chain segments. 

Additional information may be obtained from measurement of excess 
dielectric polarization, as already mentioned above, and of optical anisot- 
ropy of macromolecules. The latter has to be derived from streaming 
birefringence. It turns out to be independent of coil dimensions, that is, 
on long-range interactions, but again larger as calculated on the assumption 
of free rotation around the chemical bond.** However, the calculated data 
of optical anisotropy are still so much affected by the uncertainty of orien- 
tation possibilities of optically anisotropic groups in the monomer that 
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they can not be used for a quantitative estimate of hindered rotation in the 
main chain although they yield a good qualitative support to results ob- 
tained by other methods. 


Long-range Interactions 


In good solvents, characterized by positive value of second virial coeffi- 
cient A, (Equation 18), the square of gyration radius increases more than 
linearly with WM. Very generally one can put 

R= ke (13) 


e being 0.25 for polymethylmethacrylate”’ in acetone at 20° C., 0.33 for cellu- 
lose nitrate” in acetone at 20° C. and 0.005, 0.07, 0.22 for polystyrene” in 


O z 4 
: Ficure 4. The product Pé over ¢ = ksa for a chain with 800 monomers and cos 
~ a=0.29. Note the knee at ¢ = 2.% 
57 
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> 
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Ficure 5. Relative increase of effective link length b/by = h/boZ'/2 = Zé/? for 
non-Gaussian chains. 
cyclohexanone at 43° and 57° C., and in toluene at 20° C., respectively. 
Such a behavior cannot be explained by hindered rotation, due to any of 
the sources mentioned above with short-range interactions, because in the 
limit of sufficiently high degree of polymerization it always leads to « = 0. 
The difference may be seen best by comparing the effective link length b 
in both cases (FiGuRES 1 and 5). Valency angle and hindered_rotation 
yield curves concave to the N axis, while Equation 13 produces an upward 


NT eu 
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curvature with the increase in b becoming more marked with increased de- 
gree of polymerization. 

In order to account for this additional coil expansion which, in the same 
polymer series, increases with molecular weight, interactions between any 
two monomers must be considered, regardless of the length of the section of 
the chain by which they are separated. In spite of the fact that (with the 
exception of electrostatic forces between free charges, which are considered 
later) all interactions between monomers are rather short-ranged; they are 
not negligible because there are chain configurations that bring any pair of 
distant monomers to direct contact. The number of possible contacts in- 
creases with the square of degree of polymerization since any monomer has 
the chance to get in contact with any other monomer. The exception to 
this is the case of closest neighbors, which are prevented from that contact 
by restrictions imposed on them by valency angle and limited rotation, that 
is, by short-range interactions. On the other side, the volume of polymer 
coil increases as R’, that is, as P*°*°”’, so that the probability of a contact 
per unit volume (contact density) increases as P*%9? As long as e is 
smaller than one third, the contact density increases with molecular weight, 
and the corresponding contribution of interaction energy does likewise. 

The energy of interaction between two monomers in a given solvent must 
express the impenetrability of monomer and eventual difference of affinity 
between the monomers and between the monomer and the solvent. Its 
precise dependence on distance r between the monomers does not matter 
very much, at least in the present stage of investigation. Therefore, one 
prefers to choose functions permitting the easiest calculation of interaction 
energy averages. Bueche” and Yamakawa and Kurata” assume (FIG- 
URE 6) 


u(r) = +0 r <.2at 
— uy exp(—3r’/2d”) r > 2a* 


a* radius of monomer, d range of interaction <A (segment length). In 
good solvents the monomer prefers to be surrounded by solvent molecules, 
the forces between any two monomers are repulsive (uw < 0). In poor 
solvents one has the opposite situation, uw > 0. At the © point the at- 
traction forces favor a just sufficiently close approach of segments to com- 
pensate for the effect of monomer impenetrability. 

The exclusion of configurations with overlapping segments yields a more 
extended coil because the excluded configurations are chiefly those leading 
to small end-to-end distance. This increase of average dimensions is even 
greater with repulsive forces, uo < 0; but it can be overcompensated by 
strong attractive forces (uw) > 0). Lattice model investigations of “ex- 
cluded volume” effect on coil dimensions were done by Wall and Mazur,” 
who by the use of an electronic computer derived the average end-to-end 


(14) 
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distance of chains and its dependence on molecular weight and average 
energy of interaction. A great advantage of this treatment is that it is not 
limited to small deviations from ideal solvent. 

A similar interaction between free segments occurs in oligomer solution 
As a consequence of interaction the osmotic pressure is not more propor- 
tional to concentration. The deviations are considered by a power expan- 
sion (Equation 18) in which the two-segment interaction is represented 
by the second virial coefficient A.. By a rather rough procedure Peterlin”™ 
was able to show proportionality of « and Az, which seems to be supported 
by experimental evidence. “However, his treatment did not yield quanti- 
tatively the proportionality factor; it is also not rigorous enough. 


oe 
Ficure 6. ae energy between two monomers according to Equation 14. 
Flory” assumed a uniform coil expansion measured by a factor | 
a= RiRe= h/he (15) 
and derived the following basic relationship 
oe — oe = 20m,(1 — 0/Pyw” (16) 
Cu = (27/2°?n') (0° /N a1) (Mo/bo)” (17) 


» = specific volume of polymer, N 4 Avogadro’s number, v; molar volume of 
solvent, M» molecular weight of monomer, ¥ an empirical constant ~Y. 
_ The same parameter enters the second virial coefficient: A» in osmotic pres- 


sure 
a/e = NakT(1/M + Ao-c + -:-) 


(18) 
(162/3°) (NaR*/M’) In [1 + (2"?/2)(a* — 19), 


caf 
T 
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leading to the power expansion 

a = 1+ (27/2)(Mo/2nbe?)*?(A2/Na)M? + -:- eC 
= 1 + -1.2992"-- ==5 


where z is defined in Equation 21. By uniform expansion the Gaussian 
distribution of intrachain distances remains unchanged. 

Nonuniform coil expansion was considered by a long series of investiga- 
tors.!" They all closely agree with the following series expansion of h 
and KR 


We = he'(1 + 42/8 +--+) = ho (1 + 1.832 + ---) 
R? = Ro'(1 + 1342/105 + ---) = Re (1 + 1.2862 + ---) (20) 
6R’/h? = 1 — 0.0472 + --- 
with 


z = 2(3Mo/2mrbe’)*?(Ao/N 4) M/F (z) 


(21) 
F(z) =1—28652+--- 


The end-to-end distance increases faster than the gyration radius. The 
expanded coil is definitely non-Gaussian. The value of a” from Equation 
19 lies exactly between those given by the first two Equations of 20 when 
only the factor 27/2 in the linear term of Equation 19 is replaced by 27/4: 
1.333, 1.299, and 1.286 are the corresponding numerical coefficients at z. 
The higher coefficients of series expansions in Equation 20 are rather un- 
certain. 

Further modifications in the treatment of excluded volume effect were 
made recently””*’ in order either to obtain higher terms in power expansions 
for R, h, and A, or to get better functional expressions valid in a greater 
range of 2, which would agree with experimental data and the lattice model 
investigations. Both yield an increase of (a° — a’)/M"” with M, indicat- 
ing a nonuniform coil expansion. A more direct proof would be highly de- 
sirable. 

Combination of X-ray and visible light scattering may provide a means 
for such a check. The latter yields the mean square gyration radius and 
the former the persistence length. At uniform expansion a and R must 
increase by the same factor when going from ideal to-good solvent, while 
nonuniform expansion requires a net increase of R and nearly constant a. 
Indeed, when the short-range interactions are unaffected by the change of 
solvent and temperature, a ought to be strictly constant. 

The scattering envelope with visible light yields an additional means of 
investigating the distribution function of intrachain distances. When the 
mean square of distance between any two chain elements obeys the same 


— 
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power rule as the end-to-end distance (Equation 13) 
(rin) = KM (22) 


with M ;, molecular weight of the chain section between the jth and the kth 
element, the ratio between gyration radius and end-to-end distance de- 
creases to™ 


6R’/h? = 1/(1 + 5e/6 + &/6) = 0/7 (23) 

25; 

20- €®0/ N/a 2 

I5- 

4 

10+ 
e | 2 
rd ; v=(ksR) 
% 6) 40 e213 OO 
# : lar d d 7 of scattering intensity P-! = cKM/R(&%) over 
3 v Bey toon Sri x ppachans Fe eae according to Equation 22. The calcu- 
4 lation is based on a Gaussian distribution of distances between any two monomers. 
4 with v = (ksR)’ and r parameter wanted in P(#), and the scattering func- 
 tion®” becomes (FIGURE 7) 
4 P(8) = (2/1 + @)ir "= e/(1 + 6),71! ae 
E- 
- Se otal = e)/(1 ae €),7]!} 
a 
4 : g . . 
4 [yz]! = / w’ exp (—w) dw incomplete factorial function 
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Its main characteristics are a persistent downward curvature of P™ plot over 
sin’(3/2) in contrast with the persistent upward curvature in the case of 
¢ = 0. The downward curvature becomes still more marked in the case 
of coils with finite number of segments. 

This check, however, is not unambiguous because polydispersity distorts 
the scattering envelope in exactly the same sense. Moreover, chains with 
¢ = 0 but with sufficiently small Z yield scattering curves that lie between 
the curve for a rigid rod and that for an infinite coil with e = 0 (FIGURE 3). 
Although the curves are definitely different from those corresponding to a 
nonvanishing ¢, the difference is not so striking that from the scattering 
envelope P ’ it could be decided with certainty whether the coil structure 
corresponds to e ~ 0, or to « = 0 and small Z, or whether the whole effect 
is due simply to polydispersity. 

A special case of long-range interaction occurs with charged polymers. 
The electrostatic repulsion potential energy between two chain elements 
having each a charge e reads 


u = e exp (—«r)/er (25) 
with 
c= 4a Dd) ce: /ekT 


Dielectric constant of solution ¢,; c; concentration; and e; charge of ions 
of the 7th kind. As a consequence of electrostatic repulsion the coil 
strongly expands with increasing ionization, that is, increasing number of 
charged chain elements and decreasing ionic strength >. c,e,”. The expan- 
sion is not uniform, the inner parts of the coil expanding more than the 
ends.” At the same time, however, the charges also influence the organi- 
zation of next neighbors; there, as already mentioned with short-range in- 
teraction, under special steric conditions, configurations leading to a con- 
traction of chain segment may be favored when this configuration increases 
the distance between neighboring charges. This possibility of simultaneous 
contraction of short segments and expansion of the whole coil requires a 
detailed study of positional and sterical particularities of chain structure in 
order to obtain the final effect of charges on coil dimensions.”*” On the 
other hand, at finite concentration of polymer solution there is always a 
marked increase of gegenion concentration inside the coil, providing a more 
efficient shielding of charges on the macromolecule.“’"* Therefore, the 
condition of constant ionic strength during dilution can be guaranteed 
only by adding low molecular ions in proper amount, for example, by dialy- 
sis with a standard solution of constant ion concentration.” Until now 
there has been neither enough experimental data at sufficiently well defined 
conditions nor a satisfying theory” available that could be checked by 
experiments, although there is a very rough qualitative agreement of rather 
rudimentary theoretical consideration with observations. 
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Solution in Laminar Flow 


Hydrodynamic interaction. Hydrodynamic interaction is a typical long- 
range interaction that decreases with inverse distance of interacting chain 
elements. Two cases are of interest: (1) the macromolecule moving with 
constant velocity through the solution at rest (pure translation: sedimenta- 
tion, electrophoresis), (2) the macromolecule in a linear flow field (rotation 
and deformation: viscosity, streaming birefringence). In both cases the 
effects are determined by the hydrodynamic resistance coefficient Ao of 
monomer and by their interaction, which in first approximation is propor- 
tional to inverse distance and also depends on the angle between velocity 
and distance of interacting units. 

Let us first consider velocities and gradients so low that the deformation 
of coil by the flow may be neglected. This is nearly always the case with 


_ translational motion and with intrinsic viscosity at zero gradient. In a 


very rough but quite satisfactory approximation, the polymer coil behaves 
like a semipermeable sphere with two characteristic parameters: diameter 
and permeability. Both may be determined from combination of trans- 


- lational resistance and intrinsic viscosity data. The sphere diameter yields 


_ determinations by other methods. 


~ coil in ideal solutions. 


the gyration radius of coiled macromolecule in fairly good agreement with 

However, one may also start with average coil configuration in studying 
the effect of hydrodynamic interaction. In order to do that, the averaged 
values of all intrachain distances must be known. Most calculations are 
based on Equation 3, and therefore their results are applicable only to the 
*8.79-8 By taking into account the unavoidable poly- 
dispersity, it was possible to use the results for determination of hydrody- 


a . . . . . . . 83 
namic resistance Ao in satisfactory agreement with geometric expectations. 


Treating the case with nonvanishing excluded volume effect™’” has also 


proved successful, so that the hydrodynamic interaction at translation and 


rotation of undeformed coil is fairly well understood. The results were 


62,78 


extensively reviewed quite recently so that it seems justified here to 


omit a treatment of this case. Let us consider in more detail, however, 


the effects when chain deformation and orientation are of primary concern, 


' that is, the gradient dependence of intrinsic viscosity and streaming bire- 


fringence. ee 
With increasing gradient G, the observed intrinsic viscosity decreases in 


close agreement with the series expansion™ 
[n] = [nlo(1 — B’G? + ---) = [nlo(1 — BBo’ + ---) 
Bo = M[nlonG/N skT 


(26) 


4 and it is very likely that it approaches a new limit [ny]. < Ino , although no 
- experiment has yet demonstrated it. Due to the rather difficult extrapola- 
‘tion to zero concentration, there is also some uncertainty as to whether the 
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initial drop of intrinsic viscosity is indeed proportional to the square of the 
gradient, although no clear-cut case demonstrates a different type of de- 
pendence. 

Theory always yields a quadratic term in the series expansion Equation 
26, but there is no final agreement about the value of B, and still less about 
higher coefficients. Moreover, all calculations are based on Gaussian coil 
and therefore can be applied only to ideal solutions. 

Without any hydrodynamic interaction, intrinsic viscosity” 


[n] = (NaR*/6Mo)A (27) 
turns out to be independent on gradient (FIGURE 8). Isotropic hydrody- 
namic interaction changes Ao into Ao*, which, in first approximation, may 
be written 

Ao* 


Ao/(1 + e) 
(2/27 n°)? AgP!?/b 


More detailed treatment of hydrodynamic interaction” using Rouse’s 
model”’ slightly modifies Equation 28 without changing the functional con- 
nection between [n] and A as the average value of resistance coefficient of 
monomer in laminar flow.”'*’ This reduction is the consequence of two 
opposing effects. In laminar flow the ends of coiled macromolecule have 
average velocities that are equal in amount but opposite in direction. 
Monomers of one half of the chain therefore move in the same direction 
and by hydrodynamic interaction reduce Ao of any monomer of this same 
half of the chain. Since monomers of the other half are moving in opposite 
direction, they increase Ao, their effect being smaller because they are 
farther away from the monomer under consideration. 

With increasing 6 the coil gets more and more extended 


h® = ho (1 + 26°/3) (29) 


However, the extension is not uniform,” since the distances between re- 
mote elements increase more rapidly than those between closer neighbors. | 
As a first consequence, the influence of the more distant parts of the coil 
rapidly vanishes so that the reduction of Ag to Ao* becomes still larger. 
Therefore, intrinsic viscosity decreases as” 


[n]* = [nlo*[1 — 0.152¢e6°/(1 + e) + +5] (30) 


At still higher coil extension, the contributions of the close-by neighbors 
also markedly decrease and yield a rather slow increase of [n]* (FIGURE 8). 
The effect is limited by the model to such 8 values that h < Pb). When h 
reaches the extended length Pb no further expansion can occur so that [n]* 
has to decrease and to approach a limit [n]..*, roughly corresponding to that 


of rodlike particles of similar length-to-thickness ratio at maximum orienta- 
tion.® 


(28) 


é 
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Hydrodynamic interaction is not isotropic, and its average strongly de- 
pends on space distribution of all chain segments. As a consequence Aj, 
for longitudinal (radial) movement of chain segments is smaller than AZ, 
for transversal (rotational) motion. Intrinsic viscosity at zero gradient 
remains nearly unchanged, the coefficient Ao* merely being replaced by 
Aj,. But the increase in rotational resistance results in better orientation, 

_ yielding smaller extinction angle x of streaming birefringence. The de- 
"crease in longitudinal resistance, however, results in a smaller extension of 
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< Figure 8. Relative intrinsic viscosity [n]re1 = [n]/[n]o for free draining coil, coil 
~ with hydrodynamic interaction according to Equation 80 with e/(e + 1) = 1({m]re1), 
and for coil with finite inner viscosity ([n]**) according to Cerf.** The broken line 


represents Copi¢’s measurement on polystyrene.*° 


* the coil in flow. ‘Therefore, even in the case that the direct mfluence of 
coil expansion on hydrodynamic interaction is neglected, [n] decreases 


[nlo* = [nlo(1 — BoB’ + ---) (31) 


The coefficient B, is 0.148 according to Copié and Peterlin*””""* and 0.0196 
according to Ikeda.®” The difference is due to a different way of averaging. 
In order to include the effect of coil expansion on hydrodynamic interac- 
tion, it is necessary to multiply Equations 30 and 31, that is, to add the 
respective coefficients B of the linear term. The resulting initial drop in 
[n] is greater than according to either Equation 30 or 31. 
In order to explain the steady decrease of [7] with G, it is necessary to 
- make the additional assumption that in laminar flow the coil does not ex- 
_ pand as much as given by Equation 29. Such a behavior may be easily 
understood. During one full rotation with the period T = Ar /G, the coil 
~ must change its shape four times from maximum expansion to maximum 
- compression or vice versa, and in so doing, a great many changes in angular 
orientation of adjacent monomers must take place in a rather short time. 
Any restriction in free rotation around the valency bond opposes such a 
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sudden change. This short-range interaction effect can De described as a 
force proportional to the rate of shape change of the coil ~ 


—F dh/dt (32) 


The coefficient F of “inner viscosity’’ has the dimension of viscosity times 
length. In order to account for the change of extinction angle of streaming 
birefringence with the gradient” the factor F must be independent of degree 
of polymerization. With such an assumption alone and neglecting the 
change of hydrodynamic interaction with the still remaining coil expansion, 
one obtains the curve [y]** (FIGURE 8). For a complete fit or comparison 
with experiment one will have to include the effect of anisotropy of hydro- 
dynamic interaction and the influence of actual coil extension which, al- 
though smaller than according to Equation 29 is not negligibly small. 

Some experiments with extremely viscous solutions, lubricating oils at 
low temperature, 7 ~ 10° P, containing small amount of polymer added as 
improver,”’ show a gradient dependence that is strikingly similar to that of 
[n|*. The experiments, of course, were performed at too high concentra- 
tions, so that they do not yet represent the properties of isolated macro- 
molecule, and it is always hard to estimate how much of the observed effect 
is due to interaction effects. Nevertheless, these experiments are an indi- 
dication that at special conditions, that is, extremely viscous solvent and 
low gradient, one may find a situation where the macromolecular coil mark- 
edly expands under the influence of laminar flow because the high solvent 
viscosity seems to overshadow completely the influence of “inner viscos- 
ity.” At 6 values, high enough for marked gradient dependence of vis- 
cosity, the gradient @ is still so low that the macromolecular coil has no 
special difficulty in changing shape, and the forces exerted by the solvent 
so high that the ‘inner viscosity” resistance may be completely neglected. 
As a consequence, the intrinsic viscosity first decreases and then increases 
in qualitative agreement with the curve [y]*. Further experiments at lower 
polymer concentration will be necessary before the soundness of this pre- 
liminary extrapolation can be evaluated. Similar experimental evidence » 
was presented by Merrill’’ on polyisobutylene solutions in benzene, a poor | 
solvent just a few degrees above the 9-point. Because in this case the vis- 
cosity of the solvent is not high, it must be concluded that PIB molecules in 
this solvent must be extremely soft, so that they easily expand. Again, 
the solutions were too concentrated and it could not be ascertained how 
much of the observed effect is due to the contribution of a single macro- 
molecule and how much to intermolecular interaction. 

The “‘inner viscosity”? concept satisfactorily explains the observed initial 
change of extinction angle x with gradient (FIGURE 9). With sufficiently 
viscous solvents, Equation 33 is obtained?” 


— (dx/d@) g-0 = B+ Cno (33) 
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The intercept B is proportional” to molecular weight of solute (FIGURE 10). 
When, however, the solvent viscosity is too low or the macromolecule too 
rigid the gradient dependence of extinction angle becomes the same as in 


_ the case of rigid particle: the coefficient B vanishes, and the curve 


(dx/dG) a—o 
over mo goes through the origin, the slope C being much higher than in 


' Equation 33 with B ¥ 0. Such a case occurs with deoxyribonucleic 


99 ,100 


acid (FIGURE 9), which is nearly rigid in water and relatively flexible 


in water-glycerin mixtures (7 > 2cP). Theory” yields for not too low sol- 
vent viscosity 


~ B=0.0062FMb'/M kT, C =0.2M[nho/NakT small [eae 2 34) 
0.0045 0.1 strong |interaction 
and for nearly rigid particles (mo > 0) 
B =0,C = 0.9M[nh/N skT small a on 
0.705 strong} interaction 


Optical interaction. In an external electric field, the induced dipoles in 


- the macromolecule interact and change the field value inside the coil. The 
interaction is proportional to the third power of inverse distance between 
- the interacting chain dipoles. The effects are so small that in most cases 


ee 
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they can be neglected, for example, in dielectric polarization and in light 


scattering when depolarization of scattered light is not in question. The 
effects, however, do matter in birefringence measurements where they in- 


~ fluence the relatively small difference of main refractive indices in the ori- 
ented and deformed macromolecule. 


The average over-all macromolecular configurations leading to the same 


end-to-end vector 7:2 yields a nearly ellipsoidal monomer distribution in- 
~ stead of spherical one that results merely at vanishing end-to-end distance. 
As a consequence, the macromolecule exhibits a shape birefringence, 
- which according to Copi 


101,102 


6”! reads 


r { 9M \[m(n? — w) 0 
(m — ¥2)5 = Osh/ho = he! (7) == (35) 


-n refractive index, p density of macromolecule. This contribution has to 


be added to intrinsic birefringence” 


3 2 
(ya — y2)s =Oer"/ho’ =F (on — aa) (36) 


3 a , 2 polarizability of chain segment in chain direction and perpendicular 


to it, in order to obtain the total optical anisotropy of macromolecular coil 
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with the end-to-end distance r. Intrinsic anisotropy may be either positive 
or negative depending on the optical properties of monomer unit. Form 
anisotropy, however, is always positive. 

The resulting birefringence in laminar flow turns out to be 


An = (4tN4/3rmM) (mo + 2)6[0.86 0, + (1 + 6”)"6,) (87) 


_ The contributions of both kinds of optical anisotropy vary with the gradi- 
- ent or 6 ina different way. The corresponding effects are most conspicuous 
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fs Ficure 12. Streaming birefringence An/cG of polyacrylic acid solutions!"® at 
 @ = 19600s— and ionization degree 0.5 as function of added Nat: (a) polymer con- 
centration 0.92 X 107%; (b) 2.36 x 10° gm./em.’ Reproduced by permission of 
' Academic Press, Inc.1 


4 with systems exhibiting negative intrinsic birefringence, for example, poly- 
E. styrene and polyacrylic acid. In dioxane with considerable difference be- 
- tween nm and n a high molecular PS™ (M = 8.10°, [m7] = 930 cm.*/gm.) 
4 yields a positive streaming birefringence at small gradient, indicating a 
greater absolute value of 6; than of 4 (riguRE 11). With increasing gradi- 
4 ent the factor (1 + 6”)”” lets the influence of intrinsic anisotropy prevail 
- over that of form anisotropy. The birefringence, after reaching a maxi- 
- mum, drops, and becomes negative (FIGURE 11). From experimental data 
Shem alues 6, =.312.10." em.*, 0, = —78.10,” em.',.0, > a2 =_=13010 

em.’ may be calculated.” 
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As a consequence of contribution of form anisotropy the initial specific 
streaming birefringence 


(gaa) - se ( 3N0 kT ( i ay a 


changes with shape of macromolecular coil or/and with the solvent. The 
first case is easily observable with polyelectrolytes in which added salt 
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Figure 13. The birefringence An/cg at G = 1000s~! for polymethylmethacrylate 
in several solvents of different index of refraction.!8 Reproduced by permission of 
Academic Press, Inc.!"! 


increases the ionic strength, by electrostatic shielding reduces the repulsive 
forces in the coil, and hence reduces its size. Polyacrylic acid in water solu- 
tion’®* changes the sign of birefringence from negative to positive with in- 
creased Na* concentration (FIGURE 12). At low salt concentration, the 
coil is rather extended so that the negative intrinsic birefringence prevails. 
With more added salt the coil contracts, 6; increases rapidly due to the 
factor ho ° and résults in positive streaming birefringence. 

The theory also satisfactorily explains the solvent effect, °""™ yield- 
ing a nearly parabolic curve for (An/enonoG) co plotted over n¢ (FIGURE 
13). In the case of negative intrinsic anisotropy, the birefringence van- 
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ishes for two distinct values of solvent refractive index. The agreement in 
most cases is fairly good. However, some effects of molecular weight on 
gradient dependence of streaming birefringence’’—the sign change in 
the system polystyrene-dioxane occurs at higher @ with lower molecular 
weight—cannot be explained by Equation 37, which does not consider the 
partial coil rigidity and the limitations of coil extensibility by the require- 
ment that the r.m.s. end-to-end distance must be always smaller than the 
extended length of the macromolecule. Exactly the same situation occurs 
with the gradient dependence of intrinsic viscosity °""™ 
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STATISTICAL THERMODYNAMICS OF COILING-TYPE 
POLY MERS* 


Frederick T. Wall and Jacob Mazurt 
Noyes Chemical Laboratory, University of Illinois, Urbana, Il. 


Introduction 


During recent years a considerable amount of work has been carried out 
on the thermodynamics of coiling-type polymer molecules."” From a sta- 
tistical point of view, one of the important considerations deals with the 
over-all dimensions of such molecules, taking into account the fact that any 
given molecule can assume numerous configurations. An appreciable 
amount of statistical work, including Monte Carlo-type calculations using 
high speed computers, has been directed toward establishing the mean 
square end-to-end separations of flexible polymer chains. This Monte 
Carlo-type work has been carried out by generating nonintersecting ran- 
dom walks in various lattice systems.“ 

The model assumed for the nonintersecting random walks, in which 
double occupancy of any site is forbidden, implies an infinite potential en- 
ergy when two elements of the chain are superimposed. It was also tacitly 
assumed, in the work carried out heretofore, that the potential energy was 
constant (for example, zero) for all configurations except those involving 
double occupancy. Thus two values for the intrachain interaction energies 
were implied, infinity and zero. 

In this paper we shall report the results of carrying out somewhat more 
elaborate calculations, assuming a potential energy function different from 
that implicit in the earlier Monte Carlo-type calculations. As before, we 
shall stipulate that double occupancy, or approaches within a certain speci- 
fied short distance, will correspond to infinite energy. We shall further 
assume that there will be a short range of intrachain element separations, 
just outside the above-mentioned excluded range, for which there will be a 
finite energy of interaction different from zero. Finally, for intrachain ele- 
ment separations beyond the range of nonzero finite interaction energies, 
we shall assume that the energies are zero. In this way our model becomes 
slightly more sophisticated than before, because we now assume three (in- 
stead of two) possible interaction energies, namely, infinity, a constant 
different from zero, and zero. This does not mean, of course, that the en- 
ergy of a given molecule cannot assume more than three values. Actually 
the energy of a molecule will range from zero (for an extended configura- 
tion) to an upper or lower bound (depending upon the sign of the finite 
interaction term) for molecules that are tightly coiled. To establish the 
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potential energy of a molecule, it is necessary only to count the number of 
intrachain interactions with energies different from zero that occur ina 
given configuration. The use of this model should introduce some interest- 
ing new features into the problem of coiling-type polymer dimensions. In 
particular, it should provide a means of learning more about the effect of 
temperature on the size and shape of molecules. 

At this juncture, we should point out that if the intrachain interactions 
were restricted to near-contour neighbors only, an analytic solution of the 
problem might be obtained through consideration of Markoff chains. A 
solution of this problem was attempted by Montroll,’ who treated random 
walks in a two-dimensional 90° lattice with exclusions limited to loops of 
four steps. Various other mathematical procedures have been adopted by 
a number of authors to deal with the problem of configurations of real poly- 
mers subject to interaction energies. However, these mathematical treat- 
ments are not only complex, but they involve serious approximations, the 
validities of which cannot be readily checked. The simplified problem 
of configurations of coiling chains subject to the excluded volume effect 


_ belongs to this category of problems for, as indicated above, the excluded 


volume concept can always be regarded as corresponding to hard sphere 
interaction potentials between the elements of a polymer chain. 

The computations were carried out in the following way. As with the 
earlier Monte Carlo work, we generated random walks in various lattice 
systems, rejecting those walks involving either double occupancy of lattice 
sites or, in some instances, approaches within a finite distance such as those 
corresponding to nearest-neighbor sites. In addition, we counted the num- 
ber of intrachain approaches lying in a specified middle-distance range. 
Having counted the number of such interactions, each of which would con- 
tribute to the molecule a finite amount of energy different from zero, it 
was possible to classify the configuration by its energy value. All this work 
was carried out using the Illiac, a high-speed digital computer, at the Uni- 
versity of Illinois. 


Theory 


In the following discussion, we shall assume that the polymer chain con- 
sists of n links or n + 1 chain elements. A given polymer configuration 
will be characterized by its end-to-end separation, 7, , and by a number, z, 


- which equals the number of intrachain interactions corresponding to the 


afore-mentioned intermediate distance range with potential energy of in- 
teraction, «. The total potential energy of such a molecule will therefore 
be ze, and the partition function for the system will be given by 


Q = DUN, exp(—2e/kT) } (1) 


where N, is the number of ways a molecule can exhibit configurations with 
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z interactions of the type specified (for our purpose, we shall let NV. , which 
also depends on n, equal the number of samples found by the Monte Carlo 
technique). ; ; 

Expressions may now be written for certain average physical properties. 
For example, the mean square end-to-end separation will be 


>> (rn’).Nz exp (—2e/kT) 
(ee (2) 
: Q : 
where (7,,). is the mean square end-to-end separation of molecules with z 
interactions. 
The average energy, (EZ), will equal (z)e and will be given by 


> zeN, exp (—2e/kT') 
(E) = @)e = =—_, = 
20 In Q 

oT 


= kT 


The free energy will, of course, be given by 
F = —kTmQ (4) 


and the thermodynamic potential associated with a chain element (or 
“partial-link” free energy) will be 
OF dlInQ 
= — = —kT 
ean on 


ee eee 


Q z On 


(5) 


exp (—2e/kT') 


Computations 


Our Monte Carlo calculations, performed on a high-speed computer, all 
involved random walks that were carried out, subject to appropriate rules, 
in several different lattice systems, both in two and in three dimensions. 
Six different sets of computations were carried out, four of them in two- 
dimensional lattices and two in three-dimensional lattices. The essential 
characteristics of these lattices and the distance ranges to be associated with 
intramolecular interaction energies are summarized in TABLE 1. The first 
two entries (a, b) deal with walks in the two-choice square lattice, the sec- 
ond involving a wider range of distances associated with the finite energy 
of interaction. The third and fourth entries (¢ and d ) both deal with walks 
in the two-choice hexagonal plane lattice, once again with interaction en- 
ergies identified with different distance ranges. The fifth entry (e) deals 
with walks in the four-choice (90° bond angle) three-dimensional simple 
cubic lattice. The last (f) deals with the three-choice tetrahedral lattice. 
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The columns in TABLE | are, for the most part, self-explanatory. The next 
to the last column, ‘Possible number of intermediate neighbors,” is essen- 
tially a type of coordination number applicable to interior points of a chain. 
Because of the constraints imposed, an interior point in a walk of type c, 
for example, can have only one intermediate neighbor, not counting those 
that must lie within the specified range because of the continuous contour 
of the chain (the end point of a c-type chain can have two intermediate- 
range neighbors). For each entry, the last column denotes the initial 
number of samples, subject to attrition, employed in the calculations. 

For all the computations we assumed the step lengths or distances be- 
tween adjacent lattice sites to equal unity. Some of the random walks 
(those designated c, e, and f in TABLE 1) were generated by using the simple 
techniques employed in the earliest studies carried out in this laboratory.’ 
This procedure consisted essentially of generating a string of unit vectors, 


TABLE 1 
Z See ees 
. No. of No. of Excluded ntermediate o. of in- 
Desig- | dimen- pene lattice range range Senienge termedi- N 
nation | ‘sions | 228'© | choices (ejj = @&) (eij = € ~ 0) ee ate neigh- 
bors 
a 2 90° 2 0< rig <1) 1< ji < V2 W?2 < rj < @ 2 200,000 
b 2 90° 2 OK te < tl Nee | 2 ey << e 4 200,000 
c 2 120° 2 0< rip <1) 1 < ij < V3 V3 < ij < © 1 50,000 
d 2 120° 2 OC <1) 1<t <2 | 2X47 < @ 5 60, 000 
e 3 90° 4 0< rey <1) 1S ey < V2 WV < Ti < @ 4 100,000 
a 3 Tetra- 3 OS rp Let <r 2 2<nj © 16 220,000 
hedral 


which in the aggregate would represent a molecular chain. The individual 
vectors were chosen by a pseudorandom number technique and, after each 
selection of an additional vector, the chain was tested to make sure that the 
last-added vector did not terminate on a site already occupied. When such 
exclusion was encountered, the sample was rejected and a new chain was 
begun. Other walks (designated a, b, and d in TABLE 1) were generated 
by the “stride” technique descriked by Wall et al.° 

In carrying out the computations it was necessary to count the number 
of intrachain interactions falling within the prescribed intermediate range. 
Once the value of the interaction energy has been fixed, this count deter- 
mines the energy of the chain. From this it is easy to establish the parti- 


‘tion function and the other derived quantities by making use of Equations 


1 through 5. oe 
The results of the calculations here reported are limited to values of 


n < 60. Collection of statistical data for polymers with a greater number 
of steps was seriously limited by time considerations for the following rea- 
sons: (1) the number of tests that the computer must perform-to check all 
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possible intralattice distances corresponding to nonzero interaction energies 
is of the order of n’; (2) the rejection of random walks with double occupan- 
cies seriously diminishes the computer yield of chains with a large number 
of steps (the procedure of Wall and Erpenbeck’ for chain enrichment had 
not been fully developed when the present work was done); and (3) the 
spread of possible end-to-end distances increases rapidly with increasing n, 
and more data are required for large n to eliminate excessive scattering of 
the data. na 

The data have all been analyzed through further computations involving 
the use of Equations 2and 4. The reliability of the data diminishes as we 
depart from values of « = 0, because the relative weights of the samples are 
changed. For positive values of e, smooth graphs can be obtained with e¢ 
less than 0.75 kT for smalln. However, for e > 0.5 kT there is serious scat- 
tering of data when n is greater than 40. 

For model polymer chains with negative values of e, we could not extend 
our calculations much below « = —0.5 kT’ because of the scattering of the 
results. The lowest value of e that still yields reliable results depends very 
strongly on the maximal number of intermediate-range neighbors and on 
the maximal number of choices in the propagation of the random walk. 
The reason for this lies in the fact that coiled molecules will be identified 
with high Boltzmann weight factors even when the coiled molecules corre- 
spond to configurations whose yield from the computer is small. 


Discussion 


The results of our computations are presented in graphical form. In 
FIGURES 1 and 2 are plotted the quotient of the mean square end-to-end 
separation, (r,’), by the number of steps, n, as a function of n. The curves 
in FIGURE 1 refer to two-dimensional walks carried out in the two-choice 
(90° bond angle) square lattice, designated a@ in TABLE 1. The several 
curves in FIGURE | correspond to different values of the parameter ¢«/kT. 
It will be observed that as e/kT increases the afore-mentioned quotient in- 
creases for a given value of n. This, of course, is to be expected since a 
repulsive potential will tend to elongate the molecules. For negative values 
of e/kT, the chains become more compact relative to the simple models that 
exclude double occupancy but for which e = 0. It will be further observed 
that none of the curves in FIGURE 1 exhibits values of the ordinate less than 
unity, for reasons discussed below. 

FIGURE 2 is a similar graph for four-choice (90° bond angle) restricted 
walks carried out in a simple cubic three-dimensional lattice. Qualita- 
tively, the considerations cited with respect to the two-dimensional system 
described above are once again valid, except that when ¢/kT is sufficiently 
negative the values of (r,”)/n can be less than unity. Asa matter of fact, 
when ¢e/kT’ = —0.50 the function (7,”)/n appears to be constant and equal 
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to unity for all values of n. This is an interesting phenomenon, correspond- 
ing to what is observed when the pressure-volume product of a nonideal gas 
is plotted against pressure at various temperatures. At low-enough tem- 
peratures, the slopes of such curves are negative, whereas at high tempera- 
tures they are positive. At a particular temperature, known as the Boyle 
temperature, pV is constant. In precisely the same way, (r,")/n is constant 
at what is substantially the Flory © temperature’ which, for our simple 


0.75 


O25 


-0.625 


20 40 60 
- N 


- Fieure 1. Graphs of (r,2)/n versus n for chains with 90° bond angles in the two- 
dimensional square lattice (designation a in TABLE 1) for various values of «/kT in- 


' dicated next to curves. 


three-dimensional model, would equal —«/0.50k. At such a temperature, 
the attractive forces just compensate the hard-sphere repulsions that are 
implied by exclusion of double occupancy. 
At this point, it may properly be asked why no theta temperature is 
observed for the two-dimensional model while it shows up in the three-di- 
mensional model. The reason is clear from simple geometrical considera- 
tions. In two dimensions, the most tightly packed molecule would corre- 
spond roughly to a spiral with an area proportional to (rn )s This same area 
would also be substantially equal to n, the number of points in the chain, 
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assuming unit distances between points in a square lattice. Hence (r,°)/n 
would be constant in two dimensions for the most tightly packed form of 
molecule. Such a molecular state could be realized at finite temperatures 
only if the attractive potential were infinite. Alternatively, for a finite 


ay) 


ie) 20 N 40 60 


Figure 2. Graphs of (r,2)/n versus n for chains with 90° bond angles in the three- 


dimensional simple cubic lattice (designation e in TABLE 1) for various values of 
e/kT indicated next to curves. 


attractive potential, the two-dimensional, closely packed spiral could exist 
only at absolute zero (one could therefore say that absolute zero is the 
theta temperature for a two-dimensional polymer). In three dimensions 
the situation is somewhat different. The most closely packed random coil 
would be like a tight ball of yarn with volume proportional to (r,)*, ° The 
volume of this sphere would likewise be substantially equal to n, the num- 
ber of lattice sites occupied in a simple cubic lattice. This would imply 
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that (r,”) is proportional to n”* and hence that (r,)/n is proportional to 
ag indicating that with an infinite attractive potential the quotient that 


has been plotted would diminish with increasing n. Hence there should be 


VAT 


Fieurn 4. Graphs corresponding to those of Figure 2 drawn on a log-log scale. 


a finite attractive potential at a temperature other than absolute zero, for 
_ which the slope of the (r,’)/n plot would be zero. This is, in fact, found 
to be the case in the three-dimensional lattice. 
3 When the logarithm of (r,)/n is plotted against the logarithm of n, rea- 
 sonably straight lines are obtained, as indicated in rigurEs 3 and 4. These 
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curves reveal that, except for low values of n, all such plots are linear. This 
means that a relationship of the form 


(rn?) = an” (6) 


provides an adequate description of random-walk polymers even with in- 
teraction energies. 

Graphs for the other kinds of walks described in TABLE 1 have also been 
prepared; they are not reproduced here since no significantly new features 
arise that would alter the general discussion or interpretation. It should 
be mentioned, however, that the values of the exponent b depend mainly 
on the number of dimensions and are little dependent on the lattice details. 
For category c, the molecular properties do not exhibit a strong dependence 
on e, except for extreme values, since the possible number of intermediate 
neighbors is only equal to 1. For the tetrahedral lattice, denoted by f, 
some minor differences in behavior were encountered, because the excluded 
range covered not only occupied lattice sites but also nearest-neighbor sites 
adjacent to the chain. This subject is being investigated further. 


Thermodynamic Properties 


By using Equation 4 one can readily calculate the free energies of coiling- 
type molecules for different interaction energies. What we actually calcu- 
lated is a difference in free energy, F — F°, choosing as our standard state 
that corresponding to « = 0. For the two lattices already discussed, 
namely, the two-choice square lattice in two dimensions and the four-choice 
simple cubic lattice in three dimensions, this difference in free energy is 
found to be a linear function of the number of bonds. This means, in ac- 
cordance with Equation 5, that the partial-link free energy, or the chemical 
potential per link, is independent of n. As e/kT' increases, the free energy 
likewise increases; for negative values of «, F — F° becomes negative. This 
behavior results principally from the contribution of the energy of the 
system, which more than compensates for any entropy effect. Regardless 
of the sign of ¢, the entropy difference, S — S°, is never positive, because 
the new distribution resulting from the introduction of an interaction po- 
tential will necessarily be less random than the old. It will be observed in 
FIGURES 5 and 6 that the lines representing F — F° as functions of n fan 
out in both directions from the line corresponding to e = 0. The behavior 
may be empirically represented by Equation 7 


F — F° = (n — 2)AkT In (1 + Be/kT) (7) 
where A and B are constants. 


From Equation 7 it is easy to establish equations for S — S° and for 
H — H° (which will be the same as E — E’), as follows: 


PM Ae 
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S— 8 = —(n — 2)Akin (1 + Be/er) + % = PAB/T (g) 


1 + Be/kT 
i eS (n — 2)ABe 
1 + Be/kT ° (9) 
Ke) 
0.875 
8 — 
0.75 
0.625 
0.5 
0.375 
Ss 0.25 
0.125 
bK 
ink 
tis Z 
©) (@) 
oF 
AS 
-0.125 
a -0.25 
-0.375 
-05 
-0.750 ~0625 
hy 20 40 60 


Figure 5. Graphs showing how the free energy of chains in two dimensions 
(designation a) depends on n for different values of ¢/kT. 


It is evident that both S and H are linear functions of n. H — H° was 
also calculated directly on the computer by using Equation 3, and the linear 
dependence of H on n was confirmed. The aforementioned linear depend- 
ence of the several important thermodynamic functions of n was established 
for all six designations, a through f, including the tetrahedral system. 
These thermodynamic properties are established with considerable ac- 
curacy, even up to n = 60, since no excessive scattering of the data occur, 
as was the case with computations of the mean square end-to-end separa- 


tions. 
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Ficure 6. Graphs showing how the free energy a chains in three dimensions 
(designation e) depends on n for different values of e/kT 


TABLE 2 
Designation (s)o B 
n—2 . 
a 0.285 0.288 
b 0.382 1.131 
d 0.463 0.979 
é 0.274 1.428 
f 0.865 1.500 


Another interesting result may be deduced from the thermodynamic 
computations. Since we are dealing with a single-step type of potential 
energy function, Equation 9 may be rewritten in terms of the average num- 
ber of Sriteadtain interactions corresponding to the given value of e: | 
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_ AB(n — 2) 
OSes ay) 
If « = 0, then 
(z)o = AB(n — 2) (11) 


where (z)o is the average number of intrachain distances lying within the 
specified intermediate separation range for an excluded volume polymer 
with « = 0. 

For small values of Be/k7T, the logarithmic term in 7 can be linearized 
and, with the help of 11, we obtain 


F — F° = (z)oe (12) 


Therefore the free-energy difference changes linearly with « when the po- 
tential energy of intrachain interaction is sufficiently small. Moreover 
when e¢ is small, it is also clear that 


H — H° = (zoe (13) 
and 
S—S=0 (14) 


Actually, S — S° must be negative, but this shows up only with higher- 
order terms in e. TABLE 2 summarizes the pertinent thermodynamic pa- 
rameters. 
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THERMODYNAMICAL EXPLANATION OF LARGE PERIODS IN 
HIGH POLYMER CRYSTALS AND DRAWN FIBERS 


E. W. Fischer 


University of Mainz, Mainz, Germany 


Introduction 


Nearly all crystallizable high polymers show one ormore structural 
periodicities between 50 and 500 A, revealing themselves in small-angle 
X-ray scattering. The occurrence of long-period diffraction is always asso- 
ciated with the existence of crystallites. Therefore it is very reasonable to 
assume that the origin of the long periods is the regular placement of the 
crystallites. I shall call the process leading to this regular arrangement 
“periodic crystallization.” This paper deals with the reasons for the spe- 
cial process of crystal growth that occurs in polymer systems. For this 
purpose the experiments that led to the introduction of the concept of 
periodic crystallization will be considered. Then some experimental fea- 
tures will be discussed in order to decide whether the cause of periodic 
crystallization is kinetical or thermodynamical. Finally the fundamentals 
of a theoretical treatment developed by A. Peterlin and myself will be 
considered. This theory maintains that the phenomenon is caused by the 
anisotropy of lattice forces in chain crystals. 


Periodic Crystallization 


Often the structure of the high-polymer system causing small-angle scat- 
tering is not known exactly. In three cases, however, the associated struc- 
ture has been revealed by means of electron microscopy. These scattering 
effects are: (1) the numerous and sharp small-angle interference patterns 
of certain natural fibers such as collagen; (2) the meridional reflections of 
synthetic drawn fibers; and (3) the small-angle interferences that originate 
from polymer single crystals grown from solution. 

In the first case it was shown’” that the periods measured by X-ray in- 
vestigations correspond to a system of regular transverse stripes within the 
collagen fibrils. Probably the period is due to a certain order of the amino- 
acid residues within the protein chains.* This type of periodic structure 
in which large periods are caused in one way or another by the structure 
of the molecules will not be discussed in this paper. 

However, the large periods obtained from synthetic drawn fibers cannot 
be explained by a periodicity within the molecule itself. According to 
Hess and Kiessig* the small-angle maxima are produced in this case by a 
periodic change of crystalline and amorphous regions.* This idea is gen- 


*It must be remarked that some authors®:® believe that the large period in pro- 
tein is also caused by alternation of ordered and disordered regions. 
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erally accepted today although there are differences of opinion about the 
inner structure of the amorphous areas.’ Fiaure 1 shows the model of 
fiber polymers proposed by Hess and Kiessig* and constructed by Statton.” 
All the details of this model need not be accepted; for further consideration 
it is only essential that the fiber polymers consist of elementary fibrils 
having nearly periodic alternation of poorly and well-ordered regions. 


Figure 1. Model of fiber structure with approximate periodic alternation of 
erystalline and amorphous regions (W. O. Statton).8 Reproduced by permission of 
Journal of Polymer Science, vol. 41, copyright 1959 by Interscience Publishers, 
Inc., New York, N. Y. 


The third case mentioned above deals with a cake of polyethylene single 
crystals. Here the reason for the small-angle period is obvious. During 
the crystallization of polyethylene from dilute solution, single crystals 
grow into the shape of a lamella in which the chains are perpendicular to 
the flat surfaces.” The average length of the chains is about 50 times 
greater than the thickness of the single crystals: therefore Keller! con- 
eluded that the molecules must be folded at the outer planes of the lamella. 
Such lamellae when arranged in layers form a cake, and Keller and O’Con- 
nor showed that this cake gives a small-angle period up to the fourth 
order, corresponding exactly to the thickness of the platelike crystals. 

In the alternation of crystalline and amorphous regions in drawn fibers 
as well as in the folding of the chains that occurs during the crystallization 
from dilute solution, the structure of the crystal is subjected to periodic 
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lattice distortions in chain direction. Furthermore investigations have 
shown that in polycrystalline high polymers crystallized from melt, from 
concentrated solution, or from the glassy state the crystallites have a small 
and uniform length in the chain direction. Electron-microscopy studies 
of the surface of bulk-crystallized polymers also often show a lamellar 
structure.”"* For special reasons, which cannot be considered in detail 
here, I am rather certain that the whole polymer solid consists of crystals 
similar to those seen on the surface. For example, the replicas of fracture 


Figure 2. Fracture surface of high-density polyethylene, crystallized at 125° C. 


surfaces from polyethylene show clearly the individual lamallae (riguURE 
2).* This statement does not necessarily mean that in crystallization in 
bulk the chains are folded at the boundary surfaces of the crystallites. 
However, crystallites are always found to have the same thickness and to 
be arranged in a particular order with one another; the folding of the mole- 
cules is the simplest assumption to explain this morphological structure. 
Independently of this assumption I shall name the process that leads, at 
least in small regions, to a periodic arrangement of crystallites with nearly 
the same size ‘periodic crystallization.” I believe that this special man- 
ner of crystallization is a general feature of all crystalline polymers. 


*See also the fracture experiments of Bunn et al.*! with polytetrafluorethylene | 
and of Geil* with polyoxymethylene. 
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Kinetic or Thermodynamic Causes 


The cause of periodic crystallization is not yet known. There are two 
principal possibilities to explain the phenomenon. It may be that the 
periodic structure is kinetically caused. On the other hand is it conceiv- 
able that crystals with such a small thickness in chain direction can be in 
a state of lowest free energy for the polymer material. 

The first idea was extended in the theories proposed by Lauritzen and 
Hoffmann” and independently by Price.’® In the case of crystal growth 
from dilute solutions these authors assume that in the outer planes of the 
nucleus the molecules are already folded and therefore further growth of 
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Figure 3. Dependence of the long period of drawn 6-nylon on the temperature 
(K. Kobayashi, private communication). 


the crystal is kinetically retarded in chain direction. I shall not discuss 
this assumption because I believe that the cause for the periodic crystalli- 
zation is the same in all the previously mentioned cases, that is, that the 
cause is the same whether the crystallization is from solution or from bulk 


or whether it is in the drawn state. This assumption excludes kinetic 


causes. It should not be forgotten that during the stretching the mole- 
cules are arranged to a great extent in a parallel way from the exterior. 
If the molecules are given enough mobility by annealing, there results ex- 
tremely favorable kinetic conditions for the formation of large single crys- 
tals. However, this observation has not been made. 

The assumption of periodic crystallization is supported by similar behav- 
ior of the long periods in drawn material and of polyethylene single crystals 


- during annealing. In ricure 3 the dependence of the small-angle period 


of drawn 6-nylon fibers on temperature can be seen. These measurements 
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were carried out by K. Kobayashi (private communication). The long 
period increases with heating; the increase is particularly rapid at high 
temperatures. 

A curve similar to FIGURE 3 is obtained by heating polyethylene lamellae 
crystallized from solution. Data for FiguRE 4 were obtained by Statton 
and Geil!” and by G. F. Schmidt and E. W. Fischer (unpublished data). 
In the temperature range between 110°C. and the melting point, the fold 


Annealing temperature Kell 
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_ Ficure 4. Dependence of the long period of polyethylene crystals on the anneal- 
ing temperature. Key: @ = measured at room temperature by Statton and Geil!” 


(original erystals 104 A); O = measured at annealing temperature by G. F. Schmidt 
and E. W. Fischer (unpublished data; original crystals 120 A). 


period of the crystals changes, and the system comes to a new state of 
equilibrium. The points in rigurE 4 obtained at the annealing tem- 
perature show that the new period had already formed at these tempera- 
tures before the sample was cooled. It is obvious that smaller crystals 
must melt at a lower melting point because of the large contribution of the 
surface energy to the whole energy per volume unit. This melting process 
could also be observed by means of the X-ray wide-angle scattering with 
cakes of polyethylene single crystals held at constant temperature (G. F. 
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Schmidt and E. W. Fischer, unpublished data). According to the kinetic 
theories, however, the crystals should grow infinitely because the kinetic 
hindrance of the prefolded surfaces has been removed. Nothing should 
prevent them from growing into their most stable thermodynamic config- 
uration; that is, they should grow as large as possible. This phenomenon 
was not observed; however, with a small difference of temperature, a new 
fold period that is slightly longer is observed. In my opinion this can be 
explained only by the fact that for high-polymer crystals there exists a 
stable length, which depends on the temperature. 

It seems useful to mention two other experiments. The first experiment 
concerned the dependence of the fold period on the chain length; the second 
concerned the structure of polyoxymethylene single crystals. It was 
shown by Keller and O’Connor”™ that the fold length of polyethylene crys- 
tals is independent of the molecule length even if there is a change of the 
molecular weight by a factor of about 100. It is important for the theory 
of periodic crystallization that this independence holds true even in the 
region where the molecular length becomes comparable with the fold 
length. Other experiments carried out by Zahn™ with very short mole- 
cules of the 6-nylon type gave similar results. At first the long period 
increased with the number of the monomer units n, but above n = 6 it 


~ remained constant at about 70 A. This value corresponded exactly to the 


y 


long period measured in high-molecular 6-nylon. I cannot believe that 
there is a kinetic reason for this behavior because kinetic theories start 
from the fact that the nucleus with folded chain configuration is formed 
by only one molecule. 

Finally let us consider the structure of polyoxymethylene. Recently it 
was shown by Geil et al.” that this polymer grows from dilute solution in 
the form of lamellae about 100 A thick; hence the molecules have a folded 
conformation. It is true that polyoxymethylene single crystals—however 
different—have been known for a long time.” This was the first and, 
to date, the only polymer that permitted the formation of crystals large 
enough to obtain rotation patterns for the determination of lattice. Ob- 
viously the production of crystals as large as 1 mm. was possible in this 
case only because the crystal growth coincided with the polycondensation. 
The building elements of these crystals were not previously formed chain 
molecules; instead the growth took place in such a manner that formalde- 
hyde molecules were added to a nucleus. Nevertheless, according to the 
thermodynamic theory, the formation of single crystals consisting of chain 
molecules with such large macroscopic dimensions is highly improbable 
since the minimum of energy is within the range of about 100 A. On the 
other hand, of course, kinetic theories concerning macromolecules cannot 
be applied to these crystals. Therefore it was important to study the 
possibility that these crystals might have a microstructure. 

Indeed electron micrographs of surface replicas of macroscopic polyoxy- 
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methylene single crystals* show that the crystals are built up of many 
growth layers arranged extremely regularly (FIGURE 5). The period in 
chain direction amounts to approximately 200 A; consequently it is exactly 
within the expected order of magnitude. There is no doubt that the macro- 
molecules are longer than this period, and I believe that this experiment also 
favors a thermodynamic interpretation of the periodic crystallization. 


Nitin 5. Electron-microscopie surface replica of a polyoxymethylene single 


Basic Idea of the Theory 


What feature of a lattice built by chain molecules can cause the energy 
density of a erystal to be size-dependent? It must be a very general prop- 
erty of such lattices, and the main difference between polymer and low- 


molecular crystals is the extreme anisotropy of lattice forces in the poly-— 


mers. In chain direction the forces are diamondlike; perpendicular to 


*I thank O. Schweitzer, Frankfurt/Mainz, Germany, for preparation of the , 


sample. 
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chain. direction the forces are considerably weaker, but those forces be- 
tween neighboring chains preventing rotation around or translation along 
the chain axis are particularly weak. Stockmayer and Hecht,” for exam- 
ple, who calculated the heat capacity of chain crystals, proposed a model 
in which the forces perpendicular to the chain amount to only one twenty- 
fifth of the forces in the chain, and the force opposing the translation in 
chain direction amounts only to two thousandths of the force between ad- 
jacent groups in the chain. 

What are the consequences of this anisotropy of the lattice forces? If 
distortions are introduced into an ideal lattice, the position of a mass point 
S, in first approximation, given only by the position of the immediately 


Figure 6. Lattice distortions in high-molecular and low-molecular lattices. 


neighboring mass points within the molecular chain. The lattice forces 
are not strong enough to cause the fitting of the mass points into the ideal 
positions in the lattice. On the contrary the minimum of the potential 
energy of the lattice with an included defect is reached if the chain bonds 
and angles are only slightly deformed. This fact is demonstrated in FIGURE 
6 by the example of a branching structure. Although this pattern is very 


- much simplified it shows that in contrast to the case of low-molecular lat- 


tices such lattice defects do not heal easily with increasing distance. This 
one-dimensional defect is a case similar to that of dislocation. Therefore 


. the density of the defects per unit volume and also the defect energy de- 


‘pend not only on the number of the branchings in the crystals but also on 


the size of the crystals in chain direction. It should be recalled that in 
normal lattices the energy of a linear dislocation is infinitely large for an 
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infinitely extended crystal. In such crystals only closed dislocation rings 
have a finite energy. 

In high-polymer crystals distortions of this kind may be caused by vari- 
ous phenomena such as branchings, hydrogen bonds wrongly fitted, and 
distorted or entangled chains. It must be expected that in a crystal such 
as that shown in FrIcuRE 6 the long-range order in chain direction breaks 
down by introduction of further branchings. Therefore the crystal is no 
longer stable after it has exceeded a certain length. In agreement with 
this concept electron microscopy and X-ray small-angle scattering show 
that the sizes of the crystallites diminish with an increasing degree of 
branching” or with an increasing number of lattice failures caused by poor 
crystallization conditions. 


Lattice Vibrations 


It can be demonstrated that in an ideal chain lattice “long-range”’ dis- 
24,25 


tortions can also be initiated through heat motion. For this purpose 


Fiaure 7. Model of a linear chain within a longitudinal lattice field. 


FIGURE 7 represents a simple one-dimensional model. The groups of the 
molecular chain are represented by N mass points with the mass m, linked 
to one another by spring forces f. The simplifying assumption is made 
that the chains in the crystal fluctuate almost independently and that the 
force field due to the neighboring chains can be represented by a cos-poten- 
tial function: 


Q7rz 


go — $ COS =A (1) 


The meaning of the potential amplitude can be seen in FIGURE 7; c is the 
repeating unit of the chain. If we assume that ¢= 0, it can easily be 


shown that the mean square deviation of the free chain with its center of 
mass fixed is given by 


— N (2) 


Thus the mean square deviation of chain points from their equilibrium 
positions increases linearly with the number of the links of chain. 
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The energy of the chain embedded in the lattice field is the sum of the 
kinetic energy plus the potential energy of the springlike forces within the 
chain itself and of the intermolecular cosine field. The result is: 


Mo 2 fy 2 = 212; 

Br 5 Dyer +5 (27 = 251) + Ndo — ¢ >) cos——? (3) 
j=1 j=1 j=1 (6 

Using this setup Peterlin and Fischer” calculated the state sum of the sys- 

tem in a relatively comprehensive calculation. From the state sum the 

mean square fluctuation of the chain links in the lattice can be derived. 

Equation 4 is obtained: 


2 _ NkT Hs 0.4 
seer oaG (1 eth) 4) 


with 6= $¢/kT and A” = NkT/fc’. It is apparent that the mean square 
fluctuation in the lattice is smaller than in the free chain (2); however, it 
increases with increasing N. This dependence is a specific property of 
polymer crystals with large anisotropy of lattice forces. Although the 


_ chains are prevented from making translation in the plane perpendicular 


to the chains, the groups can still move rather easily in chain direction. 
The same applies to the motion around the chain axis, as shown later. 

In a former publication on this subject™ I was of the opinion that this 
dependence of the fluctuation amplitude on N led to an increase of the 


energy density with increasing crystal thickness. F.C. Frank and M. H. 


L. Pryce have called my attention to the fact that this is not correct if 


the potential amplitude ¢ is kept constant. This changes however, if it is 


taken into account that the potential is smeared by longitudinal oscilla- 
tions of the neighboring groups. It is assumed that ¢= ¢o for the case 
in which all neighboring groups are in state of rest. However, if all four 
neighbors oscillate independently of one another, it can easily be shown 
that: 


@ = oo exp (—2n's,/c’) (5) 


That means that the potential amplitude ¢ decreases with increasing mean 


square fluctuation s,” or with increasing number of chain groups in the 


chain direction. If the free-energy density is calculated from the state 
sum, taking into account the contribution of the surface energy, approxi- 
mately the following expression for the energy density is obtained. 


ee 1, ob 
NkT = const + N a kT os 
2 5 (0.7 + 0.5 In N(1 + 4NBA®) + 20°0/KT| 


In Equation 6 the terms independent of NV are included in the constant; 
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b? means the size of the boundary surface per chain. In this expression 
for the free-energy density, the potential function ¢ contributes a term in- 
creasing with N as a consequence of the smearing-out effect. The term in 
brackets diminishes with increasing N, the term with the boundary surface 
energy being the most important one within the brackets. Thus a term 
increasing with N is obtained as well as another that decreases with N. 
Both terms together yield a minimum of free energy with a certain N, which 
depends mainly on the temperature and the intermolecular forces. <A plot 
of this minimum Nin representing the thermodynamically stable crystal 
thickness as a function of 2b’o/kT with Bo = ¢o/kT as parameter is shown 
in FIGURE 8. The higher fp is, the stronger the interchain forces and the 
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Fiaure 8. Stable crystal thickness Nmin as function of surface ene i 
Bo = ¢0/kT as parameter. oe 


smaller the stable crystal thickness; this is in qualitative agreement with 
experimental data obtained from different polymers. 

With increasing temperature the crystal expands chiefly in the direction | 
perpendicular to the chain axis; consequently the interchain potential en- 
ergy Bo decreases. Furthermore the number of thermic failures increases 
and this also leads to a decrease of By). Since the changes in 7’ are anatl 
compared with the absolute temperature 7’ itself, there is only a small 
change in the abscissa. On the other hand, by using-a smaller parameter 
in the plot an increased thickness of the thermodynamically stable crystal 
is obtained. The thickness increases more quickly the smaller the original 
Bo , that is, with increasing temperature the thickness increases faster the 
higher the temperature; this is in good qualitative agreement with the ex- 
perimental results of Statton and Geil.’ 

The values of the parameters received from the theory are also in the 
right order of magnitude. I started from a boundary surface energy of 
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o = 50 erg/cm.’, a value also estimated by Price’® from the heat of fusion. 


_ If a value of the surface energy is chosen then the magnitude of the lattice 


potential amplitude ¢o can be obtained from the theory by means of the 
stable crystal thickness given experimentally. The result is 


(0) crystal thickness = NEN ee erg. 


On the other hand, an estimation for ¢p can be given with the aid of the 


_ well-known setup for the van der Waals potential between two CH» groups 


of different molecules. From the values for the van der Waals constants 
found in the literature” ¢ can be calculated by summation of the inter- 
action potential for the four next neighbors, giving: (d0)yan der Waals = 
0.3 X 10 “erg. This agreement, although not very good, is satisfactory 
considering the uncertainty of the constants. 


Torsional Vibrations 


A similar dependence of the free-energy density on the chain length is ob- 


_ tained if the torsional vibrations of the chain are considered instead of the 


longitudinal vibrations. Again the chain is placed into a periodic lattice 


potential, and the following expression for the energy of the chain in ques- 


tion is obtained. 


N-1 


ae N 
(9; — gin) + NE, — E ay, cos 29; 
sl j=1 (7) 


F N 
u= dio t+ 
2 jl 2 sin’ 5° 


- ; means the displacement of the jth chain link from its position of rest 
- in radian measure. Instead of the spring constant f the expression 


c/2 sin’(0/2) 


_ is obtained; according to Szigeti’’ this amounts to 0.64 X 10 ” erg/group 
~ when considering torsional vibrations of CH: groups. In this case the angle 
- @ is the C—C— bond angle, the constant angle between three C atoms. 
Contrary to the calculation of the longitudinal vibrations, this case gives 


10 “erg. It can also be calculated from the heats of transition of paraffins. 


If these values are substituted into the equations for the free energy of the 
vibrating chain, the curve in ricuRE 9 is obtained. In this case too the 
~ minimum is placed near N = 100 if the temperature T = 373° K. A cor- 
_ responding boundary face energy of 135 erg/cm.’ is obtained, a value that 


— 
Z 
4 
a good approximations for Hy). According to Eicker”® Ep is equal to 1.75 X 
A 
/ 
4 
E 
: 


_ is slightly too large. This is a consequence of the inaccuracy of the setup 
for the smearing of the potential. .The mutual coupling of the chains 


_ markedly undergoes this effect so that the smearing actually appears less 
than calculated. With a smaller smearing-out effect, the theory also gives 
~asmaller value for surface energy. Thus the incorrect value of 135 erg/cm. 
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indicates that the torsional vibrations of the chain are not independent of 
those of the neighboring-chains. 
Branching and Other Lattice Failures 


A consideration of the influences of branching and other lattice failures 
is in order here; therefore the scheme of a lattice with branchings shown in 


F, 
KR 
NkT 
0.3 
0.2 
Ol 
T=100°C. 
E= 18x!0erg 
fo) 


200 300 400 500 N 


Fiaure 9. Free-energy density as function of crystal length for torsional vibra- 
tions, 7 = 100°C.; Hy = 1.75 X 10-“ erg. 


FIGURE 6 should be reviewed. Obviously each of those failures causes a 
smearing of the potential in the same way as that shown in the case of 
longitudinal vibrations and torsional vibrations. In accordance with this 
observation Hendus et al.” found the dependence of the long period on the 
degree of branching plotted in TaBLe 1. There are good reasons to 


TABLE 1 
Lone Prriops or Vartous PoLYETHYLENES* 


Bee es iveackine Slowly cooled Quenched 
CH3/1000° C. 
1st period 2nd period ist period 2nd period 
1.5 344 165 222 86 
5 310 140 204 78 
15 280 85 184 60 
25 218 _— 154 — 


*From Hendus et al.2® Reproduced by permission of Ergebnisse der exakten 
Naturwissenschaften. 
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assume that the second period gives the size of the crystals in ¢ direction. 
The size decreases with increasing degree of branching as would be expected 
from the model in FIGURE 6. 


Conclusions 


The mean amplitude of fluctuation of the thermal motion of groups in 
chain crystals depends on the thickness of the crystals in chain direction. 
These distortions cause a smearing of the lattice potential. Therefore the 
energy density of a crystal consisting of chain molecules increases with in- 
creasing length of the chains. It seems that this dependence is an appro- 
priate explanation for the periodic crystallization observed in polymers. 
Having passed its critical length the nucleus of the crystal then continues 
to grow in chain direction until the minimum of free energy is attained. 
Then the chains fold, and the loops cause the decoupling of the longitudinal 
vibrations of subsequent parts of the chain. Another cause of this de- 
coupling is to be found in the amorphous regions of drawn fibers. 

Of course my model is not fully satisfactory in spite of qualitative agree- 
ment with the experimental results. It will be necessary to examine closely 
the chain motion in polymer crystals. Nuclear magnetic resonance 
(NMR) investigations form a suitable method for this purpose. These 
investigations have already shown that the chain mobility in the crystalline 
regions is greater than previously assumed. From their measurements 
McCall and Slichter” derived a value of the mean torsion amplitude very 
similar to the one I had assumed. Furthermore the NMR studies of Hug- 


gins et al.” force the conclusion that translational motion of appreciable 


magnitude occurs in polymer crystals below the crystal melting tempera- 
ture. 

Such investigations are suitable for reexamining the theory in a quanti- 
tative way, but I believe that the simple model gives a possible explana- 
tion for periodic crystallization. 
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LATTICE THEORY OF CHAIN POLYMER SOLUTIONS 


Michio Kurata* 
Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Mass. 


Introduction 


It is now well known that the Flory entropy of mixing for chain polymer 
solutions’ is straightforwardly derivable from the phase integral of the sys- 
tem without the aid of a special model such as the lattice model, if the cor- 
relation between segments arising from the chain connectivity is entirely 
ignored. This implies that the Flory entropy, although crude in itself, 
can be given a fundamental role in a more general polymer solution theory. 
Furthermore, the thermodynamic properties of systems characterized by 
- the Flory entropy have been fully investigated by many authors in connec- 
~ tion with phase separation and other practical problems." For these two 
- reasons, the Flory entropy may be adopted as the most convenient stand- 
~ ard for thermodynamic description of general polymer solutions, and the 
entropy of dilution S and the heat of dilution Ho of the system may be 
_ written in the following forms: 


S = —R{In(1 — 6) + [1 — (1/n)]¢ + o¢} (1) 


and 

Hy = RT Ke (2) 
where R is the gas constant, 7 is the absolute temperature, ¢ is the volume 
fraction of polymer, and n is the number of segments per molecule. More 
explicitly, n is defined as 
where M is the molecular weight, dis the (partial) specific volume of poly- 
mer and V» is the molar volume of solvent. The quantities o and x are 


respectively an entropy and a heat function, defined by 1 and 2 in a purely 
phenomenological sense. The chemical potential of the solvent is then 


written in the form 


uo = wo + RT{In(1 — 6) + [1 — (1/n)]6 + x) (4) 


with 
X=ote (5) 


and the chemical potential of the polymer is obtained from po with the 
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Gibbs-Duhem equation. Here yo* represents the chemical potential of 
pure solvent. 

At infinite dilution, the osmotic pressure derived from 4 should tend to 
the van’t Hoff limiting law. This condition shows that both o and « ap- 
proach a finite value at zero concentration. At reasonably low concen- 
trations, therefore, o, x, and x can be expanded in series form: 


o = Dod’ = oo tap + ox + °°: (6) 

r= Dinh’ = co + md + mop +: (7) 
and 

x = Dix’ = x0 + xb + xh + > (8) 


and the following expressions can be obtained for the second and the third 
osmotic virial coefficients: 


A= (P/V) (% — a0 — ko) (9) 
A; = (8°/Vo) (1% — a1 — m) (10) 


It is now well established that the second virial coefficient A» is a de- 
creasing function of molecular weight. This indicates that in this forma- 
lism 9, the oo and xo should depend on the molecular weight in contrast to 
the prediction of the existing lattice theories." It is also obvious that 
these functions should depend on the structural properties of polymer 
chains, such as flexibility and branching. 

The theta temperature where the second virial coefficient vanishes may 
be defined by 


14 — oo 0, M) = xo(O, M) (11) 


where oo(©, M) indicates the value of oo(7, M) at T = @,and soon. Ac- 
cording to the recent experiment on polystyrene-cyclohexane system by 
McIntyre et al.,° the © value is found to increase slightly as the molecular 
weight is decreased. For the same system, Krigbaum and Geymer’ have 
observed that the third virial coefficient A; also vanishes at the theta tem- 
perature. A similar result has been obtained for polyisobutylene-benzene 
system by Flory and Daoust.” Thus I am almost sure that 


1(0, M) = 0(0,M) + m(@, M) = % (12) 


The establishment of the dilute polymer solution theory in this decade is 
due largely to the experimental data obtained at the theta temperature or 
in its vicinity. For the same reason, the two facts mentioned above, that 
is, the molecular weight dependence of the @ value in the region of rather 
low molecular weight and the simultaneous vanishing of A, and A; at O, 
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will furnish a key for the determination of the functional form of ¢ and x 
as functions of 7, M, and ¢. 

Concerning the thermodynamic behavior in the more concentrated range, 
I refer to the interesting analysis by Zimm” and Zimm and Lundberg.” 
Based on the general solution theory of Kirkwood and Buff,” they derived 
an equation relating the pair distribution function F2(7, 7) to the activity 
of solvent a). This is 


Cale ) ao uke 
f= (5)- a 


with the definition 


Goo = (1/V) ff We, 4) — Maca) aC) (14) 


_ Since F2(7, 7)d(2) d(7) gives the probability of finding a solvent molecule 7 


in the volume element d(7z) at the position (7) and another solvent molecule 


— jin d(j) at (7), the quantity (1 — $)Goo/Vo represents the mean number of 


solvent molecules in excess of the random expectation in the neighborhood 


of agiven solvent molecule. According to the analysis of the activity data 
- of rubber-benzene and polystyrene-toluene systems by Zimm and Lund- 
berg, the value of Goo/Vo exceeds unity over the whole range of concentra- 


tion, indicating a strong clustering tendency. This observation is impor- 
tant, because the combination of the Flory or Huggins-Guggenheim-Miller 


theory with 14 leads to 


eh rN 


‘el 
4 


Yh ee 


PRC? Eee 


Goo/ Vo 0 (Flory ) 


(15) 
2/(2 — 2) (Huggins-Guggenheim-Miller ) 


which can never exceed unity unless an unreasonably small value, for ex- 
ample, 3, is assigned for the lattice coordination number z. Accordingly, 
even in concentrated solutions the clustering effect of solvent molecules, 


- which is undoubtedly in the closest connection with the chain connectivity, 
can no longer be neglected. 


In this paper I shall develop a lattice model theory of concentrated poly- 


mer solutions and investigate the relation between the thermodynamic 


properties and the chain connectivity in the light of the presented theory. 
From a purely statistical mechanical point of view, there are at least three 
fundamental problems to be considered in the lattice model treatment of 
liquids or solutions: (1) the communal entropy, (2) the free volume, and 
(3) the combinatory problems.””’ The first problem is, of course, very 


_ important for detailed investigation of properties of pure liquids, but gen- 


erally not for rather rough treatment of solution properties. Concerning 


the free volume problem, Prigogine™ and his co-workers have recently de- 


veloped an approximate treatment called the ‘average potential method” 
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and have applied it to liquid mixtures involving polymer solutions with con- 
siderable success. In the case of polymer solutions, however, the chief 
defect of existing theories seems to arise from the incomplete treatment of 
chain connectivity, as mentioned above. This is, of course, a typical con- 
figurational problem or, in the language of the lattice theory, the combina- 
tory problem. Thus in this paper I am principally concerned with the 
combinatory problem, but not with the free volume problem. 


Basic Formalism 


Polymer solution as a special type of regular assembly. Consider a polymer 
solution composed of No solvent molecules and of N; polymer molecules of 
the same molecular weight M. Each polymer molecule consists of n 


© 


@ 
Q @ © © 
© © © @ 
@© @—Z @) 
© O-2 © 
© © © 
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Freure 1. Schematic representation of polymer solutions. 


identical segments, which are linearly connected by n — 1 chemical bonds. 
A segment occupies one lattice point as well as a solvent molecule does. 
Later, I shall need to distinguish between two kinds of segment contact, 
that is, the intramolecular contact that is formed between two segments 
belonging to the same polymer molecule and the intermolecular contact. 
For this purpose, let us label the same index number 7 for all segments of the 
jth polymer molecule. In Figure 1 is presented a schematic illustration 
of the system. From rigure 1 it will be easily understood that the given 
system is nothing other than a special type of the so-called regular assem- 
bles. In other words, the system is a regular solution composed of No 
molecules of type 0, n molecules of type 1, n molecules of type 2, --+ , and 
n molecules of type N;, all having the same size. The number of com- 
ponents is thus counted as N; + 1. There is, of course, a remarkable speci- 
ficity in the present system that distinguishes it from the regular solution 
of ordinary sense. This is the strong tendency to clustering of molecules 
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of the same kind or, in the usual language of the regular assembly statistics, 
the existence of an order of extremely high degree. This kind of peculiarity, 
however, will not obstruct the formal application of regular assembly theory 
to the present system in a decisive way. 


Short sketch of the regular solution theory. According to the pseudocon- 
figuration method that was first proposed for treating the Ising ferromagnet 
by Kikuchi” and Kurata et al.,"* the combinatory factor G of the pair ap- 


_ proximation is formally written as 


G = [N!/TI (x,N)!]-(Q/P) (16) 
with 
Q = 4zN1/ II et (yinbazN )! (17) 


for the regular assembly of multiple components. Here z is the coordina- 


~ tion number of the lattice, x; 1s the probability of finding a molecule of type 
jon a given lattice point, and y; is the probability of finding a molecule of 


\ 


MR 


type j on a given lattice point and another molecule of type / simultane- 
ously on another given point adjacent to the former. By this definition 


- Equation 18 is obtained: 


Lj = es Yik (18) 


The meaning of 16 is as follows. The first factor, N!/]]; (z;N)!, is the 
combinatory factor of the point approximation, which obviously represents 


~ the number of distinguishable ways of distributing a given set of molecules, 


for example, z,N molecules of type j and N molecules in total, on the lat- 


See oy te 


tice, irrespective of the number of molecular pairs. To obtain the combina- 


tory factor corresponding to a fixed number of molecular pairs, it is neces- 


‘sary first to distribute the given set of molecular pairs, YixvQ2N pairs of 
type j = k and éN pairs in total, on the lattice. This is done in Q dis- 
 tinguishable ways. However, the pairs are here independently distributed 
on the lattice; accordingly, these Q configurations include a large number of 


wrong configurations, in the sense that a certain lattice point is simultane- 
ously occupied by more than one molecule. Ficur® 2a illustrates such 


“wrong configurations in a one-dimensional assembly of two components, 


while FIGURE 20 illustrates a right configuration. ‘Thus the combinatory 
factor G is obtained by picking out the right configurations from Q or mul- 


tiplying Q by a selection factor T. 


To evaluate I, we define the pseudoconfiguration P as the configuration 


consisting of the same number of molecules as that involved in the parent 


configuration Q. These are x;2N molecules of type 7 (Equation: 18) and 
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2N molecules in total. Thus the number of pseudoconfigurations P is 


readily counted as 


P = 2N1/TI (z2N)! (19) 


These P configurations also include a large number of wrong configurations 
of a nature similar to those included in Q in addition to the right configura- 
tions, whose number is obviously given by N!/ II; (z;N)!. These are il- 
lustrated on FIGURE 2c and d. Thus if we use the selection factor of the 


e ° t e °o ce) ie) ° ® @ ° 
ee 
e @ ° e ° e ° e ° e 


(c) "wrong" configuration of P 
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(d) "right" configuration of P 


Figure 2. Parent configuration Q and its pseudoconfiguration P. 


right configurations from P instead of I’ we can obtain 16 for the combina- | 
tory factor G. This approximation has been proved to be correct for such 
a singly connected lattice as the one-dimensional lattice or the dendritic 
lattice shown in FIGURE 3.7° 


The potential energy of the system is easily written as 


E= > (yinbG2N ) & jx (20) 


where &j, is the energy of each 7 = k contact. This energy F is then re- 
written in the form 


Ho >. 1 eN Ej; + 2, Y jn2N AS jx (21) 
J j 
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with the aid of 18 and the definition that 


AS jn = jx — 44( 855 + xx) (22) 


The most probable value of each y jz is, of course, determined by minimiz- 
ing the free energy of the system F = EH — kT In G with respect to all in- 
dependent: variables, for example, all y;,s but not yii8. This minimization 
yields 


Vix/Ysiyex = exp(—2AEy,/kT) (23) 


which is usually called the equation of quasi-chemical equilibrium. Here 
k is the Boltzmann constant. If A&js are all negligibly small compared 


Fieure 3. A dendritic lattice (g = 4). Reproduced by permission of The Journal 
of Chemical Physics .38 


to kT, the following solution of 23 is obtained: 
Yik = Ujlk (24) 


_ which corresponds to the case of ideal or perfect solution. In general, 23 


leads to a more complicated expression for y;,, which is not reproduced 
here. Nevertheless, in the case of regular solutions, the most probable 
value of y; for solute-solute contacts tends to zero at infinitely high dilu- 
tions (x — 1). This fact, while rather obvious, is important, especially 


in comparison with the following case of polymer solutions. 


Combinatory problem for polymer solutions. ‘Turning to the polymer solu- 


_tion defined above as a special type of regular assembly, the singlet proba- 


bility 2; can be written as 
zo = No/(No + 2M1) = No/N, (25) 
x; = n/N for =-1,2,-°: ,.N; 


It is also obvious that each polymer molecule involves many intramolecular 
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segment contacts: at least n — 1 corresponding to chemical bonds and, gen- 
erally, more. Thus it may be stated 


yiihgeN = (n—-1) +C€ forj = 1,2, ---, Wa, (26) 


where C denotes the number of intramolecular segmental contacts in excess 
of the chemical bonds. 

Concerning this quantity C, three properties are to be noted: 

(1) C must have a very close relation to the average statistical dimension 
of polymer chain, which plays an essential role in the dilute solution theory. 
The value of C is naturally expected to vary with expansion of chains. 

(2) The number of this kind of segmental contacts never tends to zero 
even at infinite dilution but reaches a finite value, depending on the nature 
of the given chain, such as flexibility or branching. In other words, to de- 
termine the equilibrium value of C’,, an equation other than the quasi-chemi- 
cal equilibrium Equation 23 is needed. In this connection it is to be noted 
that in the case of rigid polymer solutions C is definitely a structural con- 
stant that depends on the geometrical shape or the compactness of the 
molecule. These facts indicate that in the present formalism the quantity 
y;; for each molecule should be regarded as a structural quantity rather 
than a variable. 

(3) The existence of C in excess of n — 1 chemical bonds generally de- 
stroys the single connectivity of the chain; in other words, the chain cannot 
be separated into two fragments by removal of one of the n — 1 + C bonds 
unless C = 0. Since the pseudoconfiguration method, 19, is rigorously 
applicable to singly connected systems alone, these C excess bonds must be 
excluded from the pseudoconfigurational treatment. 

Keeping the above three comments in mind, let us proceed to calculate — 
the combinatory factor G of the polymer solution. The number of dis- 
tinguishable ways of distributing a given set of molecular or segmental con- 
tacts on the lattice is straightforwardly written from 17 and 26 in the form 


O% YozN ! 
Ni Ny Ny : 
(yo42N)1 TT (yoid42N)! TT TT (inde) 1Cu(n — yyy (27) 
jxk 


Here the chemical bonds were treated separately from the excess intra- 
molecular segmental contacts. Because of the identity of the N, poly- 
mer molecules Equation 28 can be written: 


Yor= Yr and Yr = Yu (28) 
and then Equation 27 rewritten into a more simple form: 
YozN ! 
Q= Za (29) 


(yoo'%eN ) !(yont4eN ) it (yn%eN ) pa IC I(n ra 1) hee 


ey Sa tie $e Wi act 


Aa a aaa 


LR NT RY Nh RT RS TS 
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As for the mutual relation of variables, Equation 30 is obtained from 18, 
25, and 26 


Dyn = 23 — Y53 = (No/N) = Yoo ifj = 0, 


kAj 


(30) 
= q/N Wie Gms Uh Deetin s. N 
with the abbreviation that 
g=n — (Q/z)(n —A + C) (31) 
Then, using 28, we get 
Yoo = (No/N) — Niyo , 
(32) 


(Ni — lyn = (g/N) — yo 


The number of molecules or fragments of the pairs used in the above Q 
configuration is as follows: zNo of the solvent molecules, zq of the fragments 
of the intermolecular segment contacts for each polymer molecule, 2(n — 1) 
of the fragments of the chemical bonds, and 2C of the fragments of the ex- 
cess intramolecular segment contacts. Thus applying the pseudoconfigura- 
tional method to the remaining part of Q except C term, we obtain 


L6zN! (zN — 2C)! 


is (C!1)"*(Y4eN — C)! 2No!(zq!)™[2(m — 1)!]" (33) 
then, substituting 29, 33, and 25 into 16, we obtain 
N! 2No! (zq!)™? 

© Nol(nty WENN HEN (34) 


_ OeeN =O) 2(m — 1) 1 1 
(n — 1)!" (2N — 2C)! M! 


where the last factor 1/N,! represents, of course, the indistinguishability of 
N, polymer molecules. 
The energy of the system is obviously given by 


E = yoobzN &00 + Z (Yor + Yiv)2N &o; 


Ni M1 (35) 
ar ape Y ne VO2N & jx = N,C&;; 
ee 
oe 
Because of the identity of all segments, we can simply put 
&0; =) Gop and Six = GF = &1 (36) 


and rewrite the 35 as 


ee 14zN 800 ote [Leen ae (n = 1)|Mi8u + you2N Ni AE j (37) ; 
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with 
AS = &: — 14(&00 + &u) (38) 
The Equations 34 and 38 give the expression for the free energy of the 


system in terms of an independent variable ya. Thus, minimizing the free 
energy with respect to yo , we again find the equation of quasi-chemical 
equilibrium, 

yon/yoyn = exp(—248/kT) (39) 
which indicates, as usual, that the larger value of A& is followed by the 
higher clustering tendency of the same kind of molecules. This kind of 
clustering tendency, however, is of minor importance in the polymer solu- 
tion theory when compared with another extremely strong clustering tend- 
ency due to the chain connectivity. In other words, the adoption of the 
approximate solution of 39 corresponding to the Bragg-Williams approxi- 
mation may be acceptable as a good approximation. Thus we obtain from 
39 and 32: 


Yo = Nc /NNg; yn = Nog/NN,; and yn = g/NN, (40) 
with the abbreviation; 

No = No + Qi (41) 
Then, substituting these ys into 34 and 37, we obtain 


qa NL Den = lt ONG ae 
Nol(n!)™ (YazN = CN)! N;! 
and 


E = WeNobw + [Yen — (n — 1)INi8u + (2NogMi/Nq)AS (48) 


As is easily seen, the first two terms in E represent the energies of pure sol- 
vent and polymer. The configurational entropy and the energy of mixing 
of the system are then written in the form 


S = —k{Noln(1 — ¢) + Niln} 
— [MeN — (n— 14 C)NJnf[l — 2271 — 2 + Cn™)¢) 


— Ni(n — 1) Ine "(1 — n *)d] (44) 
+ (44zN — CN) In{l — 2C(zn)¢}} 
and 
B = (Nae) ha 1 = n+ Cn 101 — @) (45) 


1 — 227(1 — n+ Cn“)¢ 
where ¢ is the volume fraction of polymer, 


= nN,/(No + nN) (46) 
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The parameter C is generally expected to be a function of ¢; therefore 
the entropy and the heat of dilution are respectively obtained from 44 
and 45 as 


So = —R Ee Se (re Le . 


do 
=} =1 ‘es (47) 
lars (1 —n* + Cn 2) 
1 — 221Cn“o 
and 
(ies hie we On yi 
+227 — 9) EE (ag) 


[IT — 224(1 — nt + Cn) oP 


ignoring the difference between HE and H. These two are the basic equa- 
tions of the present theory. Here N.« represents the Avogadro number. 
As is easily seen, these equations reduce to the Huggins-Miller-Guggenheim 
equations when C' = 0. 


Systems Involving a Constant Number of Intramolecular 
Segmental Contacts 
Lattice solutions of rigid polymers. According to the definition, the quan- 
tity C must be a simple geometric constant in the case of rigid polymers. 
In such a case, therefore, we can always neglect the concentration depend- 
ence of C and obtain 


a 1 — 227*(1 — n+ Cn')¢ 
Ses Ee — ) + We in( 1 — 2z7Cn—o a 


fi — 22711 —w + Cn) Pe" 
[1 — 224(1 — m1 + Cn™)oP 


Accordingly, the functions o and « defined by 1 and 2 are written as 
1 1 — 22"Cn “¢ ) ( ‘) | 

See ie! nC a 1 

B 2 E i ¢ — 221(1 — n+ Cn“) : n Salen bt) 


. ogi =n Cn) ) 
K=\"pp JM — 2°00 — 1 + Cn) 


respectively. Hence we have 


= Se ST eee ae (53) 
Al n n n 


aM) 


Hy = (2N.A8) (50) 


and 
(52) 
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On the other hand, we already know the rigorous value of oo of several 
lattice chains, which were obtained by elementary geometric computation 
of the excluded volume of the chain.”””’ These are shown in TABLES 1 and 
2. The shape effect on oo is, in fact, dominant for both short and long 
chains. 

In the present theory this kind of shape effect is expressed in terms of C, 
which can be easily evaluated from a molecular figure. As shown in TABLES 
1 and 2, the values of o» calculated from 53 are generally in good agreement 
with the rigorous values, in spite of the simplicity of the equation. There 
are, however, some systematic deviations of the calculated values of oo from 
its rigorous values. As is easily seen from TABLES 1 and 2, the present 
theory often predicts an overestimated value for oo , and the overestimation 
becomes large for molecules of high compactness. Thus, for example, in 
the case of a cubical molecule of high molecular weight, the theory predicts a 
value, 0.833, for oo , which is considerably higher than the rigorous asymp- 
totic value, 0.500. A value of oo higher than 0.500 is also obtained for the 
third example in TABLE 2: that is, a ribbon molecule with three-segment 
width in the two-dimensional square lattice. Since the second virial coeffi- 
cient of athermal solutions is necessarily positive, such predictions of high 
value of oo are obviously wrong. In this meaning, the range of applica- 
tion of the present theory is restricted within the range that 0 < C/n < 
(z — 2)/4 or, more seriously, within the region of small value of C/n as 
0.2 or 0.5 depending on z. It is also obvious from 53 that the present theory 
cannot predict the shape effect of bent molecules unless they involve at least 
one intramolecular segmental contact. The third and the fifth examples 
in TABLE | and the fourth example in TABLE 2 illustrate this effect. To 
treat this kind of bent effect, it is necessary to construct the combinatory 
factor of the system in a more detailed manner. For such treatment the 
reader may refer to theories presented by Guggenheim and McGlashan” 
and by Kurata e¢ al." However, apart from the mathematical interest, 
the bent effect is practically small, as shown in the tables, and so this defect 
of the present theory is probably not serious for most practical purposes. 


Theta solvent systems. Another interesting example of a system with 
constant C is the so-called theta solvent system. According to the theory 
of excluded volume effect,”””* polymer molecules generally take a more or 
less expanded form in dilute solution, owing to the volume effect of the seg- 
ments. ‘This expansion effect, however, vanishes at the theta temperature, 
at which the second osmotic virial coefficient tends to zero, and the polymer 
molecules reduce to the so-called unperturbed Gaussian chain. It is also 
known that the volume effect generally diminishes with increasing concen- 
tration and most probably vanishes in the pure polymer state. These 
facts suggest that at the theta temperature the statistical dimension of the 
polymer chain is kept constant at its unperturbed dimension over the whole 
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TABLE 1 
RELATION BETWEEN THE HNTROPY PARAMETER op AND THE MOLECULAR SHAPE 
Short Chains 


Two-dimens. square lattice Simple cubic lattice 
Mol. shape n @ Fie are 
Eq. 53 Rigorous Eq. 53 Rigorous 
5—? 2 0 0.063 0.063 0.042 0.042 
eeey 3 0 0.111 0.111 0.074 0.074 
| 3 0 atid) Werle iis 0.074 0.079 
= 4 0 0.141 0.141 0.094 0.094 
| : 4 0 0.141 0.172 0.094 
: : j 4 1 0.235 0.219 0.156 0.156 
A ee 60-0 
ae 5 0 0.160 0.160 0.107 0.107 
2 ih. 5 1 0.240 07235 0.160 
3 eee 8 0 0.191 0.191 0.128 0.128 
ze | a Suey 15 ~ as 0.310 | 0.289 
e a, 12 0 0.210 0.210 0.140 0.140 
a Pt 12 9 — — 0.369 0.337 
TABLE 2 
RELATION BETWEEN THE ENTROPY PARAMETER op AND THE Mo.LecuLAR SHAPE 
Long Chains 
Two-dimens. square Simple cubic lattice 
C lattice zs = 4 2= 
Mol. shape (n = ©&) = 
Eq. 53 Rigorous Eq. 53 Rigorous 
O—0—-O— - ---- —-o—o—0 0 | 0.250 0.250 0.167 0.167 
ao, = “OO 
: : i ; ; A 0.500 0.4388 0.333 0.333 
o—do—0- - - 
gel 
Saag ae - 24| 0.583 | 0.472 | 0.388 | 0.388 
> - — ---9—O 
AN /\, ENO o| 0.250 | — | 0.167 | 0.208 
— - 
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range of concentration because of the lack of the excluded volume effect. 
Thus, at least in an approximate sense, theta solvent solutions may be char- 
acterized as systems of constant C. 

For the sake of convenience, let us now define the coiling probability g by 


C = Weng (55) 


The g gives the probability that, if a segment of a molecule is on a given 
lattice point, another segment of the same molecule will be found on a near- 
est-neighbor site of the former. Here we neglect a minor end effect. The 
g is, of course, a function of 7’, M, and 4, in general; but at T = 0, it would 
be a function of M alone. We denote it by go(Q, 1) or simply go . 

Using this probability g we can rewrite 51 and 52 into 


x(0,M,¢) = o(0,M,¢) + «(0,M,¢) 


= £ | me In (ra - 6| (56) 


422) (Laat = oY 
k @\1 — (221+ + go)o 


Then we find: 
¥o(O, M) = % — ap(O, M) = (1 — 22 — 2) (57) 
xo(9, M) = (2A8/kO)(1 — 227 — go)” (58) 


Hence remembering the definition of the theta temperature (11) we obtain 


_ 22s (2 — 2\ [1 — (z/z — 2)gol” 
6 = ES) See 


k 2 1 — 2(z2/z — 2)go (68) 


This equation indicates that the molecular weight dependence of @ is a 
second-order effect of go(@, M) and is generally small, in accordance with 
the experimental knowledge. 

For the corresponding monomer solution, we obtain from 11, 53, and 54 — 


© = 2zA&/k (60) 


by putting» = land C = 0. This implies that the © of very low molecu- 
lar substances is some 20 to 50 per cent (depending on the value of 2) 
above the © of very high polymers of the same chemical species. If the 
latter is near room temperature, about 300° K., this difference ranges from 
60° to 150°C. Thus even if 1/n is as small as 0.01, the end effect on @ 
value still counts 1°. In the case of real polymers the molecule usually has 
an end group that is different from intermediate ones, and this end segment 
effect may enhance the variation of © as the molecular weight is decreased.° 
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Now, substituting 59 into 56 for elimination of © and expanding 56 in 
the form of 8 we obtain 


x(0, M) = (22 — 34)/z + [(z — 4)/zlgo — go’ (61) 


x2(O, M) = (6z — 10)/z* + [(62 — 20)/2"Igo 
3% 2 , (62) 
+ [(342 — 12) /z]go — 2g0 
TABLE 3 
THEORETICAL VALUES OF THERMODYNAMIC QUANTITIES AT THE THETA TEMPERATURE 
z 80 vo xo x1 x2 x(d = 1) |x( @ = 1)/«o 
0 0.250 0.333 0.219 1.386 4.00 
0.05 0.200 0.330 0.227 1.482 4,94 
.4 0.10 0.150 (0.500) 0.323 0.227 1.559 6.25 
0.15 0.100 0.310 0.216 1.590 8.16 
0.20 0.050 0.293 0.193 told intial 
0 0.333 0.259 0.121 0.967 2.25 
0.10 0.233 0.288 0.158 1.115 3.11 
6 } 0.15 0.183 (0.500) 0.287 0.170 1.180 ant 
0.20 0.133 0.286 0.173 1.229 4.59 
0.25 0.083 0.280 0.169 1.243 5.76 
0 0.375 0.208 0.074 0.817 1.78 
0.10 0.275 0.248 0.116 0.953 Deol 
ge 0.15 0.225 (0.500) 0.261 0.133 1.018 2.78 
0.20 0.175 0.268 0.146 1.077 3.31 
0.25 0.125 0.271 0.152 1,122 4.00 
0 0.500 0 0 0.500 1.00 
0.10 0.400 0.090 0.013 0.605 1.23 
00 } 0.20 0.300 (0.500) 0.160 0.044 0.718 1.56 
0.30 0.200 0.210 0.081 0.836 2.04 
lL 0.40 0.100 0.240 0.112 0.944 2.78 
Complete vanishing of virial 0.500 0.333 0.250 00 — 
coefficients 


~ It has been established in the past several years that the polystyrene- 
cyclohexane system is well specified by the following parameter values say 


© 307° K; and Yo = 0.19 — 0.23 (63) 


TaBLE 3 shows the values of Yo, x1, and x2 calculated by 57, 61, and 62, 
together with the used values of z and go. From TABLE 3 the following 
combinations of parameter values may be adopted as the best choice for 


the polystyrene-cyclohexane system: 
2=A4 and go(®, M) = 0.02 — 0.07, 
z=6 and go(9, M) = 0.10 — 0.15, (64) 
a= 8 and go(9, M) = 0.15 — 0.20 
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It is to be noted here that the predicted values of (14 — x1), hence the third 
virial coefficient A; is in fact very small at theta temperature, although it 
does not completely vanish. 


LS 
PS-Cyclohexane 
(34°G) 


° fraction I Mn= 25,100 
e fraction IT Mn= 72,000 


— 
. 
—"s 


e fraction Il Mn=440,000 


X (0) 


O 
ie) 


Free energy function, 
Oo 
NI 


O 0.2 0.4 06 08 1.0 > 
Volume fraction, ¢ 


FiaureE 4. Free energy function x fo ly : c mi 
Pier ee x for polystyrene-cyclohexane system at 34° C. 


After the values of z and go are once determined, the values of o, «, and x 
can easily be calculated over the whole range of concentration, with the aid 
of 56 and 59. Ficures 4 and 5 show the results in comparison with the 
recent experimental data obtained by Krigbaum and Geymer.’ Agree- 
ment between theory and experiment is satisfactory over a wide range of 
concentration, but an appreciable discrepancy is observed in the concen- 


OR TE eT Ne TT ee TE Pe 
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trated region higher than 60 or 80 per cent of volume fraction, especially 
in the case of the heat function x. In this region, however, the system is 
probably in the glassy state or in a transition state, as pointed out by Krig- 
baum and Geymer. 


1.4 


PS-Cyclohexane 


rc (@) 
O 


Heat function , 


O 0.2 0.4 0.6 0.8 1.0 


Volume fraction, ¢ 
Ficure 5. Heat function x for polystyrene-cyclohexane system at 34° C. and its 
theoretical values. 


Nature of the Coiling Probability 


Short-range coiling and long-range coiling. Let us now investigate the 
nature of the coiling probability g(7, M,¢). As may be easily understood 
from FIGURE 6, the coiling of chain occurs in two different modes, the short- 
range mode and the long-range one. The short-range coiling corresponds 
to contacts between two neighbor segments, as shown in FIGURE 6a and b, 


~ and this may be determined chiefly by such short-range effects as the restric- 
tion on bond angle and the steric hindrance for rotation around bond axis. 
On the other hand, the long-range coiling corresponds to contacts between 
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far distant segments, as shown in FIGURE 6c, and this must be closely con- 
nected with the problem of excluded volume effect. Thus it may be ex- 
pected that these two modes of coiling are considerably different from each 
other in their dependence on M, T’, and ¢. 

Let us first consider the molecular weight dependence of these coiling 
probabilities. Denote by y’ the probability that a sequence of four seg- 
ments finds a coiled form or square as given in FIGURE6a. This probability 


(a) (b) 


(c) 


Ficure 6. Short-range and long-range modes of chain coiling. 


y’ is essentially identical with the probability of appearance of cis configura- _ 
tion in the usual czs-trans rotational problem and may be independent of 
the molecular weight of chain. As there are n — 3 four-segment sequences 
in a linear chain, we have 


loeng,’ = y'(n — 3) & y'n (65) 
Similarly, for the six-segment sequences shown in FIGURE 6b we have 

loeng.” = y"(n — 5) S y"n (66) 
and soon. Therefore, adding up all these contributions, we obtain 


Moeng. = (y' ty” + -++)n = yn (67) 
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Equation 67 indicates that the coiling probability corresponding to the 
short-range mode, g,, is approximately independent of the molecular 
weight. 

On the contrary, the coiling probability corresponding to the long-range 
mode, g;, can be calculated in the following way. Assume an equivalent 
sphere model, for instance, in which the polymer chain is replaced by a uni- 
form distribution of unconnected segments within a sphere whose dimension 
is chosen to give the correct radius of gyration. This is, of course, 


(S’) = na’/6 (68) 


where a is the bond length or, more precisely, an effective bond length, which 
is determined by the short-range interferences. As the volume of the 
equivalent sphere is 


: V* = (49/3) (5/18)? (na)* (69) 
the number of segment contacts in this sphere can be obtained as 
Legeng, = rn’/V* = (34m) (18/5)°"(r/a')n\? (70) 


where 7 represents the volume of an effective segment. This calculation 
can be refined by using a segment distribution of Gaussian form instead of 
_ the uniform distribution” or by using an equivalent ellipsoid instead of a 
sphere, or both.” Such modifications, however, do not alter the propor- 
tionality of g: ton *”. 
According to the theory of random chain statistics,”’” the probability of 
- finding a segment, such as 7, in a volume element dv neighboring another 
- given segment, such as 7, is written as 


4 P(O;;)dv = (34n)*?a *(j — i)? dv (71) 


The total number of segment contacts in a molecule is then approximately 
evaluated by integrating 71 over all possible pairs of segments, that is, 
eng = av-[ [ POs) adidj 
2 vaeng (ae cells 1 t+2 ; 


(72) 
"n— 2 
= (z- 2)v9:(34na")* | an 


where v is the volume of a lattice point (or the volume of a solvent mole- 
cule). Hence we have 


(ASG » 


- The first term in the last parentheses represents the short-range mode of 
coiling, while the second term proportional to n represents the long-range 
‘one. Accordingly, we may conclude that the method of coupling of two 


TR Te ee PN 
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modes of coiling, gs and gz, is not addition but subtraction, in contrast to 
simple presumption. This indicates that the segment contact counted as 
the long-range mode, 70, isa correction for an overevaluation of the contact 
number in the short-range-type calculation, 67. 

In the presence of the excluded volume effect the polymer chain takes a 
more or less expanded form from its unperturbed dimension, depending on 
temperature and concentration. We define a linear expansion factor a by 


of = (L*)/{L*)y = (L*)/nae (74) 


where (L’) represents the mean square end-to-end distance of real molecule. 
If a uniform expansion of the polymer coil is assumed for simplicity, we 
obtain 


lézngi = const X (v9/a°)a on™” (75) 


instead of 70 for the long-range mode of coiling. On the other hand, the 
short-range mode of coiling is essentially independent of this kind of chain 
expansion. Accordingly, we obtain 


g(T, M,¢) = c(z — 2/z)(v0/a*)[1 — (¢’/ain'*)] (76) 


instead of 73. Here cand c’ are numerical constants. The coiling proba- 
bility g(T, M, ¢) is thus expected to increase as n and a are increased. 


Second virial coefficient, A... Substituting 57 and 58 into 9 and neglecting 
a term of order of g’ we obtain 


ee = ) 


where @ is the specific volume of polymer and Vo is the molar volume of sol- 
vent. As the coiling probability g is an increasing function of n, as men- 
tioned above, the second virial coefficient A» is decreased as the molecular 
weight of polymer is increased. Similarly, the quantity A./(1 — T”'®) 
is also a decreasing function of a. This behavior of 77 is favorably com- 
pared with that of Equation 78, which is obtained by the distribution func- 
tion theory of dilute solutions: 


me a h(z) (1 = 7 (78) 

with 
h(z) = 1 — 2.8652 + 18.512" — --- (79) 
z= (3¢r)""(r/a’)[1 — (O/T) In” (80) 


Here 7 and mp represent the volume and the molar weight of effective 
segment, and N4 is the Avogadro number. In this formalism (78), the 


ee eC 


ap \) 


A a8 


> | 
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function h(z) is unity at theta temperature, which monotonously decreases 
with increasing temperature, or a(z). 


Zimm’s clustering function, Goo. As mentioned earlier, the coiling prob- 
ability g consists of two parts, and the long-range part diminishes with in- 
creasing n. In a crude approximation, therefore, we may neglect the long- 
range mode of coiling and assume the quantity g to be a constant parameter 
independent of 7’, M, and ¢. 

With this assumption we can easily calculate Zimm’s clustering function 
Go of an athermal solution” by using 13 and 56. The result is 

Goo 2 + 229 — g(2 + 29)¢ 


Vo 2-2-2 + 92+ zg) ay 


which yields 


6-0 Vo 2— 2 — 2zag 


(82) 


at infinite dilution. The substitution of the parameter values, z = 6 and 
g = 0.15, into 81 predicts the following values for Goo/Vo : 


Go/Vo = 1.72 (at ¢d = 0); 1.48 (at¢@ = 1) (83) 


- The quantity Goo/Vo is slightly decreased with increasing concentration, 
but still kept above unity even at ¢ = 1. This result is in good accordance 
with Zimm and Lundberg’s observation.” 
For a more complete description of the solution properties we shall need 
- to investigate the concentration dependence of the coiling probability in 
detail. This is the problem of the excluded volume effect in concentrated 
solutions; it will be a subject of the next step of the present work. 
Finally, a theory similar to mine has been presented by Staverman,” 


which leads to 
by} 
m=!(1-149) (84) 
zZ n n 


corresponding to 53. However, the Staverman equation (84) also stands 
an overestimation of oo , as well as the present theory (53) does. 
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CONCENTRATION DEPENDENCE OF POLYMER 
DIMENSIONS 


Marshall Fixman 
Mellon Institute, Pittsburgh, Pa. 


Introduction 


Although several attempts have been made to generalize the theory of 
the excluded volume effect’ to include the dependence of the excluded 
volume factor a on polymer concentration,’* none has been free of such 
gross oversimplifications as, for example, a complete severing of the links 
of the polymer chain or the limitations of a perturbation calculation. 

A calculation of the radial distribution function on the basis of the real- 
istic, if approximate, model of Flory and Krigbaum’ was previously used 
to compute the osmotic pressure. The hope was expressed that the use 
of the osmotic pressure to compute a free energy and the minimization of 
the free energy with respect to a could yield the concentration dependence 
of a. The method has not been pursued to completion, at least not for 
high-polymer concentrations in “good” solvents, but such progress as has 
been made is here reported. Alternative calculations that avoid some of 


_ the present defects will be suggested. 


The Model 


The physical model is a slight generalization of that used previously in 


’ ; % c 7,8 : ; 
_ discussions of the osmotic pressure. A particular polymer molecule is 


represented as a continuous distribution of segments, such that the number 


~ of segments in a unit volume is a Gaussian function of the distance to the 


wa 


ee Te 


center of mass of the molecule. The total interaction energy of a pair of 
molecules is obtained by integrating the interaction energy density over all 
space. ‘The energy density is assumed to be proportional to the product 
of the two segment densities at the given point. The assumptions result 
in an intermolecular potential that is a Gaussian function of intermolecular 
separation. This essentially completes the specification of the model 
previously used, inasmuch as the parameters that enter into the potential 
were taken to be the same for all molecules and configurations of the sys- 
tem. Now, however, it is desired that the size of a polymer molecule be 
subject to fluctuations; the possibility is introduced that the mean dimen- 
sions can vary with polymer concentrations. A rather simple mechanism 
is at once available for the description of size fluctuations. In the Gaus- 
sian function that gives the segment density distribution of each molecule, 
there is a single parameter proportional to the square of the radius of gyra- 


- tion or the end-to-end distance. This parameter can then be taken equal 


to its mean or most probable value multiplied by an expansion factor a; , 


- for the 7th molecule. 


657 
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The probability distribution of the parameters of the system is readily 
given, in a purely formal way, as a Boltzmann distribution of the energy. 
The total interaction energy of the system is a sum of intermolecular en- 
ergies Vi; = V(Ri;, ai, a;) and intramolecular energies V;(a;). The 
parameter a; may be taken, for concreteness, as the ratio of the actual 
radius of gyration of molecule 7 to the most probable radius, and R;; is the 
distance from i to j. It follows from the hypothesis of Gaussian segment 
densities that'” 


Vi; = Kila? + a?) *? exp [—KsRij(ai + a;') |] (1) 
Vi = Ksai* (2) 


The parameters K, , K2 , and K; are independent of configuration, although 
they may depend on molecular weight, 7’, P, and the chemical composition 
of polymer and solvent. With 1 and 2 and an a priori probability factor 


Il exp [— Kia; |da; (3) 
da; = Ana; da; (4) 


arising from the entropy of expansion of the molecule, the probability den- 
sity of a given configuration is 


f(Ri ++: Ry, a1: +: aw) 


5 
= Q”" exp [—K: 2 ai — (2 Vis + a V;)/kT| e 
where Q is a normalizing factor. 

It is hardly necessary to state that 5, although it already involves a 
considerable oversimplification of the internal structure of a polymer mole- 
cule, provides no easy access to f(a). If the polymer concentration is 
low, a cluster expansion of a” would be useful and no doubt could be car- 
ried out through the first power of concentration. An approximate method 
of calculation will be sketched, which is feasible at all concentrations; how- | 
ever, it has not been carried through because so much more numerical work 
would be required than in the method actually used. The proposed method 
is to isolate one molecule from the collection of N, say 1, and to suppose 
that, in the calculation of the probability of a1, all the other (NV — 1) 
molecules may be assigned the same a, the mean or most probable value. 
Evidently this assumption would be most reasonable at high polymer con- 
centrations. A further approximation, reasonable in the same spirit, gives 
for the probability of ay 


f(a) mg, Q’ exp| — Kua 7s IG > Vi) +; Vateu) |/ wr (6) 
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where 


(= Vit) = 10) f V (Rx, a, a)9(R, a1, «) dR. (7) 


Here p is the number density of polymer molecules, and g(Riz, a1, a) is 

the radial distribution function of molecules of expansion factor a around 

the center molecule, which has a different expansion factor a,. It could 

be obtained from the Born-Green-Kirkwood equations by much the same 

method as we have previously used® to obtain g(Ri, a, a). A maxi- 

mization of ay f(a) with respect to a then gives an equation for a, in 

g( Riz, o, aw), and therefore f(a,) would be completely known. This 

method can be carried through when molecule 1 is a chemical species differ- 

ent from the other polymer molecules and does not require the integration 

over p that is required to compute a free energy (see under The Free Energy, 

below). Moreover, it should be stated that only the two-polymer one-sol- 

vent mixture is likely to yield unambiguous light-scattering determina- 

tions of dimensional changes ((a;’) versus ps). 

A somewhat different point of view is actually taken here. Only one a 

will enter into the calculation, and it is used to characterize the ensemble 

- rather than a particular molecule. Fluctuations of the dimensions of in- 

- dividual molecules will be suppressed during the calculation of the free 

energy of a fictitious system, namely, a system of polymer molecules all of 

which have the same Gaussian segment density. The scale factor a that 

enters into the density distribution appears also in the intermolecular po- 

__ tential (Equation 1), and thereby in the free energy. After the free energy 

of interaction per molecule is computed, there is added to it the intramo- 
lecular free energy, and the whole is minimized with respect to a. 


The Free Energy 


> 

g A few equations will be required to establish the free energy of the solu- 

tion in terms of the osmotic pressure of the solute, the latter being known 

_ from a previous calculation.” By assumption, a is to be adjusted to its 

~ most probable value at constant 7 and atmospheric pressure. Conse- 

- quently, the Gibbs free energy G rather than the Helmholtz free energy is 
~ to be minimized. 


Mow 


a G(T, P) = Nua + Neue (8) 


where y; and pw: are the chemical potentials of solvent and solute, respec- 
4 tively. Let P = Po + 7m, where P» is the fixed atmospheric pressure and 
“a the osmotic pressure. G(7', Po) is required for the minimization with 
respect to a. From the definition of z it follows that uw. = fie the poten- 


tial of pure solvent at 1 atmosphere. 
Now a change in 7 may be imagined to result from a change in the num- 
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ber density of polymer molecules. The Gibbs-Duhem equation, which 
states in general that 


Nidu + Nedus = —SdT + VdP (9) 
reduces under osmotic conditions to 
Nodpo cas Vdr (10) 
Therefore 
p 
a(T, Pop) = | pt de + wa(T,Po, po) (11) 
“Po 


where py is a concentration so small that the interaction energy of polymer 
molecules becomes negligible. Equation 8 becomes, with di; 


Pp 
G(T, P) = Nu ae Nope(T, Po, po) oe Na | p- dr 
Po (12) 
{2 = BP 0 +- tis! 
Suppose now that the osmotic conditions are removed (the semipermeable 


membrane is sealed) and the closed system is expanded at constant 7 to 
the pressure Po. 


dG = —SdT + VaP + >. nd; 
(13) 
= Var 
Therefore 
G(T ,P) — G(T, Po) = Vr (14) 


if the compressibility of the system be neglected. Equations 14 and 12 
give 
Gy; Po) =a Ny = Nouo( T, Po; po) + Nilp ‘x =e po (po) 
p (15) 
+ N, | p dp —Vr 
Po 


after a partial integration of | p dr. That part of G(7, Po) which de- 


pends on the intermolecular potential is, per molecule, 


ehces 
GC, = p dp (16) 
Po 


This is the result_one might have expected since, were a the hydrostatic 
pressure of a one-component system, G2 would be the Helmholtz free energy. 


The intramolecular part of the free energy per molecule is, according to 
Flory, 


veh) Un i at hee LL 
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G; = kT[Ka*Z — 3ina + 34(a’? — 1)] (17) 


where K is a numerical constant (to be taken as 134/105 when a determines 
the radius of gyration, in order to gain agreement with the perturbation 
theory’” of a), and 


Z = (3/24R,)*6n? (18) 


for a polymer chain of n segments, mean square end-to-end distance R,? 
in the random flight approximation, and segment-segment excluded vol- 
ume 6. Thus the quantity to be minimized with respect to a is AG = 
G. + G;. 


p 
AG = / p «dp + kT | Ke’ — 3lna =f : (a? = )| (19) 
Po 
The Osmotic Pressure* 


The osmotic pressure has been computed,® approximately, on the basis 
of a potential of form 1: 


V(R) = kT(4/x'”)A exp (— BR’) (20) 
A = 7.18a °Z 


(21) 
B = 9.61(a’Ry)* 


where the numerical constants in 21 have been chosen to agree with per- 
turbation calculations” of the osmotic second virial coefficient. An 
approximate radial distribution function was used 


g(R) = 1 — aexp (—éBR’) (22) 


and, by a variational procedure, the parameters a and 6 determined in 
terms of A, B, and p. Equations 20 and 22 gave an osmotic pressure 


a = pkT — (p'/6) J RV'(R)g(R) dR (23a) 
or 
ae = pkT + QnkTPAB 1 — a(1 + 8)? (23) 
Substitution of 23 into 19 gives 


(kT) AG a if cde as 2nAB? [ ee re eel See aR 
PO BG oe 


=~ | Ka‘ =] 3ine + : (a 1)| 


* In this section, 7 designates the osmotic pressure only on the left-hand sides of 
Equations 23 and 23a. Elsewhere z is the pure number. 
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The equation of d(AG)/da to zero now determines a. Note that AB *” 
is, according to 21, independent of a. 


p 
(nr) 1 2D) _ — apa? 3 | alt + 0)*" dp 
Qa 0 


— 3[KZa* +a — a] (25) 
= 0 


where po has been allowed to approach zero. 

The procedure that has been used here will naturally raise several ques- 
tions as to self-consistency. In part, this doubt results from the approxi- 
mations’ being lumped together and thereby obscured, relative to the pro- 
cedure suggested above under The Model. Some of the difficulties and, 
occasionally, a resolution of them can be pointed out. Equation 23a re- 
quires the assumption of pairwise additivity of the potential V(R) of 
average force between solute molecules. It also requires, almost as a 
matter of definition of the potentials, that V() not depend on concentra- 
tion, as here it eventually will, through a. It nevertheless turns out that 
the basic question of self-consistency is favorably resolved for a very simple 
reason. The reason is that a is not treated as a function of concentration, 
but rather as an arbitrary parameter, up to the point where AG is fully 
known, 24. This means that the a that occurs in 24 explicitly, and im- 
plicitly in A, B, etc., will contribute nothing to the first derivative of AG 
with respect to p(@AG/da = 0). Therefore the (osmotic) pressure ob- 
tained from AG by differentiation with respect to p will agree with the virial 
theorem pressure, 23a and 23. For a more detailed discussion of a similar 
problem, see Mayer and Careri.”® 

A solution of 25 is quite difficult to obtain. Only for the case of weak 
interactions (A < 0.2) is the solution reasonably straightforward and the 
result for every A accessible with one calculation. Only one other situa- 
tion will be investigated, namely, the case of fairly large A(A > 1.2), at 
not too large p. It turned out here® that a & 1 up to such large p as to 
yield considerable changes in @ or, for very large A, to exceed the a priori 
limits of the chosen potential V(R) (a polymer volume fraction of ca. 
0.1). 


Small A 


When A is small, a is small and becomes smaller as p increases. The 


equations previously given® for a and 6 may then be linearized with respect 
to a, and take the form 


atz + of2 — pahe 
ats + ofs — pah, 


(26) 


Fe yeas EASY 


és 
c 
“ 
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where ¢ = a ° and 


ty = =a x” exp (—2a”) dx (27) 
“0 
fo = AN(5/1 +8)"; fa = 34(8/1 + d)fr (28) 
he = 9 ABT (2 + 8) 7; hg = 36(1 + 8)(2 + 8) "hy (29) 


A' = ofA; B' = o'B (30) 
The definitions of A’ and B’ make them independent of variations in a. 
Although 25 apparently requires a and 6 in terms of a and p, it turns out 
to be simpler to obtain a and p in terms of 6 and a, and change the p in- 
tegration to a 6 integration. From 26 
Pre Ban «AU (lg id) (2.423) (31) 
a = —2fr[i6(1 + 5)]* (32) 


and therefore 
n =a(1 +5)" = (—2A4’/i2)g8"7(1 + 8) (33) 


Notice that the dependence of 7 on a is a simple proportionality to ¢ = 
a ’, as follows from 33 and 31. That is, 31 shows that 5 is independent of 
a; for a given polymer (A’, B’), 6 is a function only of p. 

Equation 33 gives in 25 


p 
a —a = KZ — 471, A”"B i 8'"(1 + 8) dp (34) 
0 
A transformation of variables p — 6, through 31, gives 
5(p) 
5 PP 3 me ee A” et eee 9 5 3/2 
a a O.1777A 7 ‘ (2 + 6) (35) 
(7a — 1) + 8) ds 


where A’ = 7.18Z (Equation 21) and K = 134/105 have been introduced. 
It proves convenient to tabulate the solution for a not against p but against 
the dimensionless concentration 


y = ded Cyaan (36) 

or 
y= (2/2) (16A') (6 — 1){2 + 8)-” (37) 
as follows from 31. The integration in 35 was performed numerically, 
and then a simultaneous tabulation of (a*° — a’) /A’ and »'A’ against 6, 


35 and 37 respectively, gave as a result a tabulation of (a — a’) /A’ against 
yA’. To an order of approximation consistent with the previous lineariza- 
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tion with respect to a, (a° — a°)/A’ & 2(a — 1)/A’, and FicuRE | 1s 80 
computed. 

Although the actual dependence of @ on p will be governed by the scale 
factors A’ and B’, one may conclude from FIGURE 1 that as p increases from 
zero, a at first decreases rapidly, and then more and more slowly. It is 
also clear from FIGURE 1 and Equations 35 and 37 that a approaches asymp- 
totically a value less than unity, as p> %. This was originally a rather 


0.! 
0.06 
0.02 

=0,02 


-006 


(a —1)/A! 


-0,10 
-0.14 


-0.18 


-O22 
Oo O8 1.6 24 32 40 48 56 
v'N 
Figure 1. The expansion factor a plotted against the dimensionless measure of 


concentration v’. The scale factor A’ is an intermolecular potential parameter; 
A’ < 0.2. 


surprising result, as we interpreted the behavior of g(R) at high concen- 
trations—namely, g(R) — 1, all R, as p > »—to mean that the effect of 
segment-segment interactions on a would become unimportant. That is, 
we expected the Flory result, a > 1 as p > », tobeconfirmed. However, 
some pondering over the details of the calculation will make it clear that 
only by a rather fortuitous cancellation of two tendencies could the result 
be a—1. The intramolecular forces expand the molecule, and the inter- 
molecular forces cause the molecule to contract. Furthermore it is not 
true that g(#) — 1 implies unambiguously a random distribution of seg- 
ments. If the center of one molecule is held at a fixed position, it 7s im- 
plied by g(R) — 1 that the segment density contributed by other mole- 


~~. ee 
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cules becomes uniform in space as p > 0 ; nevertheless, the fixed molecule 
still maintains its Gaussian density distribution. Although the fixed cen- 
tral molecule contributes only a small part of the total segment density at 
large p, that small part may and evidently does play a disproportionately 
large role in the determination of its own degree of expansion. 

A reminder is in order that the Flory-Krigbaum derivation’ of V(R) did 
not include, in the local energy density, entropy of mixing terms beyond 
the second power of volume fraction. Quite possibly these higher terms 
will make a — 1 as p > ~, in a more complete theory. 


Large A 


Another and different simplification of the equations that determine 
a.and 6 is possible if the solvent is very good. For large A (roughly A > 
1.2), the value of a was shown to be essentially unity for values of p as 
large as need be considered. However, it turns out that quite large changes 
ini a occur as the concentration is increased, and the intermolecular poten- 
tial energy, through its dependence on a, becomes concentration dependent 
to an extent that cannot be taken as a small perturbation. The process of 
compensating for these changes in a requires considerable although straight- 
forward numerical labor. Only a few results are available at present; a 
more extensive tabulation of a (and 6) against »’, for a sequence of values 
of A’, will in due course be submitted elsewhere. 

For a = 1, 6 satisfies® 


Leola apg m0 (38) 


where 


Il 


ie | In (1 —e¢™)a’e” dx = —0.1450 (39) 
0 


= 34A4'Bo*? P71 + 8)(2 + 8)? — (1 + 26)(1 + 8) 
OR ese) aes 


ay 
I 


(40) 


and f, is given in 28. 
As before, the p integration in 25 will be converted to a 6 integration. 
The derivative with respect to a or ¢ in 25 requires 


(2) = filo(ag/a8) — $(afx/98)\" (41) 


from 38. Also from 38, 


(2) = —glp(dg/d5) — o(dfs/d5)] (42) 
Op) ¢ 
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(#) = (-wa) (2), (43) 


The derivative with respect to a in 25 may be eliminated with 43, and then 
the p integration converted to a 6 integration with (05/dp)4dp = dé. The 
result is 


and therefore 


5(a,p) 
a? sa) 41347105) ee [ I(6) dé (44) 


6(a,0) 
where 


d BY siren 
(i + 8)(2 + 8)? — (1 + 25)(1 + 8) + 35) 


The limits on the 6 integral are determined as functions of p and a from 
38. In principle, the method of solution would be to choose A’ and B’, 
and the concentration p at which a is desired. Then 38 would provide the 
values 6(a, 0) and 6(a, p) if a were known. One could then guess at a 
until the proper solution of 44 were found, for the given A’, B’, and p. In 
practice the procedure has been to cast 38 into a form that gives 6 as a 
function of A = A’/a’ and »y = va’. A division of 44 by a’, and the sub- 
stitution Z = A’/7.18 gives 


I(6) = (45) 


6(A,r) 


a? — 1 = O1777A (1 — 1587 | 168) is) (46) 


5( 4,0) 
For a sequence of given A’s, there were tabulated (1) v versus 6, from 38 
and (2) a versus 6(A, v) from 46. The result is a tabulation of a and 6 versus 
vforagiven A. It is easily seen that to keep A fixed as 6(A, v) and v change 
means that A’ = a’A is changing. In principle this creates no problem; 
one need only tabulate A’ alongside (A, 5, v, a), and then search the tables, 
interpolating where necessary, for a preassigned value of A’. The result 
will be a series of entries (A’, 5, v, a), all having the same value of A’ but 
successively increasing values of 6 and vy and successively decreasing values 
of a. The computation of »’ = v/a’ is then readily made for each entry, 
and one finally has for an assigned A’, 6 and a tabulated against the good 
measure of concentration v’. Unfortunately, as I indicated, my initial 
tables were not extensive enough to carry v’ to adequately large values. 
However, even the fragmentary results here given are quite sufficient to 
permit qualitative inferences of the concentration dependence of a in 
‘“‘go0d”’ solvents. 

I have tabulated a against v for a sequence of A’s; some of the results 
are given in FIGURE 2. They probably are reliable only for small » or small 
decreases in a. That is, the abscissa could be related to p by taking B 
equal to its value at infinite dilution, if p is small. 


| 
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Figure 2. Alpha plotted against » = ay’ for various values of A = A’ /a3. 
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Fiaure 3. Solid curve: a plotted against » for A = 1.8 held constant. Dashed 
curve: a plotted against »’ for A’ = 1.8 X [a(v = 0)]* held constant. 
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For one value of A, A = 1.8, I have gone through the procedure of “re- 
normalizing” the A and » arguments to A’ and 1’, to as large a »’ as was 
possible. Every point on the dashed curve in FIGURE 3 refers to a single 
polymer-solvent system (fixed A’ = 1.8 X [a(0)]*) and the ordinate v’ 
is proportional to p with a fixed proportionality constant. It is seen that 
only very small changes in the curve are induced by this “renormalization” 
other than a simple reinterpretation of the abscissa. Evidently, for qual- 
itatively unknown reasons, the two changes in A and » tend to cancel each 
other. There is, however, one very significant effect at the higher concen- 
trations. ‘The renormalized curve begins to flatten out, and it is to be ex- 
pected that a will asymptotically approach a limiting value at high con- 
centrations. The qualitative behavior of a is, accordingly, the same in 
poor (small A’) or good (large A’) solvents: an initially rapid decrease to 
values below unity and a subsequent leveling off at high concentrations. 

There is one very useful although tentative conclusion that may be 
drawn from the dashed curve of FIGURE 3. Up to concentrations in which 
appreciable overlapping of polymer domains is forced to occur (v’ ~ 1) 
and a is significantly lowered, the dashed curve is quite close to a straight 
line. This means that a virial expansion of a in powers of p could also 
give information about experimentally detectable decreases in a. 


Summary 


The mean dimensions of chain polymer molecules in dilute solution will 
usually be expanded beyond their random flight values because of the 
mutual repulsive forces exerted by chain segments. These repulsive forces 
may, however, be screened by the presence of other polymer molecules in 
solution, and accordingly the mean dimensions depend on polymer concen- 
tration. I have discussed the evaluation of this concentration dependence 
by a minimization of the total free energy of the solution. To the free 
energy of expansion of an isolated polymer molecule is added the free energy 
of interaction with other polymer molecules. The latter part of the free 
energy is known from a theoretical expression for the osmotic pressure and 
its dependence on polymer dimensions. The well-known dimensional ex- 
pansion factor a, which is greater than unity for a polymer molecule in 
dilute solution in a good solvent, will dip below unity when the concentra- 
tion is increased. At a sufficiently high concentration beyond the reach 
of the present theory, a may again approach unity. 
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THEORY OF THE NON-NEWTONIAN VISCOSITY OF 
HIGH POLYMER SOLUTIONS 


Bruno H. Zimm* 
School of Science and Engineering, University of California, La Jolla, Calif. 


Peterlin in an earlier paper in this monograph has already discussed the 
possible theoretical sources for the non-Newtonian viscosity of high poly- 
mer solutions. I, too, have been investigating methods for calculating 
numerically the non-Newtonian effect of the hydrodynamic interaction. 
Since the calculations are still in progress, this note attempts merely to 
summarize some of the more striking preliminary results. 

The exact form of the hydrodynamic interaction tensor, which gives the 
velocity of flow of the solvent produced at point j by the application of a 
unit force at point 2, is 


(i/|r)G.— rn/9), 


where r is the vector from i to 7. This form, originally used by Burgers 
for the case of simple single particles, raises serious algebraic obstacles in 
the way of a solution of the equations of motion of any system of more than 
two or three centers. Kirkwood and Riseman? introduced a simple ap- 
proximation to avoid the difficulty. They replaced the tensor with its 
average value in the molecule at rest, which depends only on the average 
distances in the molecule, a constant for any pair of chain segments, 
With this approximation the viscosity of high polymer solutions turns out 
to be Newtonian.® 

I have tried two different ways of improving this approximation. Nei- 
ther leads to exact results, but one seems much more promising than the 
other. The first and less promising method is also the more obvious; it is 
to use the Kirkwood-Riseman approximation, but with the average dis- 
tances pertaining to the molecule in its steady state at a given rate of shear 
rather than at rest. In shear the molecule is anisotropically expanded, so 
that different averages are needed for the different components of the 
tensor. The result of the rather involved calculation is puzzling. At low- 
to-moderate shear rates the viscosity drops, as in the experiments, as a 
result of a complicated interplay of different terms; as the shear rate rises 
further, however, a monotonic rise in the viscosity occurs. This rise is 
due to a general weakening of the interaction as the molecule expands and 
as its interior becomes more and more exposed to the flow of solvent. In so 
far as I know, this rise has no counterpart in experiment. 

* This work was done in part while the author was associated with the General 


Electric Research Laboratory, and in part while he held a Visiting Professorship at 
Yale University, New Haven, Conn. 
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It therefore seems in order to search for another method of improving 
the Kirkwood-Riseman approximation, and the following approach has led 
to encouraging progress. We seek to represent the interaction tensor as a 
polynomial in r and its components. In principle the representation could 
be made as exact as desired by using a sufficient number of terms in the 
polynomial, but the computational difficulties increase rapidly, so that we 
are forced to limit the number drastically. From this point of view, the 
Kirkwood-Riseman approximation consists of using a one-term (constant) 
polynomial; I have now extended the polynomial to two terms, including 
one in 7’, as well as a constant. This form of interaction can be expressed 
in terms of the normal coordinates of the problem, and the equations of 
motion can be solved. 

The result is gratifyingly simple, and it has at least the qualitative fea- 
tures of the experimental data. The viscosity if found to drop from its 
initial value as the shear rate is increased, at first falling rapidly, then more 
and more slowly, but never rising again as in the first method of approxi- 
mation. The drop in viscosity first becomes appreciable, about 5 percent, 
when the shear rate reaches the reciprocal of the longest relaxation time of 
the chain;? something like this might have been anticipated on physical 
grounds. 

This apparent agreement with experiment is encouraging, and it indi- 
cates that further refinements of this method should be examined. The 
obvious next step is to include anisotropic terms in the polynomial, and 
then higher powers of r._ While the amount of labor involved in carrying 
out these steps may be considerable, there do not appear to be any funda- 
mental difficulties. 
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RECENT RESULTS IN THE CONTINUUM THEORY OF 
VISCOELASTIC FLUIDS* 


Bernard D. Coleman 
Mellon Institute, Pittsburgh, Pa. 


Walter Noll 
Mathematics Department, Carnegie Institute of Technology, Pittsburgh, Pa. 


INTRODUCTION 


This article is concerned with phenomenological aspects of the mechan- 
ical behavior of viscoelastict fluids. The discussion is limited to those 
special circumstances in which all nonmechanical influences can be neg- 
lected; that is, no consideration is given to thermal, chemical, and elec- 
tromagnetic phenomena. 

Our procedure is as follows. After reviewing in Section 1 some kine- 
matical prerequisites, we give in Section 2 a mathematical definition of the 
concept of a simple fluid. This definition is sufficiently general to include 
perfect fluids, Newtonian fluids, and viscoelastic fluids as special cases. 
In Section 3 we discuss the behavior of a simple fluid in steady simple 
shearing flow. It turns out that for incompressible fluids several steady 
flow problems can be solved by direct appeal to the definition. Some 
properties of these solutions that may be of interest to experimenters are 
discussed in Section 4. Recent results on the behavior of simple fluids in 
periodic motions are given in Section 5. In Section 6 we formulate a 
general “smoothness assumption,” which makes it possible to prove a 
theorem that seems to justify the intuitive notion that for most fluids the 
theory of the Newtonian fluid should be a first-order correction to the 
theory of perfect fluids in the limit of “slow motions.” A rigorous pro- 
cedure for determining the form of higher-order corrections is sketched in 
Section 7. 

It is not our intention to supply here formal mathematical proofs of all 
our assertions, for this article is not intended primarily for mathematicians 
working in continuum mechanics. Our goal is to summarize in an article 
of reasonable length certain recent results in the mechanics of fluids that 
may be of interest to polymer physicists. The omitted proofs not already 


published elsewhere will appear shortly in the Archive for Rational Me- 
chanics and Analysis. 


* The work reported in this paper was supported in part by grants from the Air 
Force Office of Scientific Research, Washington, D.C., under Contract AF 49(638)541 
with the Mellon Institute, Pittsburgh, Pa., and from the National Science Founda- 
tion, Washington, D.C., under Grant NSF-G 5250 to Carnegie Institute of Technol- 
ogy, Pittsburgh, Pa. 

_ | The term viscoelastic does not have a precise meaning except in linearized theo- 
ries. Here we use the word to indicate that we are dealing with materials that do 
not obey the classical laws of fluid mechanics or elasticity. 
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Although we omit some proofs, we endeavor not only to state each result 
precisely but also to make explicit the mathematical hypotheses from which 
it was deduced. 

Before proceeding to our main topics and becoming involved in mathe- 
matical symbolism, we first try to give in this introduction a rough idea as 
to what we mean by a simple fluid. 

The word fluid, as used in everyday speech, does not have a precise 
meaning. It might be said that an essential property of fluids is absence 
of preferred configurations; thus a fluid may be defined as a substance with 
the property that every configuration is an undistorted configuration. 
Physicists sometimes say that a fluid is a substance that cannot support a 
shearing traction without giving way to it. Some dictionaries define fluid 
asa substance capable of flowing. In order to give these definitions exact 
meaning and to see whether they are related, a general mathematical theory 
of the mechanical behavior of continuous media must be presumed. Such 
a general theory now exists. 

Since the days of Euler and Cauchy,* almost all phenomenological 
theories of the mechanical behavior of continuous media have rested on 
constitutive assumptions that may be seen to fall as special cases within a 
general hypothesis framed only much later and called the principle of 
determinism for the stress. The stress S(¢) at a material point X at time ¢ 
is determined by the past history of the motion in an arbitrarily small 
neighborhood of X. 

For many materials} the principle of determinism is not needed in such 


ee a, general form; in most theories that have been proposed, S(t) is deter- 


_ mined at X by the history, up to time ¢, of only the first spatial gradient F 


at X of the displacement function. JF is called, for short, the deformation 


gradient. 


} 
= 


The deformation gradient (7) for X at time 7 is a second-order tensor; 

it is defined as the gradient of the position vector = A ff, 2), at time 7 

of X with respect to the position vector, q = fonda pot X in whatever 

configuration we have taken as reference. If we are using a Cartesian 
coordinate system, the components F"; of F(7) are given by 

ne we 

F 7 = ag (1) 

Of course, F(r) depends not only on the configuration at time 7, but also 

on the configuration taken as a reference. By the history of the displace- 

ment gradient up to time ¢ we mean a specification of F'(7) for all 7, —% < 


oY llent review of the work done up to 1952 see Truesdell.1? A com- 
teeteacive. afticls by W. Noll and C. Truesdell entitled ‘‘The Non-linear Field Theo- 


ries of Mechanics’’ will appear soon in Fliigge’s Encyclopedia of Physics. 


+ Among the materials excepted are dilute gases with large velocity gradients* 4 
and fluids with very large density gradients® (for example, fluids in situations in 
which capillarity effects cannot be neglected). 
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+ <1. In rough terms we can say that the information contained in the 
history of the displacement gradient at a material point consists in a specifi- 
cation of all the rotations and stretches (both their magnitudes and direc- 
tions) experienced by the material point in the past and, furthermore, the 
specification of the times at which these rotations and stretches occurred. 
Substances for which the history of the deformation gradient determines 
the stress have been called simple materials,’ and the general phenomenolog- 
ical theory of their mechanical behavior has been developed considerably 
in recent years.” 

The functional* &, which establishes the correspondence between the 
history of the deformation gradient and the stress tensor will in general 
depend on which configuration M is taken as a reference in computing the 
deformation gradient. For certain materials such as diamond and cross- 
linked rubber{ there are certain special configurations which, when taken 
as reference configurations, will give the functional & » a particularly simple 
form. These special configurations have the property that if a material 
point X has been held in one of these configurations for all times 7, — 2° < 
r < t, then S(t) at X must bea hydrostatic pressure. We call such con- 
figurations undistorted. 

For those materials for which only certain special configurations are un- 
distorted, the functionals  y have the following property: Suppose we are 
given the history up to time ¢ of the deformation gradient F relative to some 
general configuration M’ which may not be undistorted. Then we cannot 
use ® y to compute S(¢) unless we also know the displacement gradient of 
M’ computed with the configuration M taken as a reference. This is not 
the case for all materials, however. For some substances a knowledge of 
the history of the deformation gradient /’, computed relative to an arbitrary 
configuration M’, completely determines the stress, provided that we know 
merely the mass density corresponding to M’ or, equivalently, provided 
that we know merely the determinant of the deformation gradient relating 
M’ to M.t For such a material there exists a single functional &, which 
gives S(t) as a function of the density p of a configuration M’ and the 
history of the deformation gradient F relative to M’: 


S(t) = & (F(r); p). (2) 


Here F is a tensor-valued function of 7, —%« < 7 < #t, and p is a number. 


* A functional here is simply a function whose arguments are functions and whose 
values are tensors; the stress tensor at time ¢ is taken as a functional of the history of 
the deformation gradient up to time ¢, which may in turn be regarded as a tensor-val- 
ued function of a time variable r, —» < 7 < t. 

+ We assume here'that the temperature and our time scale are such that a negligible 
fraction of chemical bonds can break during any mechanical experiment. 

t It turns out that such a substance must be isotropic; the condition on Ry stated 
here is, however, much stronger than isotropy. In the language of Noll,® our eondi- 
tion is equivalent to the assertion that the isotropy group of & is the full unimodular 
group. 
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If Equation 2 holds at the material point X, then we say that X is a ma- 
terial point of a simple fluid.* 

Since the functional ® in 2 does not change if the configuration M’ is 
changed, any configuration may be used as a reference. In certain prob- 
lems it may be convenient to take the present state as a reference, then p 
is just the density at X at time t, and F isa function which has the value J 
(the identity tensor) at time t. 

If a simple fluid has remained in a single configuration M for all times r, 
—° <7 X< #, then the stress S(t) is a function of only the density cor- 
responding to M. Furthermore, it can be shown by appeal to the principle 
of material objectivityt that S(t) must then reduce to a hydrostatic pres- 
sure. Since this result holds for any configuration M, it follows that every 
configuration of a simple fluid is undistorted. A simple fluid that has 
been at rest forever obeys the laws of classic hydrostatics. It follows from 
these remarks that if the stress in a simple fluid is not a hydrostatic pres- 
sure, then the fluid cannot have always been at rest in one configuration. 
Thus a simple fluid cannot support a shearing stress indefinitely without 
flowing. 

The concept of a simple fluid is a very general one. Simple fluids can 
exhibit those phenomena that are now attracting the main attention of 
experimental rheologists: stress relaxation, non-Newtonian viscosity, and 
normal stress effects. In particular, we believe that the theory of simple 
fluids can be successfully applied to rationalize the treatment of mechanical 
data on pure amorphous, not cross-linked, polymers and polymer solu- 
tions. f 

We do not wish to give the impression that the theory of the simple fluid 
is so general that it is the only theory that can give precise meaning to the 
intuitive notion of fluidity. Consider the dictionary definition of a fluid 
as “a substance which is capable of flowing.”’ There is, of course, no unique 
way of making this definition precise. We can suppose that what is really 
meant by the condition of being “capable of flowing” is the condition of 
sustaining a steady simple shearing flow in at least one direction when 
subjected to appropriate shearing tractions of bounded magnitude.§ A 

* The materials that we call here simple fluids were defined by Noll,* and were then 
called fluids. That article contains a rigorous presentation of the ideas roughly 
sketched in this introduction. What we now call simple fluids have also been called 
general fluids!!! because they not only include perfect fluids, Newtonian fluids, the 
Reiner-Rivlin fluids,” and the Rivlin-Ericksen fluids,” as special cases, but 
because they are also much more general than all these special fluids in that they are 
capable of representing long-range hereditary effects such as stress relaxation. _ 

} This principle,® which is also called the principle of material indifference, is a 
mathematical statement of the idea that the response of materials is independent of 


e ob that is, invariant under changes of frame. : 
ee fies possible that the theory of simple fluids cannot be successfully applied to the 


e physics of liquid crystals. My colleague and I, however, are not sufficiently familiar 


with the mechanical behavior of such substances to hazard a categorical statement 


ag $ Note that we do not require that the shearing tractions be constant in time. 
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substance that meets this condition might be called a hypofluid. Every 
simple fluid is a hypofluid but not vice versa. One can even imagine hypo- 
fluids that may remain at rest indefinitely without obeying the laws of 
classical hydrostaties.* 

The materials defined by Ericksen in his recent articles “"" on anisotropic 
fluids are special hypofluids. Ericksen calls these substances anisotropic 
fluids because they have the property that the stress at a point X depends 
on a vector n, which can be interpreted as determining a preferred direc- 
tion at X. In each case, the vector n depends in a specified way upon the 
history of the motion of X. In some of these fluids, n reduces to the zero 
vector when the fluid has been at rest forever.t Depending on the choice 
of values for the material constants occurring in Ericksen’s constitutive 
equations, the hypofluids considered by him may or may not be simple 
fluids. Ericksen’s work illustrates the kind of behavior that hypofluids 
may exhibit when they do not reduce to simple fluids. 

Notation. We denote vectors and points in Euclidean space by boldface 
Latin or Greek miniscules: v, x, —, --- . Second-order tensors are denoted 
by lightface Latin majuscules: A, B, Q, I, F, C.(r), ---. For the trace 
of a tensor F we write trF and for the determinant of F we write detF. 
The transpose of F is denoted by F”. We say that a tensor S is symmetric 
if S = S". The identity or unit tensor is written 7. The symbol Q is 
used only for orthogonal tensors; that is, tensors such that QQ’ = Q’Q =I. 
The physical components of the tensor 7 are denoted by 7;;. We use 
both || F:(7) || and || Fc»';(7) || to denote the matrix of components of the 
tensor F;(7) relative to a preassigned basis. The gradient operator and 
divergence operator are denoted by V and div, respectively; we sometimes 
add a subscript to V to emphasize the independent variable with respect to 
which we differentiate: V, . 


16,17 


1. Somes ELEMENTARY KINEMATICS 


Consider the material point X which occupies the position x (in Euclidean 
space &) at time ¢t. Let us suppose that at time 7, X occupied{ the posi- 
tion £in &. For the dependence of & on x, t, and 7 we write 


E = x(x, 7). (1.1) 


The function x, may be called the displacement function relative to the 
configuration at time t. 


If the velocity field v is given as a function of the position & in & and the 
time 7, 


v = v(&, 7), (1.2) 


* Furthermore, a hypofluid need not be isotropic; that is, for a hypofluid there need 
not exist a reference configuration M such that the isotropy group® of ®y4 contains 
the full orthogonal group. : 

} Ericksen’s use of the word anisotropic is not the direct antithesis of the use of 
the word isotropic by Noll.§ 


t The letter ¢ will be regarded as the present time and 7 a time earlier than ¢. 
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then the position &(7) of each particle at time T; 
E(r) = x(x, 7), (1:3) 


can be found if its position x at time tis known. For E(7) is simply that 
solution of differential equation 


de(r) _ 


ae vlE(r), 7] (1.4) 
which satisfies the end condition 
E(r) ear = x1(X, t) = X. (1.5) 


Thus a knowledge of the position x of each material point at time ¢ and 
the velocity field (1.2) for all 7 is equivalent to a knowledge of the dis- 
placement function defined in 1.1.* 

The gradient F.(7) of x:(x, 7) with respect to x, 


F(r) = Vxxe(X, ae (1.6) 


is the deformation gradient at time 7 relative to the configuration at time t. 


This second-order tensor F',(7) can also be called the deformation gradient 
at time 7 computed with the present state taken as reference. The follow- 
ing equation follows from 1.5: 


Fit). =f, Clee) 
Here J is the unit or “metric” tensor. 


Let p(7) and p(t) ke the mass densities at the material point X at times 
7 and t, respectively. It follows from a well-known result} in kinematics 


_ that these densities are connected by the formula 


ee IT 


Renee = Pat (1.8) 


We now define a tensor C,(r), which is called the right Cauchy-Green 
tensor at time 7 relative to the configuration at time t: 


Cir) = Fy" (1) Fe(r). (1.9) 


It follows from its definition that the right Cauchy-Green tensor is sym- 


metric. Furthermore, since we always assume that F’,(r) is invertible,t 


it also follows that C,(7) is positive definite. 


It follows from 1.7 and 1.9 that 
Cer = 1, (1.10) 


* We make these remarks because flow problems are usually defined in terms of 
properties of the velocity field, whereas in the theory to be presented here the stress 


_ is regarded as being determined by the history of the spatial gradient of the displace- 


ment function. ‘ 
+ See Truesdell,! p. 140, Equation 13.5. 
t That is, det(F:) + 0. 


678 Annals New York Academy of Sciences 


We shall often find it convenient to regard C,(r) as a function of the 
time lapse r between 7 and ¢: 


r=t—r. (1.11) 


We shall then assume that r > 0 (that is, that 7 < ¢). 
The derivatives A, of C,(r) at r = ft, 


= (- = Cit — r) 


=n (1.12) 


r=0 


os Clr) 


are the Rivlin-Ericksen tensors.* The tensor A, is familiar in classical 
hydrodynamics. It is simply twice the symmetric part D of the velocity 
gradient tensor V,v(x, t). To see this we write 

ied 
eer ( 


= Ly"(t)Fi(t) + T(t) F."(t) = Ly"(t) + L(t) (1.14) 


where L;(¢) is defined to be the time derivative of the displacement gradi- 
ent F,(7) at time ¢: 


Ar = £ Cur) PE (r)F) (1.18) 


5 


L(t) = © F(z) (1.15) 


tt 


By 1.6 we have 


ate © vee, es (1.16) 


tt 


Assuming x; is smooth in x and 7, we can interchange the order of differ- 
entiation in 1.16, and we get 


==. Wav (X shoe (1.17) 


Tt 


I(t) = Ve SE! (x,7) 


Thus 1.14 becomes 


A, = (Vxv(x, t))” + Vyv(x, t) = 2D. (1.18) 


To compute A,, for arbitrary ” one can differentiate F;7(r)F,(7) n times 
with respect to + and then put + = ¢. In this way we get 


A, = 2 (7) Ly "Lye (1.19) 


_ *See Noll,* Sections 8and9. Rivlin and Ericksen" were the first to recognize the 
importance of the tensors A, for n > 1 and used them extensively in their theory of 


isotropic materials. ; 
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where 
L,,(t) a g oe a— i 2; Og 
z aah (1.20) 
y(t) — Le 
For example 
Ape Ig tale Oly GIoae 
From 1.6 and 1.20 we have 
& 
Le = dr? Vx e( X, T) a (1.22) 
We assume the velocity field to be smooth in x and r._ Then, 
2 
Dee ye (xs) ee ve (1.23) 
dr T=t aot 


and, if we remember 1.17, it follows that A, is determined by the spatial 
gradients L, and Le of the velocity field v(x, t) and the acceleration field 
V(x, ¢) at x at time ¢. 

The function C; plays a central role in the theory of simple fluids, and 
the two tensors A; and A» occur frequently in special applications. 


2. SIMPLE FLUIDS 


Compressible Fluids 


For the purposes of the present article we can define a simple fluid as 
follows. 


Definition. The material at X isa simple flucd if the stress tensor S(¢) at 
X at time ¢ depends on the history of the motion through an equation of 
the form 


S(t) = $ [Cult — 7); ol] (2.1) 


. Here, is a functional which has for its argument the history of the right 


Sith lh Debt le ele ENE 
ere ‘ 


Cauchy-Green tensor relative to the configuration at time ¢ and the density 
at time ¢t. § is assumed to obey the following identity for all histories 
C(t — r) and for all constant orthogonal tensors Q: . 


QS [Cult — 7); pHIO” = GOCE HO"). (BB) 


In general, § will depend on the material point X, but we do not make this 


dependence explicit since our present, discussion is limited to homogeneous 
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bodies, that is, bodies that have the same constitutive equation at each 
material point. : 
Although this definition looks different in form from that given by Noll 
and discussed in our introduction, it has been shown* to be equivalent to 
the original definition. ; : :' 
Let us temporarily restrict our attention to situations in which the simple 
fluid has been at rest for all times 7, —2© <7 <¢. Then, 


C(t—r) =], O0<r< _o, (2.3) 


and the functional  [C.(t — 7); p(t)] reduces to a function of p(t): 
r=0 


S() = § [Ct = 7)5 pO = G Us a0] al 
= H[p(t)). 
Furthermore, since for all orthogonal Q 
vf oad QQ’, (2.5) 
Equation 2.2 becomes 
Q Us e(H1Q" = & Ws ald); (2.6) 
that is, 
QS(t)Q" = S(t). (2.7) 


Now, when a tensor S(t) obeys an identity of the form 2.7 for all orthog- 
onal Q, then S(t) must be a scalar multiple of the identity. Hence we 
have proved the following: If a simple fluid has been at rest for all times 
in the past, then the stress in the fluid is given by a hydrostatic pressure p, 
which depends only on the density, that is, 


S = —p(p)l. (2.8) 


Incompressible Fluids 


If the homogeneous body under consideration consists of an incompres- 
sible fluid, then we need not make the dependence of S(t) on the density 
p(t) = po = const. explicit. Furthermore the stress in an incompressible 
fluid is determined by the history of the motion only up to a scalar pressure 
p. Thus in the incompressible case we replace 2.1 by 


T(t) = © [Cut — r)} (2.9) 


* See Noll® (Sections 21 and 22). Equation 2.2 follows from the isotropy of simple 
fluids and the principle of material objectivity. A functional § that obeys Equation 
2.2 for all orthogonal Q is said to be isotropic. 


4 
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where 
T=S+o! (2.10) 


is called the extra stress. The tensor-valued functional © is now deter- 
mined only up to an arbitrary scalar-valued functional of Cit —r). We 
remove this indeterminacy by the normalization 


trl (4) = tr © Che = 1) |= 0, (2.11) 


It follows from 2.10 and 2.11 that p reduces to the mean pressure: 
p= —lWtrS, (2.12) 


and T is the deviatoric part of S. The functional § is of course still iso- 
tropic in the incompressible case; that is, 


Q © [C(t — r)]Q* = 8 [QCu(t — r)Q"), (2.13) 


for all histories C;(¢ — r) and all constant orthogonal tensors Q. 
By the same argument used in the compressible case, it follows from 2.13 
that when the motion reduces to a state of rest, that is, when C,(¢ — r) = T, 


T(t) = ) fare eat ee © (I) (2.14) 


and thus reduces to a scalar b times the identity, 
TG) = ot: (2.15) 
In the present case, however, it follows from 2.11 that tr” = 3b = 0. 
Hence, 
© = 0. (2.16) 
Since the motion of an incompressible fluid is always isochoric (volume 
preserving) it follows from 1.8, 1.7, and 1.9 that 
detC.(7) =4; —2o <7<i. (2.17) 


Hence for an incompressible fluid the functional § need only be defined 
over a domain of functions whose values are unimodular* symmetric posi- 


tive definite tensors. 
We remark that a useful, necessary and sufficient condition for the 


motion of a fluid body to be isochoric is that 


AT eA eee ee 


divv = 0 (2.18) 
- throughout the body at all times ¢. 


* That is, with unit determinant. 
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When we consider special flow problems in incompressible fluids we shall 
always assume that the body forces g (per unit mass) havea single-valued 
potential y: 


g = —grad y. (2.19) 
We shall also make use of the modified pressure ¢: 
d= pt py. (2.20) 


It coincides with the ordinary pressure p if there are no body forces. The 
dynamical equations will then take the form 


div T — grad ¢ = pv. (2.21) - 


3. SIMPLE SHEARING FLOW 


We now illustrate some of the concepts introduced in the previous 
sections by considering a simple shearing flow of an incompressible fluid. 
We define this flow by postulating that in some Cartesian coordinate sys- 
tem 2’, x’, x* the velocity field v(x) = {v', v”, v'} has the form 


v = 0, 
v = r(2'), (3.1) 
v = 0. 


Since v is assumed to be independent of ¢, simple shearing flow, as we 
have defined it, is a steady flow. Furthermore, since v clearly satisfies 
Equation 2.18, simple shearing flow is isochoric and thus compatible with 
incompressibility. 

Let é'(7) be the coordinates of the position at time 7 of the material 
point whose position at time ¢ has coordinates x’. Then 


f(t) = 2°. (3.2) 
On substituting 3.1 into 1.4 we get 
ae 
dr 
dg 
dr 
dg 
dr 
It is easily seen that the solution of the system 3.3 subject to the end con- 
dition 3.2 is 


= 0, 


= y(t), (3.3) 


= 0. 


t(r) oa 
B(r) = 2" + (r — t)v(2'), (3.4) 
B(r) = a. 
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It follows from the definition 1.6 that, in our present Cartesian coordinate 
system, the components F(»';(7) of the deformation gradient tensor F 1(T) 
are given by 

fe ok 
| oat (3.5) 
Equations 3.4 and 3.5 yield the following expression for the matrix of the 
components Fyy';(7) of Fi(r): 


1 WoW) 
; ] 
[| Per) |] = I Peo's(r) || = |} (re — the 1 OF, (3.6) 
0 Oval: / 
where 
dv(a') 
paee (3.7) 
Equations 1.9 and 3.6 yield 
litre —re 0 
Leroi are 2.1 +0 (3.8) 


0 0. a 


for the matrix of the components of C;:(t — r). Hence, in simple shearing 
flow, C.(t — r) is quadratic in r and has the form 

Cit —r) =I — 7A, + Wr'Ar, (3.9) 
where J is the unit tensor and A; and A, are the first two Rivlin-Ericksen 
tensors defined in 1.12. By 3.8 the matrices of the components of A, and 
A, are 


0 20 1ar02..0 
fArpP Seri 0 0], | Aol] = 2°10 0 Off. (8.10) 
OEs0n'0 Omtler0 


Whenever C(t — r) has the form 3.9,* Ci(¢ — 1) is completely deter- 
mined by A; and A,. When this is the case in an incompressible simple 
fluid, the functional § of 2.9 must reduce to a function } of the two tensors 
A,and A; : 


EG) =D [Ckb— 1) = Wri Aa)s (3.11) 
Furthermore, 7'(¢) must be independent of t. It follows from 2.13 that 


* There is a latge class of steady flows, including Poiseuille flow, concentric pipe 
flow, and Couette flow, for which 3.9 holds. ; 
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the tensor function } must be isotropic; that is, 
§(QA1Q”, QA.Q") = QH(A1, A2)Q" (3.12) 
must hold identically in A, , Az, and Q, where Q is orthogonal. From 2.11 
and 2.16 we have 
tr[h(A1, A2)] = 0, (3.18) 
h(0, 0) = 0. (3.14) 


It follows from 3.11 and 3.12 that whenever C,(¢ — r) is quadratic in r 
and the matrices of the physical components of A; and A: have the special 
form shown in 3.10, the matrix of the physical components of T(t) must 
have the form* 


Il 


Vi oT pa 
| T(t) || =|} Ta Te Of], Te = Tx (3.15) 
0 0 Ts 


where the 7';; are functions of x alone. By 3.13 we have 
Ty + T 22 ar T 3s = 0. (3.16) 


It is clear from 3.15 and 3.16 that the components 7;; of the extra stress 
may be expressed in terms of three independent functions of x: 


Tie 78), Ty, — T33 = o1(k), Tx — T33 = o2(x). (8.17) 


From 2.10 we see that for the components S,; of the stress tensor S(t) we 
also have 


Sie = Soi = 7(k) 
Su — S33 SS o1(k) Soo — S33 = o2(k). 


Because of the isotropy relation 3.12, the functions 7, 1, and o2 do not 
depend on which Cartesian coordinate system is used to obtain the simple 
form 3.1 for the velocity field. In other words, these three functions 
depend on only the material and not on the direction in which the material 
is being sheared. For this reason, 7, 0; , and o» are called material func- 
tions. 


The functions 7, o;, and o are not completely arbitrary. It follows 
from 3.12 thatt 


(3.18) 


t(—xk) = —7(k), oi(x) = o(—k), (3.19) 


* See Coleman and Noll!® p. 294; the proof given therein is the same as the proof 
of 5.20 in Section 5. That 3.15 follows from 8.11, the isotropy of }, and 3.10 rs also 
be proved using the Rivlin-Ericksen representation theorems for isotropic functions." 

t Coleman and Noll! pp. 294, 295; the proof is the same as the proof of 5.26 in 


Section 5. Equation 3.19 can also be proved using the Rivlin-Ericksen representa- 
tion theorems. 
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that is, that 7 must be an odd function, while o; and a are even functions. 
No further restrictions are imposed on 7, 01, and oz by the isotropy condi- 
tion 2.13. 

Equation 3.14 implies that the material functions r, o;, and o2 must 
vanish for x = 0: 


7(0) = o1(0) = o2(0) = 0. (3.20) 


A further restrictive condition on 7 follows from the assumption that 
the dissipation of energy must be positive. This condition is 


Kr(x) > 0, if x ~ 0; (3.21) 


that is, 7(x) must have the same sign as x. If it is assumed that 7 is twice 
continuously differentiable, then 3.20 and 3.21 imply that r(x) must be a 
strictly increasing function of x in some interval —K < «x < +x9 around 
x = 0. In this interval 7 will have a strictly increasing odd inverse 7’. 
Whenever 7 occurs in the remainder of this article we assume that we are 
in the range in which 7 has a single-valued inverse. 

The “‘shear-dependent viscosity,” used in the rheological literature can 
be identified with r(x) /k: 


n(x) = ee (3.22) 


It follows from 3.19 that 7 must be an even function. Since 7(0) = 0, 
if we assume that 7 is twice differentiable at x = 0, then 7 is differentiable 
at x = 0 and 


7/(O) = 0. (3.23) 
If we assume that co, and o» are differentiable then 
o1'(0) — o2'(0) = Q), (3.24) 


If approximations based on Taylor expansions of 7, 01 , and oz in a neigh- 
borhood of x = 0 are used to fit experimental data, only even powers of x 
can occur. 

Fr Let us examine the special forms taken by the material functions 7, 1 , 
and a» for some of those particular materials whose behavior is covered by 
the theory of incompressible simple fluids. 

For incompressible perfect fluids 


r(x) = oi(k) = oo(x) = 0 (3.25) 
4 for all x. 
For incompressible Newtonian fluids 
; = = const. > 0.- 
T(k) = km, 7 =C Ase 


oi(k) = o2(k) =0, 
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for all x; that is, the viscosity function 3.22 is a constant and the normal 
stress differences shown in 3.18 are zero. ees . 

In the incompressible case the theory of the Reiner-Rivlin fluids 
places no restrictions on 7 other than those which hold for all simple fluids. 
When we consider the normal stresses in simple shearing flow, however, it 
turns out that the Reiner-Rivlin theory yields the result 


Su = Soo . (3.27) 


12,18 


Thus by 3.18, Reiner-Rivlin fluids have the special property 
oi(k) = o2(k). (3.28) 


The theory of incompressible fluids of the differential type’ which was 
first considered by Rivlin and Ericksen,“ yields results in simple shearing 
flow” that are equivalent to those obtained in the theory of simple fluids; 
that is, there are no restrictions on the material functions 7, 0; , and o» for 
Rivlin-Ericksen fluids beyond those which we have stated for simple fluids. 
The greater generality of the theory of simple fluids over that of fluids of 
the differential type is not exhibited in either simple shearing flow or in the 
steady flows which we shall consider in Section 4.* 

There is a famous conjecture, due to Weissenberg,’* to the effect that in 
simple shearing flow one should have 


Su = S33, (3.29) 
which, by 3.18, is equivalent to 
a(x) = 0. (3.30) 


No such result follows from the theory of simple fluids without the addition 
of new special hypotheses. 


We complete our discussion of simple shearing flow by illustrating the | 


solution of an elementary boundary value problem in terms of the material 
functions 7, 0; , and o2. Let us consider a simple shearing flow between an 
infinite plate I at rest and an infinite plate II moving with constant speed V 
parallel to plate I. We denote the distance between the plates by d. The 
a” and x° axes are chosen to be in the plane of plate I and such that plate II 


moves in the 2” direction. The 2! axis is then perpendicular to the plates — 


with x’ = 0 at plate I. Assuming that the fluid adheres to the walls, we 
have the boundary conditions 


v(0) = 0, v(d) = V. (3.31) 


* It is indeed a remarkable fact that each of the steady flow solutions obtained 
by Rivlin" for Rivlin-Ericksen fluids are valid for all simple fluids. The greater gen- 


erality of the theory of simple fluids is best seen by considering not steady flows but 
rather stress relaxation experiments. 
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We further assume that the modified pressure ¢, defined in 2.20, has no 
gradient in the direction of the flow: 


— = 0. (3.32) 


We have pointed out that for simple shearing flow the extra stresses 7';; 
are functions of x alone. Thus by 8.7, the 7';; can vary with the 2’ co- 
ordinate only. It follows directly from this observation that for simple 
shearing flow the dynamical. Equations 2.21 must reduce to 


aT , 0 _ 4 

Ox! azt 

oT 

ae + ss = 0, (3.33) 
Op 
0a3 


| _ A simple analysis of 3.32 and 3.33 shows that these differential equations 


si) i i 


a ek 


are satisfied only if 
T 2 == b (3.34) 
@¢=Tu—-¢, (3.35) 


where 6 and ¢ are constants. 
From 3.34 and 3.17 we observe that « must be independent of x: 


xk = 7(b) = const. (3.36) 


oun 3.7 and the conditions 3.31 then imply that 


v(x!) = (3.37) 
« = V/d. (3.38) 

According to 3.18 the shearing stress is 
| Soe (3.39) 


Since Siz is equal to the tangential force per unit area that must be applied 


- to the moving plate II in order to produce the flow, Equation 3.39 gives an 
~ easily grasped physical meaning to the material function 7. 


It follows from 3.38, 3.17, and 3.16 that all the components 7;; of T 
must be constant. Hence by 3.35, the modified pressure ¢ is a constant. 


Using 2.20 and 2.10 we find : 
Sip Cou py: (3.40) 
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Hence if the potential y of the body forces is known, then a measurement 
of the normal force per unit area on the moving plate determines the value 


of the constant c. 
Equations 3.38 and 3.18 yield 


Su = S33 —= oi(V/d), Soo a S33 ar o3(V/d), (3.41) 


and thus give a physical meaning to the material functions o; and o» . 
On combining Equations 3.39 through 3.41 we find the following expres- 
sion for the matrix of the stress tensor: 


c—py 7(V/d) 0 
| S|] =|} r(V/d) c—py + o2(V/d) —0,(V/d) 0 . (3.42) 
0 0 c—py—oi(V/d) 


4. Some OTHER STEADY FLows 


In the previous section we introduced three material functions 7, o1, 
and a2, which determine the extra stress 7 in simple shearing flow. It is 
clear from Equations 3.39 and 3.41 that if one could carry out in the labora- 
tory a simple shearing flow between two parallel infinite plates then, by 
observing the relative speed of the plates and making appropriate force 
measurements, one could experimentally determine 7, «1, and o2. Of 
course, such an experimental program does not appear practicable at the 
present time. It is a noteworthy fact, however, that in addition to simple 
shearing flow, there are several other steady flow problems that can be 
solved exactly for incompressible simple fluids, and the resulting stress and 
velocity profiles are related by the same three material functions, 7, 0, 
and o2. This class of solvable flow problems contains* Poiseuille flow” 
and also helical flow,” of which concentric pipe flow and Couette flow are 
special cases. If experimental measurements are combined with any one 
of the solutions to determine 7, o; , and a2 , then complete stress and veloc- 
ity profiles can be predicted for the other flows. Experimental rheologists 
do believe that they can achieve in the laboratory and can experimentally 
study flows that are a close approximation to Poiseuille flow, concentric 
pipe flow, and Couette flow. In this section we list, without complete 
derivation, some results of the theory of simple fluids that bear on experi- 
mental studies of these flows. 


_ *Other steady flow problems whose solutions for incompressible simple fluids 
involve only 7, o; , and g¢ are torsional flow, circular flow between concentric spheres, 
and circular flow between a cone and plate. We do not go into details on these flows 
because they are not compatible, under reasonable body forces, with the dynamical 
equations unless inertia is neglected. For the reader who is interested in these par- 
ticular flows, we remark that the solutions of Rivlin, Langlois,” and Markovitz 
and Williamson?® for Rivlin-Ericksen fluids, and the statements made by Markovitz2! — 
about his ‘‘theories of class H,”’ are just as valid for simple fluids. 
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Poiseuille Flow 


We consider a steady flow of an incompressible simple fluid in an infinite 
circular pipe of radius R and use cylindrical coordinates* r, 2, and 6 with 
the z-axis coincident with the axis of the pipe. We say that we have 
Poiseuille flow if in such a coordinate system the contravariant components t 
of the velocity field v(x) = {v’, v*, v’} have the form 


use 0; 
v = v(r), (4.1) 
yo = 0, 
and if the fluid adheres to the walls of the pipe: 
HORS (4.2) 


An important measurable quantity in the experimental approximations 
to Poiseuille flow is the applied force a per unit volume in the direction of 
flow. This quantity a, which is also called the driving force, may be pre- 
cisely defined as follows. Consider the total applied force f in the direction 
of flow exerted on the column of fluid with radius R, cross section @, and 
lying between two planes z = 2; and z = 2y;: 


2II oy 
= = — mea. -- —pdAdz, (43 
f if Sr a] if S fe t @ 02 e é ( ) 


that is, 
j = / [(S.2 = pW emery ih Sz i pv) 21] dA; (4.4) 
@ 
a is then defined as 
Eee? ; 
on wh? (21 — 21) ; ( B 


It can be proved that a is a constant; that is, a is independent of the choice 


ee of zy; and z;. If there are no body forces, then y = const. and a reduces 


to the pressure head per unit length. If the only body force is gravity 
acting in the direction of flow, and if there is no pressure head, then a is 
simply the specific weight of the fluid. 

We shall express the other measurable quantities in terms of the driving 
force a, the radius R, the body force potential y, and the material functions 


T, 01, and op. 


* Throughout this section (Section 4) the symbol r always refers to the radial 
coordinate of a cylindrical coordinate system; in all the other sections of the paper r 
denotes a time lapse. No confusion should arise. 

+ The form of 4.1 is such that the contravariant components of v(x) are the same 
as the physical components, although this is not true for an arbitrary velocity field 


in cylindrical coordinates. 
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It can be shown that in Poiseuille flow the velocity profile is given by 
R 
v(r) = | #Gar) ae, ae 


and from this it follows that the volume discharge A per unit time through 
a cross section of the pipe, 


R 
A= an | v(r)r dr, (4.7) 
0 
is given by 


aR/2 
It = el v7 (v) dp. (4.8) 
0 


The rheological literature contains several derivations, based on special 
assumptions, of formulas equivalent to 4.8;* the point of interest here is 
that 4.8 can be shown to be valid for all incompressible simple fluids. 

It can also be shown that in our cylindrical coordinate system the matrix 
of the physical components of the stress tensor has the form 


See pi Sex he 
[Sl] =||S. S. 0 |}, (4.9) 
O-. Ob See 
where 

Srz = Sx = —}ar, 

gs. =a +p +f talFGa)ld +o (4.10) 
r v 

Srr — See = or [7'(4 ar), 


Sieg oo Soo == (09 [7 (4ar)]. 


Because of the assumed incompressibility of the fluid, when only a, R, 
and y are given, the normal pressures are determined only up to a constant 
hydrostatic pressure of magnitude c. Assuming y to be known, a single 
measurement of S,, fixes the value of this constant c. 

It is worth remarking that according to these equations the normal 
thrust S,. in the axial direction 


S.2 = Sp + o2[7(4ar)] — [7 (ar)] (4.11) 


depends, in general, on r. On comparing 4.11 and 3.28, we see that for 


* See Hermans,”?? pp. 249 and 250. 
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Reiner-Rivlin fluids S,, = S.,. The tendency of a stream of an incom- 
pressible viscoelastic fluid to swell upon efflux from a pipe is most likely 
related to both a dependence of S,, on r and a lack of equality between 
S,, and S,.. Either effect can be used to make the swelling plausible. 


Special Cases of Helical Flow 


Flow between concentric pipes. We consider a steady flow of an incom- 
pressible simple fluid between two fixed coaxial circular cylinders of radii 
Ry, and R, (R, < R.). We again use cylindrical coordinates f,<2, and. 6; 
with the z-axis parallel to the axis of the pipes, and again assume that the 
velocity field has the form 4.1. However, we now impose the following 
boundary conditions: 


UR = o(R:) = 0. (4.12) 


We again use the symbol a for the applied force per unit volume in the 
axial direction: he 5 


x al 
sf a(R — Ry) (21 — 21) 


a (4.13 )* 


Here f is the total applied force in the z direction exerted on the annulus 
of fluid which lies between two planes z = z; and z = 2;;. It can be 
shown that a is independent of z; and 21. 

We wish to express the other measurable quantities in terms of the vari- 
ables a, R, , and R,, y and the functions 7 , 01, and o.. To this end we 
define a function a as follows: 


a(r) = br’ — Lar’; (4.14) 
the constant b is chosen so that 


[Hate ar = 0. (4.16) 


It follows from the known solution to the problem of helical flow” that 
the velocity profile in concentric pipe flow is given by 


v(r) = [ Fla(v)] dy, (4.16) 


1 
and from this it follows that the volume discharge A per unit time, 
Rg 
Av= 2 sf v(r)r dr, (4.17) 
Ry 


* Equation 3.7 of Coleman and Noll contains misprints; it should read as our 


Equation 4.13. 
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is given by 
Re 
were if Hla(r) |r? dr. (4.18) 
Ry 


It can also be shown that in our cylindrical coordinate system the matrix 
of the physical components, the stress tensor for concentric pipe flow, has 
the form 4.9, but now instead of 4.10 we have 


rs = Dis = a(r), 


r 1 2 

Se = az + pb — J ZosfFlal>)} de + o er 
S,, — See = oi{7 la(r)]}, 

Siz — Se = o2{7[a(r)]}. 


When there are no body forces in the radial direction, y is independent 
of r and the second of the equations 4.19 yields an interesting formula for 
the difference in the radial thrusts per unit area on the inner and outer 
cylinders: 


S(Ra) — So(Ri) = — fb oul Plato) a (4.20) 


This difference S,,(R2) — S,,(2,) in the normal stresses on the two cylin- 
ders should certainly be measurable. We note that according to 4.20, 
Sir(Re) = Sy-(R1) if and only if the Weissenberg relation 3.30 holds. It 
appears to us that the measurement of S,,(R:.) — S,,.(R1) in concentric 
pipe flow furnishes the most direct and convenient way available for ex- 
perimenters to check Weissenberg’s conjecture for particular substances. 
We emphasize that there has been no “neglect of inertia” in the derivation 
of 4.20. 


Couette flow. We again consider a steady flow of an imcompressible 
simple fluid between two coaxial cylinders of radii R; and R., (Ri < Re), 
_and we use the same cylindrical coordinate system as before. Instead of 
assuming that the cylinders are at rest, we now postulate that the inner 
and outer cylinders rotate with constant angular velocities Q, and Q», 
fa Seinen The contravariant components of the velocity field v(x) = 


{v", v’, v'} are assumed to have the form 
=e 
v = 0, (4.21) 
v = w(r); 


that is, for Couette flow the motion is in circles. In 4.21 w has the dimen- 
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sions of angular velocity.* Since we assume that the fluid adheres to the 
walls, we have the boundary conditions 
w(R,) =  , w( Rp) = Qs * (4.22) 


An important measurable variable in the experimental approximations 
to Couette flow is the torque M per unit height required to maintain the 
relative motion of the bounding cylinders. 

It can be shown” that in Couette flow the velocity profile obeys the 


equation 
a z (#4): (4.23) 


dr 


thus for this flow, a knowledge of M and the angular velocity of one of the 
bounding cylinders completely determines the velocity profile. Integra- 
tion of 4.23 with the boundary conditions 4.22 yields the following formulat 
relating the observable angular velocity difference 2, — Q, to the torque 
M through the material function r: 
1 M/(27R?) 1 ar : 
Qo i Qy i x. (vy) dp. (4.24) 
2 Juj(2nR?) V 
The matrix of the physical components of the stress tensor for Couette 
flow has the following form in our present coordinate system: 


Srp 0 Syo 
|S|| =| 0 S. O |, (4.25) 
a Se O Seo 
_ where 
M 
Sor ae S,o = Dept: 


3 eit 
4 Sr Sa = a ce ’ 


sec SVE 
‘ Sw — Su = 04 | 4 (25 a: 


Cc 


const. 


* The physical components 0‘ of v(x) are #* = 0" = 0, 0? = Tw. 
a fi Bp ttzations. piaet on special assumptions, of formulas equivalent to 4.23 have 
been known for some time (Mooney,?? Section 3). 


peas br en Ve ies 
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If the potential y of the body forces is independent of 7, the second equa- 
tion of 4.26 yields the following expression for the difference of the normal 
stresses at the outer and inner cylinders:* 


S,r( Re) = S,r( Ri) 


- ff foto fof) tafe)» 


When (R; — R,)/R: < 1, this formula is approximated by 
S,r( Re) a S,r( Ri) 


=" Mo[e(ai)]-fr el): 


The error in this approximation is of the order of [(R: — R1) /Ril. The 
inertial term prw(r) in the integral on the right side of 4.27 makes a con- 
tribution to S,,(Re) — S,r(R1) of the order of [(R2 — Ri) /R,]’, and there- 
fore does not appear in 4.28. 

Padden and DeWitt” have described an apparatus for measuring S,,(R2) 
— §,,(Ri) in Couette flow. On comparing 4.27 and 3.28 we see that for 
Reiner-Rivlin fluids only the centrifugal forces contribute to S,,(R2) — 
S,,(Ri).t Markovitz,” in his analysis of available normal stress data on 
concentrated polyisobutylene solutions of high molecular weight, has re- 
marked that the results of Padden and DeWitt™ show that solutions of 
that polymer do not obey the Reiner-Rivlin theory. In fact Padden and 
DeWitt found that S,,(R.) > S,-(R:), whereas centrifugal effects alone 
would yield S,,(R2) < S,-(R1). 


5. SomE OsctnuAtory Morions 


Here we consider the behavior of incompressible simple fluids in certain 
oscillatory motions. For simplicity, we confine our discussion to a class 
of flows that we call periodic shearing flows. These flows are defined to 
have the property that in some Cartesian coordinate system 2’, 2”, x° the 
velocity field, v(x) = {v', v”, v'}, has the form 


v = 0, 
v = v(a’,t), (5.1) 
v = 0, 


where v is periodic in ¢ with period @; that is, we assume that there exists 
a finite positive number 6 such that 


v(az',t + 0) = v(2", t) (5.2) 
* A result equivalent to 4.27 was obtained by Rivlin for Rivlin-Ericksen fluids 
(Rivlin, p. 186). 


} This fact was pointed out by Rivlin in his early work on Reiner-Rivlin fluids 
(Rivlin,!* p. 279). 
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for all x’ and all t. This flow clearly satisfies Equation 2.18 and hence is 
compatible with incompressibility. 


For a periodic shearing flow, Equations 1.4 and 1.5 take the form 
a , 2 
eS (yes ee ee (5.3) 
dt r 
E(t) = a. (5.4) 


Here the £'(7) are the coordinates of the position at time 7 of the material 
point whose position at time ¢ has coordinates 2’. 
The solution of the system 5.3 with the end condition 5.4 is 


E(t) aie, 
E(r) “3 = oe Uay(2, Sims T); (5.5) 
B(r) = 2°, 


_ where the function w,)(2’, r) is defined by 
un(z,r) = i v(x',t — s) ds. (5.6) 
0 
It follows from 1.6 and 5.5 that, in our present coordinate system, the 


, matrix of the components F(‘;(t — r) of the deformation gradient tensor 
BF i)(t — r) is 


il 0 O 
A lFot —) l= lFeitt—o ll =P) 1 01, Gz 
. C01 
3 Aw(r) = ae (5.8) 


2 Since we shall henceforth keep our attention on a particular point in space, 
a we do not need to make explicit the dependence of \(1)(7) on x. 

- We note, for future use, that 5.2 implies that the following equation 
holds for all rand #: 


; Azo (7) = Aww (7); (5.9) 


that is, for each r, Ai) (7), regarded as a function of t, is periodic with period 0. 
a By 5.7 and 1.9 the matrix of the components of C(t — r) is 


7 L+AG(r) Aw(r)_ 0 
g Ck jr = Aw (7) i OAK (5.10) 
0 0 1 
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Hence in all periodic shearing flows, C:(t — r) has the form 
Ci(t = r) rf + E,(r) + G(r) (5.11) 


where J is again the unit tensor and the matrices of the components of 
E,(r) and G(r) are 


Dy a te 
E,(r) = ho@tr) | Y 0 OF, 
(Chea Ue (8. 
(5.12) 
be0.0 
G.i(r) = Nw(r) || 0 0 0 
00470 


It is obvious that whenever the history C:(¢ — r) has the form 5.11, 
C(t — r) is determined by two functions of r, EZ; and G;. When this is 
the case, it follows from 2.9 that we may regard the extra stress at time ¢ 
as given by a functional G whose arguments are the functions H, and G; . 


T(t) = © (Bur), Gr)) (6.18) 


Furthermore, it follows from 2.13 and 5.11 that G must be an isotropic 
functional; that is, the following identity must hold for all functions FE; , 
G, and all constant orthogonal tensors Q: 


QG Bur), G(r)1Q" = 8 (QBAr)Q", QGUrIQ". (B.A) 


Let us consider the particular orthogonal tensor Q whose components, 
in our present Cartesian coordinate system, are given by the matrix 


Ped Go 
|@=]0 1 O fT. (5 15) 
00 -1 | 
Then 5.12 yields 
QE(r)Q" = E(r), QGi(r)Q” = G(r). (5.16) 
Hence, from 5.13 and 5.14 we get 
QT (t)Q" = T(t). ‘ (5.17) 
Denoting the matrix of the components of 7 by 
Tu Ti. Ts 
IP = Der Tee” Toe ||, (5.18) 


T31 T's2 T'33 
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we find that the matrix of the components of Q7Q" is given by 
Tu Ty —T 
OT w= dee ash Tae = TI: (5.19) 
—T3 —Ts T'33 


According to 5.17, the matrices exhibited in 5.18 and 5.19 must coincide; 
this is possible only if Ti; = Tx = To = Ts. = 0. Hence we have proved 
that in periodic shearing flow the matrix of the components of the extra 
stress has the form 


T(t) T12(t) 0 


TG) =|! Ta Gy) *Poe(t) 0 : (5.20) 
0 0 T'33(t) 
From 2.11 we have 
Tu(t) + Too(t) + T33(t) = 0. (5.21) 


It follows from 6.12 and 5.18 that the components T.;(t) of T(t) depend 

on only the function \;). Furthermore, it is clear from 5.20 and 5.21 that 

these components 7';;(t) may be expressed in terms of three independent 
functionals with argument \(). We choose 


io} 


P(t) = ft Do], 
T(t) — T(t) = & Deo(r)h (6.22) 


Tas(t) — Tall) = 8 Deol?) 


In consequence of 5.14 the three functionals t, 8; , and 8 are independent 
of the orientation of our Cartesian coordinate system and depend only 
; on the material under consideration. We call t, 8; ,. and. 8, material func- 
 tionals. 

A The isotropy relation 5,14 places restrictions on t, 8, and &. To find 
_ these restrictions we take 


2 ea Oi 

a Qi} =\}0 -1 0 (5.23) 
z DinrurOihest 

°$ 


iS as a second choice of Q in 5.14. The Equations 5.12 then yield 


ae, 


Z QE(r)Q” = —Ei(r),  QGi(r)Q? = Gir), (5.24) 


i. 
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and from 5.20 we get 
Tu(t) —T a(t) 0 
QT(t)Q7 || = || -—Tw(t) T(t) O |}. (5.25) 
0 0 T'33(t) 


It is clear from 5.12 that the conjugation operation 5.24 is equivalent to 
replacing Ai (7) by —\y(r). Thus by 5.14 and 5.13, when we replace the 
function Xj by the function —Ai , the matrix \| T(t) || of Equation 5.20 
goes over into the matrix || QT (t)Q’ || of Equation 5.25. This is equiva- 
lent to saying that the material functionals of 5.22 obey the following 
identities: 


t [-rwlr)] =a e Dwr},  si—dAw(r)] = &PPw(7)]; (5.26) 
that is, that t is an odd functional, while & and 8 are even functionals. 

Whenever the matrix of the Cartesian components of the history C(t —1r) 
has the form shown in 5.10, the Equations 5.20, 5.22, and 5.26 are not 
only necessary, but are also sufficient, for the validity of Equation 2.13 
for all constant orthogonal tensors Q. Hence the isotropy condition 2.13 
(or 5.14) yields no further restrictions on the functionals t, 8 , and 82 be- 
yond those exhibited in Equations 5.26. 

The results obtained so far are valid for arbitrary, not necessarily pe- 
riodic, flows of the type 5.1. We can now examine the consequences of 
the periodicity 5.2 of the velocity field. 

For definiteness, let us assume that @ is the smallest positive number 
such that 5.2 is true for all 2’ and all ¢. Equations 5.9 and 5.12 then yield 
E.so(r) = Ex(r) and Gi40(r) = Gi(r) for all r and ¢, and it then follows 
from 5.13 that 


T(t + 0) = T(t) (5.27) 


for all ¢. Thus we have the not unexpected result that the extra stress is 


periodic with period 6. This does not mean, however, that @ need be the - 


minimum period for each of the components of T(t). It turns out that 
by adding to 5.2 an assumption that is usually met in experimental approxi- 
mations to periodic shearing flow and that by using 5.26, we get a result 
much stronger than 5.27. The assumption is that 


v(2', t) 7 —v(a', t+ 0/2); (5.28) 


for all x’ and t; that is, that the absolute value of v has a period equal to 
one half the minimum period of v.* 


* This assumption is obeyed, for example, by v(z!,t) = 6(2!) sin [= + (a) |. 
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It follows from 5.28 and the definition of Aw(r) that 


Aw(T) = —Nepo/y(1) (5.29) 
for all r and all ¢. Equations 5.29 and 5.26 yield 


£ Pusm(r)] = t [dol] = —t Pol) 


(5.30) 


io} 


8 [Acepo/2(7)] = & [—Aw(r)] 


& co(r)], 
7 = 1, 2; these equations hold identically in ¢. It immediately follows 
from 5.30, 5.22, and 5.21 that the following equations hold for all t: 

T(t + 6/2) = —T p(t), 

Tu(t + 6/2) = T(t), 

Pea (t 0 02 =n (2), 

Pas(t + 0/2) = Toad). 


Hence we have proved that if the velocity has the property 5.28 then the 
shearing stress T, also has this property and, furthermore, the normal 
extra stresses Ty, , Tx, and 73 oscillate with a frequency that is twice 
the frequency of oscillation of the velocity field and the displacement 
field.* 


I 


(5.31) 


6. Stow Motions anp NEwTOonIAN FLuuIps 


7 We fix the present time ¢ and use the abbreviation C(r) = C,(t — r) 
~ to emphasize the fact that the history of the right Cauchy-Green tensor 
is a function of r. Only the dependence on r matters in determining the 
value of the functional § which gives the stress. 

In Section 27 we pointed out that if a simple fluid has been at rest for 
© all times 7, —«~ < 7 < #, then the history C(r) = C(t — r), regarded 
as a function of r, reduces to that function whose value for each r is the 
_ identity tensor I: 

Coy = CS 71Ph= 7, ONS Perey (6.1) 


f * We remark that the only essential constitutive assumption required for the proof 
_ of 5.26, and thus 5.81, is the isotropy condition 2.13. In fact, results similar to 5.31 
_ can be derived for the Cartesian components of the stress tensor S in a periodic shear- 
’ ing motion in an arbitrary isotropic body that is not necessarily either a fluid or 
_ incompressible, provided only that the motion is a symmetric oscillation about an 
undistorted configuration. We further remark that just as the results of Section 3 for 
_ simple shearing flow can be shown to hold for flows that have a similar form in other 
orthogonal coordinate systems,!° so can our results on periodic shearing flow be ex- 
_ tended to flows which have the form 5.1 and 5.2 in these coordinate systems. How- 
ever, the treatment of boundary value problems for periodic flows is a much more 
- complicated matter than the treatment of such problems for steady flows. ve 
{ See Equations 2.3 through 2.8. We now drop the assumption of incompressibility: 


made in Sections 3 through 5. 
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We then showed that when 6.1 holds the stress must reduce to a hydrostatic 
pressure which depends on the density alone: 


ees (6.2) 


The dynamical theory of perfect fluids deals with those special examples 
of simple fluids that have the property that Equation 6.2 holds for all 
histories C(r) = C.(t — 1); that is, for perfect fluids the stress is always a 
hydrostatic pressure, even if the fluid has been sheared in the past or is 
being sheared at present. 

Since we know that 6.2 holds with exactitude for all simple fluids that 
have been eternally at rest, we feel that it must be possible to make precise 
the intuitive notion that if we restrict our considerations to histories that 
do not differ too much from that exhibited in 6.1, then the theory of per- 
fect fluids should be, in some sense, an approximation to the general theory 
of simple fluids. 

The theory of Newtonian fluids is based on the following constitutive 
equation for the stress: 


S = —p(p)I + [A(e)trDZ + 2n()D (6.3) 


here D is one half of the first Rivlin-Ericksen tensor A, or, by 1.18, the 
symmetric part of the velocity gradient tensor;* 7 and 2 are functions of 
p alone and are called coefficients of viscosity ;t \ plays no role in isochorie 
motions, for in such motions trD = 0. When a fluid is at rest we have 
D = 0. Thus 6.3 reduces to 6.2 whenever the fluid is at. rest, and the 
physical significance of the material function p in 6.3 is clear: p(p) is the 
hydrostatic pressure the fluid would be supporting if it were at rest at its 
present density. 

It is clear that all materials obeying Equation 6.3 are simple fluids. 
Furthermore it can be shown, using the isotropy condition 2.2, that Equa- 
tion 6.3 holds for any simple fluid for which the stress depends only on p 
and the velocity gradient, with a linear dependence on the latter. 

For some substances, such as water, the theory of Newtonian fluids ac- 
counts remarkably well for experimental measurements over a very wide 
range of conditions. Other substances, such as molten polymers and 
polymer solutions, definitely do not obey Equation 6.3 exactly. However, 
our intuition and a host of rheological measurements suggest that the 
theory of Newtonian fluids should approximate the behavior of nearly all 
real fluids in the limit of slow motions. In other words, our physical 

* D is also called the rate of deformation tensor. Truesdell, p. 126, calls 6.3 
the Newton-Cauchy-Poisson law and gives a brief summary of its history. 

t n(p) is sometimes called the shear viscosity and SMe) t Auf) the bulk viscosity. 


{ In essence this was first shown by Stokes; it is in fact an immediate consequence 
of the now famous representation theorem for isotropic tensor functions of one tensor 


pene For a recent direct proof the reader is referred to Serrin’s encyclopedia 
article. 


i heel at Ys . + 


mA 
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prejudice indicates that the logical status of the theory of Newtonian fluids 
in the theory of simple fluids should be this: when the kinematical history 
does not differ too much from the “rest history” shown in 6.1, the theory 
of Newtonian fluids should be a first-order correction to the theory of per- 
fect fluids. 

There is more than one possible mathematical interpretation for the 
statement that the theories of perfect and Newtonian fluids are “approxi- 
mations” to the more general theory of simple fluids. The particular 
interpretation. that interests us here may be roughly outlined as follows. 
We first characterize the various possible local kinematical histories at a 
material point by the tensor functions C(r) = C(t — r), regarded as 
functions of 7, and the densities p = p(t), for these are the only variables 
that enter the theory of simple fluids. We restrict our attention to a 
particular class of the functions C'(r); this is the class of functions that are, 
in some sense to be made precise later, “close to” the function shown in 
6.1. We assume that for each p(t) the functional © in the general con- 
stitutive Equation 2.1 is, again in some sense to be made precise later, a 
“smooth” functional at the function 6.1. We then show that if for any 
p(t) we consider functions C(r) that “approach” the identity function 
C(r) = I, then the values of the stress S(t) computed using the functional 
§ will approach those computed using the constitutive equations of per- 
fect and Newtonian fluids. To show that the Newtonian fluid is a com- 
plete first-order correction to the theory of perfect fluids we must show 
that as the functions C(7) approach the identity function C(r) = I the 
error made in using 6.3 in place of 2.1 goes to zero faster than the term 
(p)trD + 2nD that appears in 6.3 but not in 6.2.* 

The crucial issue to be settled before we can develop these ideas is that 
of making precise what we mean by the distance between two histories of 
Cauchy-Green tensors. We must somehow introduce a notion of distance, 
or a metric, into a space § of functions C(r) whose values are symmetric 
tensors and whose arguments are numbers r, 0 <r < o. We shall do 
this by defining a norm for the functions in $; the norm of each function 


~ can be thought of as its “magnitude” or “length.” The distance between 


two functions C’(r) and C®(r) will then be the norm, or length, of their 
difference OC (r) — C®(r). 

There is one essential property that we demand of this norm. It must 
be compatible with the fact that all real materials have an imperfect mem- 
ory: in any experiment we expect the events that happened in the recent 
past to be more important than those that happened in the very distant 

* The reader will notice that we have picked the simplest possible interpretation 


that can be given to the notion that one physical theory approximates another. 
Our discussion is limited to constitutive equations. We do not ask whether the solu- 


tions of given boundary value problems are similar in the different theories, for this 


at present appears to be an extremely difficult question. Except for extremely trivial 
ohven we do, ae even know what is a well-put boundary value problem for simple 


fluids. , Z 
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past. Indeed if this were not the case it would be exceedingly difficult to 
design experiments, because the experimenter, not being in control of events 
that happened before he was born, would never be sure that he had ade- 
quate knowledge of the history of the materials he was testing. Since 
physical prejudice suggests that there is a wide variation in the memories 
exhibited by different substances, we might expect the actual recipe to be 
used for computing the norm of a function C (r) should vary with the sub- 
stance under consideration. Any recipe used, however, must be compatible 
with the statement that all memories are imperfect and must, therefore, 
place greater emphasis on the values of C(r) for small 7 than the values 
for very large r. In the remainder of this section we shall consider a class 
of norms with such a property. By assuming that for the metric generated 
by some norm in this class the constitutive functional of each simple 
fluid is differentiable at the identity function 6.1, we shall be able to obtain 
theorems that tend to justify our speculations about the approximate 
behavior of simple fluids in slow motions. 


An Approximation Theorem 


We deal here with symmetric tensors only. We follow standard practice 
and define the norm | A | of a symmetric tensor A to be 


| A | = Vtr(A?). (6.4) 
It is easily verified that 
if A.-= 0, then |A|= 0; (6.5) 
if A #0, then|A|> 0; (6.6) 
and 
|A+B|<|A|+ |B (6.7) 


for all symmetric tensors A and B. Equations 6.5 to 6.7 justify our use 
of the word norm. 

Let us suppose that we can associate with any simple fluid a real-valued 
function h of a real variable r,0 < r < ». We call h(r) the influence 
function for the fluid. We shall assume that h(r) obeys the following 
three conditions: (a) h(r) > 0 for all r; (b) A(r) is continuous over its 
entire domain; and (c) for all positive numbers c the function 


h 
g(a,c) = sup ie, a>Q, (68)* 
satisfies 
. gla,c 
ar ie pe 


* The term sup stands for supremum of, that is, least upper bound of. 
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In rough language, we can say that the condition (c) states that h(r) 
approaches zero very rapidly as r approaches infinity. For example, the 
condition (c) is not obeyed by functions of the form ; 


hr) = +r)” (6.10) 


if m < 1, but is obeyed by such functions if m > 1. The condition is also 
obeyed by all functions of the form 


hiry =" + re", 8 > 0, marbitrary. (6.11) 


Tn fact, it can be proved that our hypotheses on h(r) are obeyed by any 
positive continuous function which for some e > 0 obeys the limit relation 


lim r'*h(r) = 0. 
monotonically for large r. 


Let A(r) be a (symmetric) tensor-valued function of the real variable 
r,O <r < . We define the norm | A(r) |, of the function A(r) to be 


PAtr ly = sup { | A(r) | h(r)}. (6.12)* 

We restrict our attention to the set $, of tensor-valued functions A (r) 
such that 

| A(r) |, < ©. (6.13) 


The set $; obviously depends on the choice of the function h(r). It follows 
from Equations 6.5 to 6.7 and the definition 6.12 that 


if A(r) = 0, then|A(r) |, = 0; (6.14) 
if A(r) ¥ 0 for some value of r, then|A(r) |, > 0; (6.15) 
and 
i |A(r) + Br) |p < [AC [ns + | BO) [a (6.16) 


for all functions A(r), B(r) in S&. Thus S, is a normed function space 
and has a metric d; defined by 


dx{A(r), B(r)} = | A(r) — Br) ha (6.17) 


d,{ A(r), B(r)} may be called the distance between the functions A(r) and 
B(r). 

Let & be a functional whose domain is the normed function space §&, 
and whose range is the set of all symmetric tensors. We say that the 


* The subscript h on | A(r) |, is a reminder that the norm is defined over a func- 
tion space and that the norm depends on the influence function h. 
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functional § possesses a Fréchet differential* 5§ at the function A(r) if 
there exists a functional 5§ such that the following equation holds for all 
functions B(r) in 8; : 


(Ar) + BOY = & An) 


(6.18) 


co 


+ 8 (4A (r); B(r)] + | B(r) bs BBO) 
r=0 | 
where, when | B(r) |, goes to zero, the remainder term | B(r) |, Reo 
[B(r)] approaches zero faster than | B(r) |, , that is, 


lm « R B))= 0: (6.19) 
|B(r) [4,20 r=0 
the functional 5§ is assumed to be both linear and continuous in its second 
variable. Here, linearity means that 


8§ [A(r); a B(r) + a®B(r)] 
ae if . (6.20) 
= a5 (A(r)s B(r)] + a9 (A (7), B(P)] 


holds for all numbers a”, a and for all functions B® (r), B® (r) in $8 ; 
continuity means that 


i] 


aR (A(r)s BO(H)] — 4 [4(r); BOY)]] = 0 (6.24) 


lim 
d;,{B“*)(r),B(r)}+0 | r=0 
for all functions B(r) in $, and for all sequences B® (r) (of functions in 
S,) that approach B(r). 

We are now ready to lay down our fundamental smoothness assumption 
about the constitutive functionals © occurring in the definition of a simple 
fluid. We assume that for each simple fluid there is an influence function 
h(r) obeying the conditions (a), (b), and (¢) and such that in the normed 
function space $, the constitutive functional § possesses a Fréchet differ- 
ential 6 at the identity function 6.1. This means that we can put for 
any history C.(t — r) in 8, 


Bir) =I -—Ci(t — 1), (6.22) 
and the stress S(t) corresponding to C,(t — r) will be given by 


i} 


S() = G [Ct = 1) 010] = 9s OI + 88 UF; Br); (0) 


r=0 


r (6.23) 
+ | Br) |n ® (Br); o(2)] 


* The theory of Fréchet differentials of functionals is developed in detail in Hille 
and Phillips.?6 ’ 


Ree oe Neti w ba Yor he ani 
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where the functional obeys 6.19 and the functional 6 is linear and con- 
tinuous in its second variable, Bre 

From the discussion of the paragraph containing Equations 2.3 through 
2.8, we see that the tensor 7. [J; p(t)] has the form 


SE ol] = — ploy (6 24) 


where p[p(t)] is the hydrostatic pressure that the material point under 
consideration would now be sustaining if it had been at rest with the den- 
sity p(t) at all times 7, —~0 < 7 < ¢. The term in 6.23 involving the 
Fréchet differential 6 gives, in a certain sense, a continuous linear correc- 
tion to the value of at the function 6.1. 

Suppose we are given a history C(r) in 8, ; we wish to have for compari- 
son a set of other histories that may be regarded as being “essentially the 
same” as C'(r) yet having been “carried out at a slower rate.” Perhaps 
the simplest way to do this is to replace the variable r by a new variable 
ar, where the range of the parameter a is 0 < a < 1. We shall denote 
the new functions of r generated this way by C(ar): 


C(ar) = Ci(t — ra). (6.25 ) 


We call the functions C'(ar) the retarded histories corresponding to C(r) = 
Ci(t = T)) 

It is evident from 1.10 that for a = 0, C(ar) reduces to the identity 
function 6.1. We are here interested in the way the retarded histories 
approach the identity function as a approaches zero. We assume that 
C(r) is differentiable at r = 0. It can then be proved, using the conditions 
on the influence function h, that for small a the retarded histories C(ar), 
corresponding to C(r) = C.,(t — r), are “close to” the function 6.1 in the 
following sense: 


C(ar) = I — 2aD + o(a,7r), (6.26) 
where o(a, 7) is a tensor-valued function of r in the space 8, , and 
ah Oe barerg (6.27) 
a0 a 


In 6.26 D is the symmetric part of the velocity gradient tensor 1.18, 1.13: 
d d 
= een i = : 6.28 
D=-—FC(r)| = -Fekt—1)|_ (6.28) 


Although the result shown in 6.26 is not unexpected, its proof (which we 
do not go into here) is nontrivial and depends very heavily on our assump- 
tion of condition (c) on h(r): Equation 6.26 would not hold if the influence 
function h(r) did not go to zero rapidly enough asr > 0, ~ 
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Once 6.26 is established it is relatively easy to prove the following 
theorem, which is a direct consequence of the definition of a simple fluid 
2.1, 2.2, the smoothness assumption 6.23, and the ‘‘lemma”’ 6.26. 


Theorem. Let C(r) be a function in §, that is differentiable at r = 0. 
For the retarded histories C(ar), corresponding to C(r) = C.(t — 1), the 
functional $, which gives the stress in a simple fluid, takes the following 
form: 


6 (Clore Te an ee (6.29) 


r=0 
Here p(p) is given by 6.24; aD is given by 
_1dC(t — ar) _ _adC: ¢ — r) 
2 


cE rs dr r=0 2 dr r=0 


- (6.30) 


G(aD; p) is a (symmetric) tensor-valued function that is linear and iso- 
tropic in aD; and o(a) is a tensor-valued function of a with the property 
that 


= 0. (6.31) 


The statement that the function G(aD; p) is linear in aD means that 
for all numbers a” and @” and for all symmetric tensors D" and D®, 
Ga? Dp” a Chen Meee p) - aG(D”: p) mt a” G(D”; p). (6.32) 
The isotropy of G(aD; p) in aD means that 
G(aQDQ"; p) = QG(aD; p)Q’ (6.33) 


for all symmetric tensors D and all orthogonal tensors Q. It can be shown* 
that these conditions of linearity and isotropy are equivalent to the asser- 
tion that G(aD; p) has the form 


G(aD; p) = [A(p)tr(aD)] + 2n(p)aD. (6.34) 

Hence our theorem states that the stress at time ¢ in a simple fluid with a 
history of the form C;(t — ra) is given to within terms of order a by 

S(t) = —plo(t)|l + {No(t)]tr(aD)}Z + n[p(t)]aD + o(a). (6.35) 


The first term on the right in 6.35, —p[p(t)|Z, gives the stress that the 
fluid would be bearing at time ¢ if it were a perfect fluid. We note that 
by 6.30 and 1.13, aD is one half the first Rivlin-Ericksen tensor (or, by 
1.18, the symmetric part of the velocity gradient at time ¢) corresponding 
to the history C(t — ra). Hence on comparing 6.35 with 6.3 we see that 


* See, for example, the proof in Serrin,?® which we have previously mentioned. 
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the first three terms on the right in 6.35 give the stress the fluid would be 
bearing if it were a Newtonian fluid. The remaining term 0(a) is of higher 
order in @ than the first three terms. Since a is a measure of the “rate of 
motion” with ‘“‘slow motions” corresponding to small a, these observations 
show that the theory of the Newtonian fluid is indeed a complete first-order 
correction to the theory of perfect fluids in the limit of slow motions. 


7. HiGHER-ORDER CoRRECTIONS TO THE THEORY OF PERFECT FLuIps 


In the previous section we showed that if we replace a given history 
of the Cauchy-Green tensor C,(t — r) by a retarded history C(ar) = 
C(t — ra), 0 < a < 1, then in the limit as a approaches zero the con- 
stitutive equation for the stress in a simple fluid agrees, to within terms of 
order one in a, with the constitutive equation of Newtonian fluids. Here 
we shall sketch a rigorous procedure for ascertaining the form of terms of 
higher order ina. At the outset we remark that it appears to us that any 
approximation theorem exhibiting these higher-order terms must rest on 
hypotheses that are stronger than the smoothness assumption laid down 
in Section 6. First, it turns out that we need a condition on the influence 
function that is stronger than our previous condition (c); second, it is 
obvious that we must now work with Fréchet differentials of order higher 
than one. 


Approximations of Order n,n > 1 


We again restrict our considerations to symmetric tensors and use Equa- 
tion 6.4 as the definition of the norm of a tensor. 

We continue to suppose that we can associate with any simple fluid a 
real-valued function h of the real variable 7,0 <r < «©. We now assume 
that this influence function h(r) obeys the following three conditions: 
(a’) h(r) > 0 for all r; (b’) h(7) is continuous over its entire domain; and 
(e’) for all positive numbers c the function 


g(a,e) = sup ae, xa 
satisfies 
Him 26% 2) = 9, (7.2) 


a>0 a” 


The conditions (a’) and (b’) are repetitions of the conditions (a) and (b) 
assumed for f(r) in Section 6; however, since n is here a positive integer 
greater than 1, the condition (c’) is stronger than the condition (c) of 
Section 6. We note that the condition (c’) is still obeyed by all functions 
of the form 6.11 and by those functions of the form 6.10 for which m > n. 
Our new hypotheses on h(r), though obeyed by a smaller class of functions 
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h(r) than the hypotheses made in Section 6, are still obeyed for a given n 
by any positive continuous function for which there exists an « > 0 such 
that the limit relation 

lim r"*‘h(r) = 0 
holds monotonically for large r. Furthermore, condition (c’) is obeyed 
for all n if h(r) approaches zero exponentially fast as r > ». 

We again use Equation 6.12 to define the norm of tensor-valued functions 
of r, and we restrict our attention to the set 8, of functions for which this 
norm is finite. 

Consider a functional whose domain is S, and whose range is the set 
of all symmetric tensors. We say that § has an nth order Fréchet differ- 
ential at the function A(r) if there exist n functionals 8’§ (of j + 1 function 
variables) such that the following equation holds for all functions B(r) in 
Sh: 


BLA) + BO) = BAO + Y 585 MAMBO), --- BO) 


r=0 


(7.3) 


L] 


+ [/B(r) hl” & (BC), 


where the functional §t obeys Equation 6.19 and, in addition, the jth func- 
tional 6’, when regarded as a functional of 7 + 1 distinct function vari- 
ables 


5 [A(r); Ba(r), -+> Buy (7), (7.4) 


is not only completely symmetric* and jointly continuous in its last 7 
variables Ba)(r), -++ Byj;(r), but is also linear in each of these 7 variables 
when all of the other variables are held fixed.t It then follows that each 
functional 6’§ is uniquely determined. Ifthe j variables Bay(r), --- By;)(r) 
are all the same, as they are in 7.3, that is, 


Bw(r) = B(r), b= 1 


then 6’§, regarded as a functional of B(r), is a homogeneous functional 
of order j. The functional 5’§ in 7.3 and 7.4 is called the jth order Fréchet 
differential of the functional §. 

We now have the mathematical apparatus required for asserting a 
smoothness assumption about the constitutive functionals § for simple 
fluids that is strong enough to enable us to obtain corrections to the theory 
of perfect fluids that are complete up to terms of order n in a. We call 

: ; ; ; ; 

Gara vathous changtig the valua Oy Los HONS OEET Se te oa aa 


_ ¢ For short, we say that the functional 67 ; : +: coke : 
waaticgs ional 6/§} shown in7.4 is multilinear in its last j 
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this assumption the augmented smoothness assumption of order n.° It as- 
serts that for each simple fluid there is an influence function h(r) obeying 
the conditions (a’), (b’), and (c’) and such that in.the normed function 
space §) the constitutive functional © has an nth order Fréchet differential 
at the identity function 6.1. In other words, our augmented smoothness 
assumption asserts that if, for any history C:(t — r) in 8,, we put Ci(t — 
r) = I + B(r) then 


S() = $ [lt 7); 0101 = $ U5 010] 
+ 259 U5 BO), + BO); a) (7.5) 


+ [| BG) hl? % (BO) 


where the functional 9 obeys 6.19 and the functional 5’, when regarded 
as a functional of 7 + 1 distinct function variables 


oD Par) iss Belts pl, (7.6) 
is multilinear, continuous, and completely symmetric in its last 7 variables, 
Bw(r) --- Boy(r), and is uniquely determined by ©. 

We are now interested again in the way retarded histories C (ar), corre- 
sponding to a given history C(r) = C.(t — r) in 8, , approach the identity 
function 6.1 as a approaches zero. The result 6.26 is not strong enough 
for our present purposes. The new condition (c’) on the influence func- 
tion h(r) now enables us, however, to prove the following statement: 
Assume that C(1r) is n times differentiable at r = 0, that is, that the kth 
_ Rivlin-Ericksen tensor corresponding to the history C(r) 


hig 


— _ Fal ‘ 
A, = (—1) a C(r) de (7.7) 
_ exists for k = 1, --- n; it then follows that 
n k 
C(ar) =I+)>> ( tt Az + ofa"; r) (7.8) 
k=1 : 


_ where 0(a”; r) is a tensor-valued function of r in §, with the property 


lim bokerehh =O}: (7.9) 


We omit the proof of 7.8 but remark that it rests heavily on our assumption 
_ of condition (c’) on h(r). 

_ The definition of a simple fluid 2.1 and 2.2, the augmented smoothness 
_ assumption of order n 7.5, and the new “lemma” 7.8 yield the following 


_ theorem without much difficulty. 
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Theorem. Let C(r) be a function in §&; , which is n times differentiable 
at r = 0. For the retarded histories C(ar) the constitutive functional 5 
takes the following form: 


© [C(ar); pl 
- : a ° 7.10 
= —plp)l $2... Gis ee ae o(a"). ag 
jal kyt---kj<n 
ky<ko<-+ +k; 
Here 0(a”) is tensor-valued function of a such that 
lim Jota") | = 0; (7.11) 
a0 (ed 
A,” is the kth Rivlin-Ericksen tensor for the retarded history C (ar), 
k d*C(ar) k aCiAt = r) k 
ee = (-—a)* —_—— = & Axs 12 
Ax ( 1) dr* r=0 ( a) dr* r=0 eae a ) 
GPR Ak a aes p), asa function of the 7 symmetric tensors Aj, , + -* 


A;,, is an isotropic, multilinear function whose values are symmetric 
tensors; and p(p) is given by 6.24. 

The statement that Gj,x,,...4; 18 an isotropic function means that the 
following identity holds for all orthogonal tensors Q and all sets of 7 sym- 
metric tensors A;,, ++ Ax; : 


QG j;ky,..4j(Any eae Mi A; ) p)Q” 
= Gin, ...bj(QAQ'; ae QA;,Q’; p)- 


When we say that Gj,x,,....; 18 a multilinear function we mean that if any 
j — 1 of the tensor variables A;,, --- Ax; are held fixed, then Gisky , te 
is linear, in the sense of 6.32, in the remaining variable. It will be noticed 
that the term Gj,2,,...4;,..4j(A%, , °°* Ak; , *+* Ax,) is, therefore, of order 
Paes k; in a. 

Perhaps it will help clarify our notation if we write out in detail Equation 
7.10 for the special case of n = 3: 


(7.13) 


© (Clan); e) = —p(p)I + Gir(Ar*; p) + Gi;2(A2"; p) 


+ Gij3(As"; p) + Gei1a(Ar", Ar*; p) (7.14) 
+ Goj1,2( Ai", Ae”; p) + -Gsy14a(Ar", Ar”, Ai; p)> 
+ 0(a°). 


To make clear the dependence of each term on a, 7.14 may be rewritten as 
follows: 


PRE ee ee 


ee ON Ck ee eb 


Coleman & Noll: Viscoelastic Fluids (it 


2 lCCar); ol = —ple)T + aGia( Ar; p) + a[Gr2(As ; p) 


F Go11(A1, Ar; p)] + a[G43(As ; p) (7.15) 
+ Gor 2(Ar, As; p) + Gsara(A1, 41, Ar; p)] 
+ ofar). 

Each function Gj;x,,...4 Hage - Ax;), ki < ke < +++ kj, because it is 


isotropic and multilinear, may ~ eprestadl as a sum su products of the 
form 


A(p)dide +++ bm(AnArn, +++ Ai, + Ay -*- AnAn) (7.16) 


where the g:8 are traces of products of some of the tensors A x, and are such 
that each A;,, occurs precisely once in 7.16.* Hence it follows from our 
theorem that for each finite n only a finite number of scalar material func- 
tions A;(p) are needed to determine the functional to within terms of 
order nina. For incompressible materials these material functions reduce 
to constants. Thus a finite number of material constants fix the constitu- 
tive equation for the extra stress in an incompressible fluid to within terms 
of order n in a. We illustrate these remarks by considering the special 
case of n = 2. 


Approximations of Order 2 
In this case Equation 7.10 becomes 
 [C(ar); p] = —p(p)I + Gir(A1*; p) + Gij2(A2"; p) 
r=0 (7.17) 
+ Ge;11(A1", Ai", p) +.0(0"). 


The tensor functions G4;1(A1 ; p), Gi,2(A2 3 p) and Gy11(A1*, A1*; p) must 
have the form 


Gin (Ax > p) = [AytrAj|7 + 2A, (7.18) 
Gi,2( Ao ; p) = [AstrAs]J + AsA2 (7.19) 


: Go. (Ay, Ay 4 p) = [As(trA;)? + AgtrAy [I + [A7trA JA, + NA: (7.20) 


The coefficients \; in Equations 7.18 to 7.20 depend on only the density p. 
It is now clear that Equations 7.17 through 7.20 yield the following expres- 
sion for the stress S(t) in a compressible simple fluid with the history 
Cit a ra): 


S(t) = pl + wy [trA,*| = H2A,* + [ustr(A,*)" + pa(trA,*)? 
+ pstrAs"\Z + [uctrAr“JA1* + v7(Ar*)” + psd2* + 0(0°) 


* Methods for further reduction of expressions of the form 7.16 have been de- 
veloped by Spencer and Rivlin.?7,8 


(7.21) 
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where the coefficients y; , being just the coefficients A; numbered in a more 
convenient order, are functions of p(t) alone. The terms involving m and 
uy. are of order lin a. Since A,* = 2aD, on comparing 7.21 with 6.35 we 


see that 


Liga Vor(p) 
2 = n(p). 
The terms involving y3, --* us are of order two in a. 


Note: The theory of Reiner-Rivlin fluids would not yield the terms with 
us and us. Hence from our present point of view, the Reiner-Rivlin 
theory does not yield for simple fluids a complete first correction to the 
theory of Newtonian fluids. 

The second-order approximation to the constitutive equation 2.9 of an 
incompressible simple fluid is obtained from 7.21 by the following modifi- 
cations. First of all, we recall that in an incompressible fluid the pressure 
p is indeterminate and not given by 6.24. The constitutive equation gives 
only the extra stress T(t). Furthermore, in the expression for the extra 
stress in an incompressible fluid all of the terms in 7.21 that are scalar 
products of the identity tensor J will be determined by the normalization 
yap 


trT = 0. (7.22) 


Also, it follows from 2.18 that in any permissible motion in an incompres- 
sible fluid we must have 


tr(A,*) = 0; (7.23) 


this means that the constitutive equation for such a fluid will have no 
terms such as those involving 1, ws, and ys in 7.21. Finally we note that 
the remaining coefficients (that is, those corresponding to ue, uz, and ys) 
must now be constants. In this way we obtain the following simple for- 
mula for the extra stress T(t) in an incompressible fluid with a history 
Cit — ra): 
T(t) = yf + BiAi* + Bo(A1*)” + BsA2* + o(a’), (7.24) 
where 8; , 62 , and 8; are material constants and y is determined by 7.22, 
y = —VWl[Botr(Ai*)* + BstrAe’l. (7.25) 


It is useful here to use a normalization different from 7.22 by absorbing 
the term J of 7.24 into the extra stress. We then have 


T = S+ pl = BA," + Bo( Ay")? + B3A0% + 0(0°). (7.26) 


Here p is no longer the mean pressure 2.12. For an incompressible New- 
tonian fluid 6. = 8; = 0, and f; is called the viscosity. For an incompressible 
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Reiner-Rivlin fluid 8; = 0; thus we again see that already in the first cor- 
rection to the Newtonian fluid, the Reiner-Rivlin fluid yields a result less 
general than that obtained for simple fluids. 

The material constants 6, , By , and 8; are related by the following formu- 
las to the three material functions 7, o; , and o, defined in Section 3: 


x>0 K 
Cor nen e 
Be - = cog z ah 2 25 (7.27) 
2 2 
2h oi(k) — oo(x) it t doy(x) S d o2(k) 
ee oa aes ldo diam (é 


In the second-order approximation 7.26 the “shear dependence of the 
viscosity” does not yet appear; however, all the ‘normal stress effects” 
are already present. 

Motivated by 7.26 we define an incompressible second-order fluid by 
the constitutive equation 


T = S+ pl = 2nD + 48,D" + BAs. (7.28)* 


Of course, an incompressible Newtonian fluid is a special case of 7.28 with 
62 = B; = 0. In many important dynamical situations} the constitutive 
Equation 7.28 leads to linear field equations just as in the Newtonian case. 
For this reason many of the problems that can be solved for Newtonian 
fluids can also be solved for these second-order fluids. 


Comparison with Dimensional Arguments 
If the kinematics are such that C,(t — r) = C(r) is an analytic function 


of r for0 <r < », then this function is determined by its derivatives at 
r = 0, which by 1.12 are essentially the Rivlin-Ericksen tensors A,. Thus 


_ the restriction of the functional $7» [C(r)] to these special analytic ‘/his- 
_ tories may be regarded as an isotropic function of all the Rivlin-Ericksen 


tensors. In those special cases in which only a finite number m of the 


_ tensors A, are important{ we are led to a relation of the form 


S(t) = G(Ai, ++: Am), (7.29) 


_ where G is an isotropic function of m tensor variables. Equation 7.29 is 
_ the constitutive equation of Rivlin-Ericksen fluids.“ By starting from 
_ the Relation 7.29 and using dimensional analysis one can obtain approxi- 


* Langlois!® used an equation of the form 7.28 to represent ‘‘slightly viscoelastic 


_ fluids.”’ We believe our derivation of 7.26 gives 7.28 a sound physical interpretation. 


- For example, unsteady shearing flows. ; 
thor Scaanpla, in the flow problems considered in Sections 3 and 4 all the A, 


- vanish except A; and A2. 
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mation formulas that are formally the same as those exhibited above. 
However, such a formal procedure does not elucidate the physical signifi- 
cance behind our equations; that is, it does not indicate the conditions un- 
der which an experimenter can expect them to be useful. The derivation we 
have outlined here assumes for C(r) only a finite norm and the existence 
of those derivatives at r = 0 that occur in the approximation formulas. 
For r ~ 0, C(r) need not even be continuous, much less analytic. Fur- 
thermore, our treatment shows that the crucial point is not smoothness of 
the history C(r) for all r but rather smoothness of the functional in a 
function space 8, whose norm accounts explicitly for “imperfect memory.” 
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